HWxtest (version 1.1.9)

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 Engels, 2009. Genetics 183:1431.

Examples

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

# }