pith. sign in

arxiv: 1603.02757 · v2 · pith:FVA6LZJRnew · submitted 2016-03-09 · 🧮 math.ST · math.CO· stat.TH

Permutation p-value approximation via generalized Stolarsky invariance

classification 🧮 math.ST math.COstat.TH
keywords valuesapproximationinvarianceestimatepointprinciplestolarskycommon
0
0 comments X
read the original abstract

It is common for genomic data analysis to use $p$-values from a large number of permutation tests. The multiplicity of tests may require very tiny $p$-values in order to reject any null hypotheses and the common practice of using randomly sampled permutations then becomes very expensive. We propose an inexpensive approximation to $p$-values for two sample linear test statistics, derived from Stolarsky's invariance principle. The method creates a geometrically derived set of approximate $p$-values for each hypothesis. The average of that set is used as a point estimate $\hat p$ and our generalization of the invariance principle allows us to compute the variance of the $p$-values in that set. We find that in cases where the point estimate is small the variance is a modest multiple of the square of the point estimate, yielding a relative error property similar to that of saddlepoint approximations. On a Parkinson's disease data set, the new approximation is faster and more accurate than the saddlepoint approximation. We also obtain a simple probabilistic explanation of Stolarsky's invariance principle.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.