I'm interested in simulating RNA-seq data using the rlsim tool. However, rather than randomly draw the expression levels from a mixture of distributions, I wish to control which transcripts are expressed at relative levels. Basically, I have an input distribution of relative abundances, and I want to simulate an RNA-seq dataset where these relative abundances are the underlying (latent) model parameters. Is this possible in rlsim? I wasn't able to find a way to do this by looking at the user manual. If it is possible, how would I accomplish it?
•
written
by
Rob
1 answer
rlsim cannot simulate from relative abundances (proportions) only, but you can easily "convert "the proportions to the absolute levels by multiplying them by a constant which is large enough for your purposes and then construct the input Fasta using those levels. You might have to play around with the multiplier using the -m flag, make sure that the "sampling ratio" reported by rlsim is small (let's say at least as small as 10^-5).
Indeed it is not that difficult. Let me know if you bump into other obstacles.
Log in to answer this question.
-
Command line or python RNA-seq read simulator
written by i.sudbery
-
Multinomial Distributions for large n
written by rtrende •
-
Calculation of tRNA abundances from different cell types
written by nkabo
-
Differential gene abundance analysis with 16s seq data
written by ic23oluk •
-
Assigning Fold Change for Polyester
written by priyankamaripuri •
-
Subtract two RNA-seq datasets
written by Milos Pjanic •
-
simNGS for small RNA-Seq
written by caner •
-
Genomic feature overlaps taking into account underlying sequence feature distributions
written by Sakti
-
How to pick simulation parameters for the RNASeqReadSimulator
written by iraun
-
Tool: Rlsim, A Package For Simulating Rna-Seq Library Preparation With Parameter Estimation
written by Botond Sipos
