[BioC] edgeR topTags output (logFC)

Susan [guest] guest at bioconductor.org
Tue Aug 27 03:28:03 CEST 2013


I'm just wanting to confirm my understanding of the edgeR topTags() output.  I'm using a design matrix with many coefficients (76 genotypes, 2 treatments[sun, shade]).  

1) In the topTags output, I am looking at logFC which represents the estimated or predicted means -- NOT coefficients, correct?  

2) I have selected one parental genotype to be my reference genotype, and I am comfortable interpreting the logFC values for the other genotypes as comparisons to this genotype.  I have selected one treatment as the reference (sun).  When I see the logFC values for 'trtshade', I assume this is also in comparison to the reference genotype in the reference treatment.  

My question is about looking at 'geno2.trtshade' values - are these values in comparison to 'trtshade' or 'geno2' or the reference genotype/reference treatment values?  

For example:  this is calling glmLRT with coef = grep("trt", colnames(Design)) to identify treatment effects

           Row.names       logFC.trtsha    logFC.genoIL1.1.trtsha   logCPM        LR        PValue           FDR
1 Solyc08g078300.2.1    0.9440264           -0.379030370 4.207443 1135.8061 1.450273e-188 1.637938e-184
2 Solyc06g060830.2.1    0.9630343            0.038988710 4.511893 1095.1525 2.539194e-180 1.433883e-176
3 Solyc06g067910.2.1    1.2468599           -0.096499016 4.795789 1040.5813 2.722511e-169 1.024935e-165
4 Solyc06g008580.2.1    1.1500683           -0.451731914 3.300982  696.1677 6.043868e-101  1.706486e-97
5 Solyc07g017600.2.1    0.5026292            0.007208458 5.209938  694.2070 1.451776e-100  3.279271e-97
6 Solyc10g083760.1.1   -0.2877568            0.045952894 7.033394  651.8629  2.232759e-92  4.202796e-89

Thanks so much!


 -- output of sessionInfo(): 

> sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C         LC_TIME=C            LC_COLLATE=C         LC_MONETARY=C        LC_MESSAGES=C        LC_PAPER=C          
 [8] LC_NAME=C            LC_ADDRESS=C         LC_TELEPHONE=C       LC_MEASUREMENT=C     LC_IDENTIFICATION=C 

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

other attached packages:
[1] ggplot2_0.9.3.1 reshape_0.8.4   plyr_1.8        edgeR_3.2.3     limma_3.16.4   

loaded via a namespace (and not attached):
 [1] MASS_7.3-26        RColorBrewer_1.0-5 colorspace_1.2-2   dichromat_2.0-0    digest_0.6.3       grid_3.0.1         gtable_0.1.2       labeling_0.1      
 [9] munsell_0.4        proto_0.3-10       reshape2_1.2.2     scales_0.2.3       stringr_0.6.2      tools_3.0.1       

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list