[BioC] edgeR and FDR values always equals 1

A Gossner a.gossner at ed.ac.uk
Thu Nov 11 11:44:07 CET 2010


Hi,

While using edgeR to analysis my Tag-seq data, no matter which way I 
analyse the data common or tagwise dispersion the FDR value is always 1. 
Typical output is shown below;

>  d$prior.n
[1] 10
>  head(d$tagwise.dispersion)
[1] 4.360269 5.192625 5.006097 5.006097 3.960243 5.192625
>  summary(d$tagwise.dispersion)
    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
   3.511   4.665   5.006   4.926   5.193   7.240
>  d$common.dispersion
[1] 4.884378
>  prior.n <- estimateSmoothing(d) prior.n
[1] 3.575805e-05
>
>  de.tagwise <- exactTest(d, common.disp = FALSE)
Comparison of groups:  INF - CON
>  topTags(de.tagwise)
Comparison of groups:  INF - CON
                         logConc      logFC       PValue FDR
CATGGGAACAATAAACTCCAC -17.88755 -10.968487 0.0004244306   1
CATGTCTGCCCAAGCACCTAC -32.20325 -35.625605 0.0009657247   1
CATGTTCCTGGTAGCACAAAT -18.89883  -9.075728 0.0011952265   1
ATTGATGTTCTACACCACATG -32.92827 -34.175568 0.0014322637   1
AAAAGGATGACTTCACTCATG -32.99637 -34.039364 0.0016256396   1
CATGATGTGACTTTTAAGTCC -32.81141 -34.409298 0.0018494476   1
AAACCCAGGGAAAGAAGCATG -32.98827 -34.055559 0.0018645461   1
CTTTTTAGATCAAAAAGCATG -33.17368 -33.684744 0.0022498203   1
CGGTCTTATTTAGGAGACATG -32.83447 -34.363162 0.0026414159   1
CATGCTCTTTATCACACCCCC -33.05201 -33.928086 0.0027765733   1
>  sessionInfo()
R version 2.11.1 (2010-05-31)
x86_64-pc-mingw32

locale:
[1] LC_COLLATE=English_United Kingdom.1252
[2] LC_CTYPE=English_United Kingdom.1252
[3] LC_MONETARY=English_United Kingdom.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United Kingdom.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] edgeR_1.6.15

loaded via a namespace (and not attached):
[1] limma_3.4.5

Would be grateful for any suggestions?

Regards

Anton




-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/bioconductor/attachments/20101111/d7fc1789/attachment.pl>


More information about the Bioconductor mailing list