Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: st: ANOVA Posthocs
From
"Ploutz-Snyder, Robert (JSC-SK)[USRA]" <[email protected]>
To
"[email protected]" <[email protected]>
Subject
RE: st: ANOVA Posthocs
Date
Fri, 22 Oct 2010 15:52:35 -0500
Yes David, we are going to run FDR through -multproc- on these t-tests as a second step. I should have mentioned this originally...
Al hit on our larger philosophical question as well. But either way I'd like to be able to capture p-values from post-hoc comparisons--whether Scheffe, Tukey HSD, Tukey-Kramer, or otherwise.
Rob
-----Original Message-----
From: [email protected] [mailto:[email protected]] On Behalf Of Airey, David C
Sent: Friday, October 22, 2010 2:52 PM
To: [email protected]
Subject: re: st: ANOVA Posthocs
.
Maybe if you are doing 40*6=240 pairwise comparisons, do you want to consider controlling the false discovery rate? Along these lines, if you are doing oneway ANOVAs, why not use -regress- with the -test- command (http://www.stata.com/support/faqs/stat/test1.html) and -post- to file as you run through the genes and then submit the p-values to one of Roger Newson's packages like -smileplot or -qqvalue-?
> Greetings Statalisters.
>
> I am working on a multiple comparison problem comparing Y by Group with 4 groups, but I have MULTIPLE Y's that I need to consider because my Y's represent fold-increases in gene expression data and I'm working with over 40 genes...
>
> So in this effort, I'd like to run Oneway ANOVAs of Y by Group for each of my genes, capture the Pairwise comparisons among my groups into a separate dataset, and repeat this with a loop for each of the genes of interest.
>
> Ideally I'd like to use the Tukey-Kramer Pairwise user-written command (thanks Phil Enders), however that does not provide p-values in the screen output or any matrix in the return list.
>
> An alternative would be to use the Scheffe posthoc option to -oneway- but I can't figure out how to steal those p-values either?
>
> Any thoughts on how to grab these p-values?
> Much appreciated...
> Rob
>
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/