用R包genetics做hardy weinberg检测

1、打开R，安装genetics包。依次输入下列内容install("genetics")library(genetics)

2、输入数据，比如一对等位基因为A和G，在case组中基因型为AA、AG和G蘅荫酸圉G的人数分别为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)