用R包genetics做hardy weinberg检测

1、打开R，安装genetics包。

依次输入下列内容

install("genetics")

library(genetics)

2、输入数据，比如一对等位基因为A和G，在case组中基因型为AA、AG和GG的人数分别为100,50和20个。

mydata <- c(rep("A/A",100),rep("A/G",50),rep("G/G",20))

3、输入以下命令，进行hardy weinberg检测

g1 <- genotype(mydata)

HWE.chisq(g1)