用R包genetics做hardyweinberg检测,自由度为1,显著性水平P小于0.01时,不满足hardyweinberg平衡。
工具/原料
电脑
R
方法/步骤
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、输入以下命令,进行hardyweinberg检测g1<-genotype(mydata)HWE.chisq(g1)