I am using the tool RNASeqReadSimulator (http://alumni.cs.ucr.edu/~liw/rnaseqreadsimulator.html), in order to simulate RNA-seq reads. In the fist script (there are 3), it is possible to specify expression values to the genes.
python genexplvprofile.py -h -e/--lognormal mu,sigma Specify the mean and variance of the lognormal distribution used to assign expression levels. Default -4,4
I'm not very good at statistics and I would like to know, which parameters of mu and sigma should be OK if I want to have all the genes expressed in my simulated data.
Any idea or suggestion..?