[BioC] RNA-seq data simulation

Mark Robinson mark.robinson at imls.uzh.ch
Tue Aug 21 17:00:36 CEST 2012


You could start with rnbinom().  The 'size' argument is 1/dispersion.

Mark

On 21.08.2012, at 16:51, Alpesh Querer wrote:

> I should mention I want to use negative binomial as underlying distribution.
> 
> Thanks.
> Al
> 
> On Tue, Aug 21, 2012 at 9:41 AM, Alpesh Querer <alpeshq at gmail.com> wrote:
> 
>> Hi all,
>> 
>> I want to simulate counts for 1000 genes for 5 RNA-seq samples (with 3
>> biological replicates each), and I want to provide the dispersion
>> parameters myself.
>> Is there a package or code that can help me do this?
>> 
>> Thanks a lot for your help,
>> Al
>> 
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

----------
Prof. Dr. Mark Robinson
Bioinformatics
Institute of Molecular Life Sciences
University of Zurich
Winterthurerstrasse 190
8057 Zurich
Switzerland

v: +41 44 635 4848
f: +41 44 635 6898
e: mark.robinson at imls.uzh.ch
o: Y11-J-16
w: http://tiny.cc/mrobin

----------
http://www.fgcz.ch/Bioconductor2012
http://www.eccb12.org/t5



More information about the Bioconductor mailing list