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

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

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

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.

Shubhodeep Mukherji <deep.mukherji@utexas.edu>

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

`HW.pval`

, `create.model`

, `test.rms`

, `test.chisq`

, and `test.gsq`

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.