[BioC] edgeR and sagenhaft

Naomi Altman naomi at stat.psu.edu
Thu Feb 19 14:08:01 CET 2009


A1                 B1               A2             B2
1376421      1577198      1948700      2448499



At 11:48 PM 2/18/2009, Mark Robinson wrote:

>What are the library sizes?
>
>
>On 19/02/2009, at 3:06 PM, Naomi Altman wrote:
>
>>Here is some of the discrepant output.   A1, A2, B1, B2 are the 4
>>samples with tag counts.
>>
>>My commands were
>>y=cbind(A1,B1,A2,B2)
>>d<- DGEList(data = y, group = c(1,2,1,2), lib.size = apply(y,2,sum))
>>alpha <- alpha.approxeb(d)
>>ms <- deDGE(d, alpha = alpha$alpha)
>>
>>Here are the first 5 cases where both samples had Fisher's exact test
>>p-values <.001 and the edgeR exact statistic >.05.  These
>>are raw p-values.
>>
>>
>>          A1            B1          A2         B2
>>Sagenhaft p-value A1 vs B1       Sagenhaft p-value A1 vs B1        ms $exact
>>          46           5           39          12
>>  7.546362e-11                          3.605042e-06
>>       0.76524756
>>          33           4           45          13
>>  1.128468e-07                          3.791248e-07
>>       0.97325389
>>           0          55           13          49
>>  1.070343e-15                          1.668152e-04
>>       0.08834922
>>        69         179           92         526
>>2.088333e-09                          3.816029e-55
>>     0.56544378
>>     109          36          122          63
>>3.988028e-12                          4.556175e-09
>>     0.86574462
>>
>>
>>
>>Naomi S. Altman                                814-865-3791 (voice)
>>Associate Professor
>>Dept. of Statistics                              814-863-7114 (fax)
>>Penn State University                         814-865-1348
>>(Statistics)
>>University Park, PA 16802-2111
>>
>>         [[alternative HTML version deleted]]
>>
>>_______________________________________________
>>Bioconductor mailing list
>>Bioconductor at stat.math.ethz.ch
>>https://stat.ethz.ch/mailman/listinfo/bioconductor
>>Search the archives: 
>>http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>------------------------------
>Mark Robinson
>Epigenetics Laboratory, Garvan
>Bioinformatics Division, WEHI
>e: m.robinson at garvan.org.au
>e: mrobinson at wehi.edu.au
>p: +61 (0)3 9345 2628
>f: +61 (0)3 9347 0852
>------------------------------
>
>
>
>

Naomi S. Altman                                814-865-3791 (voice)
Associate Professor
Dept. of Statistics                              814-863-7114 (fax)
Penn State University                         814-865-1348 (Statistics)
University Park, PA 16802-2111



More information about the Bioconductor mailing list