Computing the plain P-Value for a given lower-triangular array of genotype counts

Share:

Description

Computes the plain P-Value associated to the provided lower-triangular array of genotype counts to be consistent with the Hardy-Weinberg equilibrium model

Usage

1
HW.plain(model_dist, rms, chisq, gsq, T, n)

Arguments

model_dist

Model distribution for the input genotype count matrix determined in HW.pval(). Calculated via the function create.model()

rms

Root Mean Square test statistic determined in HW.pval()

chisq

Chi-Square test statistic determined in HW.pval()

gsq

Log Likelihood-Ratio test statistic determined in HW.pval()

T

Number of Monte-Carlo simulations desired.

n

Total count of observed genotypes

Details

Determines the plain P-value via Monte-Carlo simulation as described in Algorithm 5.1 of the referenced paper

Returns plain P-values associated to the root-mean-square, chi-square, and log likelihood-ratio statistics.

Author(s)

Shubhodeep Mukherji <deep.mukherji@utexas.edu>

References

"Testing Hardy-Weinberg equilibrium with a simple root-mean-square statistic" by Rachel Ward.

See Also

HW.pval, create.model, test.rms, test.chisq, and test.gsq