Error in as.matrix(basemeans)[, group_id] : subscript out of bounds (simulate_experiment)

[mw55309]

I basically want to simulate one group of reads:

simulate_experiment(my.fa, reads_per_transcript=nreads, numreps=1, outdir="test_sim") Error in stopifnot(is.numeric(fold_changes)) : argument "fold_changes" is missing, with no default simulate_experiment(my.fa, reads_per_transcript=nreads, numreps=1, outdir="test_sim", fold_changes=1) Error in as.matrix(basemeans)[, group_id] : subscript out of bounds

my.fa [1] "test.fa" nreads [1] 1935

sessionInfo() R version 3.2.2 (2015-08-14) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Scientific Linux 7.2 (Nitrogen)

locale: [1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_GB.UTF-8 LC_COLLATE=en_GB.UTF-8 [5] LC_MONETARY=en_GB.UTF-8 LC_MESSAGES=en_GB.UTF-8 [7] LC_PAPER=en_GB.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C

attached base packages: [1] stats4 parallel stats graphics grDevices utils datasets [8] methods base

other attached packages: [1] polyester_1.9.5 Biostrings_2.38.4 XVector_0.10.0 [4] IRanges_2.4.8 S4Vectors_0.8.11 BiocGenerics_0.16.1

loaded via a namespace (and not attached): [1] zlibbioc_1.16.0 limma_3.26.9 tools_3.2.2 logspline_2.1.9

[mw55309]

Actually I just got num_reps wrong....

[jtleek]

Ok cool!

On Thu, Aug 18, 2016 at 8:20 AM mw55309 notifications@github.com wrote:

Actually I just got num_reps wrong....

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/alyssafrazee/polyester/issues/35#issuecomment-240706402, or mute the thread https://github.com/notifications/unsubscribe-auth/ABf7WqEHjPpfEI0PbAuy0AnbakfSZnqNks5qhE39gaJpZM4JnbSq .