[BioC] Limma and negative p-values
Ovokeraye Achinike-Oduaran
ovokeraye at gmail.com
Thu Mar 15 14:50:32 CET 2012
Hi all,
Sorry...a recall of my previous mail...the p-values are fine(or not)
but not negative...the corrected p-values, however, look very shady.
Any suggestions there?
Thanks
On Thu, Mar 15, 2012 at 3:46 PM, Ovokeraye Achinike-Oduaran
<ovokeraye at gmail.com> wrote:
> Hi all,
>
> Trying to use limma(code below) to do some DE analysis...and I get
> negative p-values. Not quite sure where I went wrong here...please
> help!
>
> Thanks.
>
> Avoks.
>
> gds158dat = getGEO('GDS158',destdir=".")
> gds158eset = GDS2eSet(gds158dat, do.log2=TRUE)
> groups= pData(gds158eset)$metabolism
> groups=as.character(groups)
> groups[groups=="insulin sensitive"]= "IS"
> groups[groups=="insulin resistant"]= "IR"
> f = factor(groups, levels=c("IS","IR"))
> design_gds158 = model.matrix(~0+f)
> colnames(design_gds158) = levels(f)
> cont.matrix = makeContrasts(IR-IS, levels=design_gds158)
> fit=lmFit(gds158eset, design_gds158)
> fit2 = contrasts.fit(fit, cont.matrix)
> fit2 = eBayes(fit2)
> results = topTable(fit2, adjust ="BH", number = nrow(gds158eset))
>
> sessionInfo()
> R version 2.14.1 (2011-12-22)
> Platform: i386-pc-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=English_xxx.1252 LC_CTYPE=English_xxx.1252
> [3] LC_MONETARY=xxx.1252 LC_NUMERIC=C
> [5] LC_TIME=xxx.1252
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] ArrayExpress_1.14.0 GEOquery_2.20.8 limma_3.10.2
> [4] XML_3.9-4.1 RCurl_1.91-1.1 bitops_1.0-4.1
> [7] puma_2.6.0 mclust_3.4.11 affy_1.32.1
> [10] Biobase_2.14.0
>
> loaded via a namespace (and not attached):
> [1] affyio_1.22.0 BiocInstaller_1.2.1 preprocessCore_1.16.0
> [4] tools_2.14.1 zlibbioc_1.0.0
More information about the Bioconductor
mailing list