[R] Fwd:

[aoife doherty]

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]]