Question: Polyester's simulate_experiment man pages / vignette description
1
gravatar for Michael Love
4 weeks ago by
Michael Love26k
United States
Michael Love26k wrote:

I was walking through polyester with some students yesterday and remembered something which I seem to forget every now and then about the behavior of simulate_experiment. It's summarized in this closed issue:

https://github.com/alyssafrazee/polyester/issues/15

If you give simulate_experiment a transcript and set reads_per_transcript=0 it will turn that 0 to a 1 here:

https://github.com/mikelove/polyester/blob/master/R/NB.R#L18

I always seem to remember this after playing around a bit, so it's never ended up in one of my simulation datasets. I always remove the desired unexpressed transcripts from the FASTA instead of putting 0s, as recommended in that issue, but I noticed that it's not mentioned in the vignette or man pages of polyester. I could see someone accidentally giving simulated_experiment e.g. all the reference transcripts, and putting a 0 for a majority of them. They may be surprised to find they are expressing every transcript with one read.

polyester • 47 views
ADD COMMENTlink written 4 weeks ago by Michael Love26k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 421 users visited in the last hour