[BioC] Limma and negative p-values
Ovokeraye Achinike-Oduaran
ovokeraye at gmail.com
Thu Mar 15 14:46:06 CET 2012
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