I am currently trying to simulate one sample and I am using the reads_per_transcript as a vector, however, I noticed that the number of reads simulated is not equal to the number of reads in reads_per_transcript (ie. using c(100, 50) does result in 100 reads for first transcript and 50 reads from second transcript). It seems like you take a distribution and the reads_per_transcript is just the mean number of transcript.
I was wondering if it would be possible to have it so that the reads simulated is exactly equal to the ones specified in reads_per_transcript vector.
Hi polyester, this is a great project!
I am currently trying to simulate one sample and I am using the reads_per_transcript as a vector, however, I noticed that the number of reads simulated is not equal to the number of reads in reads_per_transcript (ie. using c(100, 50) does result in 100 reads for first transcript and 50 reads from second transcript). It seems like you take a distribution and the reads_per_transcript is just the mean number of transcript.
I was wondering if it would be possible to have it so that the reads simulated is exactly equal to the ones specified in reads_per_transcript vector.
Thank you!