[R] Fwd:
aoife doherty
aoife.m.doherty at gmail.com
Tue Dec 22 14:29:34 CET 2015
Dear all,
I am hoping to use the mt.maxT package, so I'm trying to understand how the
maxT algorithm works in the multtest package.
I have an example of data:
SNP1 p value= 0.02
SNP1 permutation p values = 0.03, 0.03, 0.03, 0.04, 0.04.
SNP2 p value =0.03
SNP2 permutation p values = 0.4,0.5,0.5,0.6,0.7
SNP3 p value = 0.3
SNP3 p value permutation p values = 0.8,0.9,0.9,0.9,0.9
I was told that for this data set, the corrected P value for SNP1 is 1/6,
for SNP2 is 5/12 and SNP3 is 1. I am so confused as to how these numbers
were reached.
I've tried to come up with an explanation in my head, I would appreciate if
someone could tell me where I'm going wrong;
1.
First, you get the smallest permuted P Value across all of the SNPs for
each permutation. In this case, it's: Perm1 smallest P val= 0.03 Perm2 =
0.03 Perm3 = 0.03 Perm4 = 0.04 Perm5 = 0.04
2.
Then I was told that the corrected p value for SNP1 is 1/6. To get the
6, is this the 5 permuted p values + my p value? So it's 1/6 chance of
seeing a p value as low as mine in the original p value + permuted p value
set.
3.
For SNP2 with an uncorrected p value of 0.03, I am totally confused as
to how the answer is 5/12. I know ties count as 0.5. So in my smallest
permuted data set: 0.03, 0.03, 0.03, 0.04, 0.04: if you add them up
(allowing 0.03 and 0.04 to be worth 0.5 each since they are both ties), the
sum of the 5 permuted values is 2.5. So then the chance of seeing 0.03,
using the logic from step 1, is that it's 2.5/5, or 3.5/6 if you add in the
uncorrected p value? Which isn't 5/12, but if you don't add in the original
p value 0.03, the answer will be 2.5/6, which is 5/12....but I don't
understand why not to add in the original p value in this step, when I did
it in step 1.
I'm obviously not understanding something, if someone could really simply
explain the calculation process/algorithm for maxT corrected p values for
this example I would appreciate it.
Thanks
[[alternative HTML version deleted]]
More information about the R-help
mailing list