limma p-value calculated how?
1
0
Entering edit mode
Kripa R ▴ 180
@kripa-r-4482
Last seen 10.2 years ago
Hi, I'm currently using LIMMA to analyze two groups, wt vs treatment group. LIMMA says gene A is significantly different between the two groups, p=10E-6, but when I do a simple t-test I am unable to get such a significant p-value. For example my normalized values for ProbeID:A may be as followed: wt - 4.8, 5.1, 5.7 treatment - 8.1, 9.8, 7 if I do >t.test(wt, treatment) I don't get a p-value as low as what is observed in the limma top table. Why would this happen? Thanks for the help, -Kripa [[alternative HTML version deleted]]
limma limma • 2.6k views
ADD COMMENT
0
Entering edit mode
@david-westergaard-5119
Last seen 10.2 years ago
Hello Kripa, Limma uses a moderatet t-test. If you are interested in the details, you can read more in Smyth, G. K. (2004). Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology, Vol. 3, No. 1, Article 3. (Copied from the limma usersguide) Best, David 2012/5/22 Kripa R <kripa777 at="" hotmail.com="">: > > Hi, I'm currently using LIMMA to analyze two groups, wt vs treatment group. LIMMA says gene A is significantly different between the two groups, p=10E-6, but when I do a simple t-test I am unable to get such a significant p-value. > > For example my normalized values for ProbeID:A may be as followed: > wt - 4.8, 5.1, 5.7 > treatment - 8.1, 9.8, 7 > > if I do >t.test(wt, treatment) I don't get a p-value as low as what is observed in the limma top table. Why would this happen? > > > > Thanks for the help, > > > -Kripa > > ? ? ? ?[[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD COMMENT

Login before adding your answer.

Traffic: 976 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6