HWtest

Test if a locus is under Hardy Weinberg Proportion

Description

The Hardy Weinberg equilibrium is a fundamental law in genetic.
This function was deeply rewrited.
If a locus is biallelic the function use the exact Hardy Weinberg test (similar to Fisher exact test) analyzing all possible tables. If you download TERNPLOT functions (ID:2299) the function plots a De Finetti's Diagram.
Often a locus shows multiple alleles with low frequency, so the application of chi square is avoided. In this case the function use a Monte Carlo algorithm to perform the exact test.
In the zip file you can find a demo.

You can visit my homepage http://home.tele2.it/cardillo
My profile on XING http://www.xing.com/go/invita/13675097
My profile on LinkedIN http://it.linkedin.com/in/giuseppecardillo

Acknowledgements

The author wishes to acknowledge the following in the creation of this submission:
Fisherextest

• MATLAB release: MATLAB 7.3 (R2006b)

[Dmitriy Fradkin]

This package has a nice demo, implements useful functionality, however it cannot be used inside other programs without further modifications to hwetest.
I would suggest adding a verbose flag to turn off all the printing, and making hwetest return some values.