exact.pval1s: p-value correction to the one-sided version of exact NNCT...

In nnspat: Nearest Neighbor Methods for Spatial Patterns

p-value correction to the one-sided version of exact NNCT test

Description

In using Fisher's exact test on the 2 \times 2 nearest neighbor contingency tables (NNCTs) a correction may be needed for the p-value. For the one-sided alternatives, the probabilities of more extreme tables are summed up, including or excluding the probability of the table itself (or some middle way). Let the probability of the contingency table itself be p_t=f(n_{11}|n_1,n_2,c_1;θ_0) where θ_0=(n_1-1)(n_2-1)/(n_1 n_2) which is the odds ratio under RL or CSR independence and f is the probability mass function of the hypergeometric distribution. For testing the one-sided alternative H_o:\,θ=θ_0 versus H_a:\,θ>θ_0, we consider the following four methods in calculating the p-value:

• [(i)] with S=\{t:\,t ≥q n_{11}\}, we get the table-inclusive version which is denoted as p^>_{inc},

• [(ii)] with S=\{t:\,t> n_{11}\}, we get the table-exclusive version, denoted as p^>_{exc}.

• [(iii)] Using p=p^>_{exc}+p_t/2, we get the mid-p version, denoted as p^>_{mid}.

• [(iv)] We can also use Tocher corrected version which is denoted as p^>_{Toc} (see tocher.cor for details).

See (\insertCiteceyhan:exact-NNCT;textualnnspat) for more details.

Usage

exact.pval1s(ptable, pval, type = "inc")

Arguments

ptable	Probability of the observed 2 \times 2 NNCT under the null hypothesis using the hypergeometric distribution for Fisher's exact test.
pval	Table inclusive p-value for Fisher's exact test on the NNCT.
type	The type of the p-value correction for the one-sided exact test on the NNCT, default="inc". Takes on values "inc", "exc", "mid", "tocher" (or equivalently 1-4, respectively) for table inclusive, table-exclusive, mid-p-value, and Tocher corrected p-value, respectively.

Value

A modified p-value based on the correction specified in type.

Author(s)

Elvan Ceyhan

References

\insertAllCited

See Also

exact.pval2s and tocher.cor

Examples

ct<-matrix(sample(20:40,4),ncol=2)
ptab<-prob.nnct(ct)
pv<-.3
exact.pval1s(ptab,pv)
exact.pval1s(ptab,pv,type="exc")
exact.pval1s(ptab,pv,type="mid")

nnspat documentation built on Aug. 30, 2022, 9:06 a.m.