[BioC] Differential expression testing for groups with unequal variances/dispersions?

Ryan C. Thompson rct at thompsonclan.org
Fri May 24 21:10:09 CEST 2013


Hi all,

I am studying a ChIP-Seq dataset (looking at gene promoter regions in 
human) where it appears that different experimental groups have widely 
different dispersions/variances using edgeR/limma. I have 4 timepoints, 
and if I use edgeR to compute the dispersion for each timepoint 
separately, I get:

0 hours: 0.407
24 hours: 0.505
120 hours: 0.115
2 weeks: 0.0531

So the dispersion seems to range from 0.05 to 0.5. I am looking to test 
for "differential modification" between these timepoints, as well as 
between cell types at each timepoint, etc., and I was wondering if there 
is any differential expression test (or dispersion estimation method?) 
that can handle groups with different dispersions/variances.

For reference, here is my experimenal design as an Excel spreadsheet: 
https://www.dropbox.com/s/3vnk4mai3dh39yv/chipseq-samples.xlsx

And here is the result of plotBCV on each group (look at the last 4 
pages for the time point groups): 
https://www.dropbox.com/s/s4caq1p0h3e4zhm/groupdisps.pdf (Warning: big 
PDF with lots of points which may bring your PDF reader to its knees.)

-Ryan Thompson



More information about the Bioconductor mailing list