Displaying 1 result from an estimated 1 matches for "hyptabl".
Did you mean:
hyptable
2013 Jul 20
1
BH correction with p.adjust
...s the original p value, n is the number of tests and i is the rank of the p value. However when I try and recreate the corrected p from my most significant value it does not match up to the one computed by the method p.adjust:
> setwd("C:/work/Methylation/IMA/GM/siteLists")
>
> hypTable <- read.delim("hypernormal vs others.txt")
> pList <- hypTable$p
> names(pList) <- hypTable$site
>
> adjusted <- p.adjust(pList, method="BH")
> adjusted[1]
cg27433479
0.05030589
>
> pList[1]*nrow(hypTable)/1
cg27433479
0.09269194
I tried...