[BioC] edgeR topTags output (logFC)
Gordon K Smyth
smyth at wehi.EDU.AU
Wed Aug 28 02:05:25 CEST 2013
Dear Susan,
> Date: Mon, 26 Aug 2013 18:28:03 -0700 (PDT)
> From: "Susan [guest]" <guest at bioconductor.org>
> To: bioconductor at r-project.org, smbush at ucdavis.edu
> Subject: [BioC] edgeR topTags output (logFC)
>
>
> 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?
Incorrect. As the abbreviation suggests, logFC gives you
log-fold-changes. logFC are indeed coefficients, not group means.
> 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.
Yes, assuming that trt is what you have called the treatment factor, then
trtshade would the coefficient estimating the log-fold-change between the
shade and sun treatments.
> 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?
Without knowing what model you have fitted, or which coefficients you have
tested for, it is impossible to tell.
Best wishes
Gordon
> 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
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list