HWxtest (version 1.0.5)

p.value: Extract just the P value(s) from a Hardy-Weinberg test.

Description

Use the p.value function to return just the P value(s) from the results of a call to hwx.test. If applied to a list of results, it will return a vector or matrix of P values. You can specify the statName as LLR, Prob, U or Chisq. You can also apply p.value to a matrix or vector and it will attempt to use hwx.test to return a P value. However, it's usually preferable to use hwx.test directly.

Usage

p.value(x, statName = NA)

Arguments

x

The result of a call to hwx.test or a list of such results. It can also be the genotype counts in any of the same formats as accepted by hwx.test

statName

can be LLR, Prob, U, or Chisq

Value

• The P value

References

The methods are described by http://dx.doi.org/10.1534/genetics.109.108977{Engels, 2009. Genetics 183:1431}.

Examples

data(HWcases)
testResults <- hwx.test(HWcases)
p.value(testResults)
p.value(testResults, statName="U")