|
|
[学习资料] handy-weinberg平衡性检验的卡方值计算方法 分享~~
请教：
handy-weinberg平衡性检验的卡方值怎么计算的吗？其中的B项怎么计算的。三个数字算的，有公式吗？
方法一、有强者厉害者已经做出了EXCELL的模板，直接输入三个数字就得出结果了。简单操作，易懂，一定深得临床或者没有学习时间的人才的欢迎。
见附件。
方法二 SAS程序：
基本思路啊：gg gt tt是实际基因型频数，以此算出算出单个等位基因（g和t）的实际频率pg和pt，按照遗传平衡定律，可以分别算出 gg gt tt 的理论基因型频数agg agt att，然后按照拟合优度检验的方法来检验gg gt tt是实际基因型频数和理论基因频数的分布是否存在差异，卡方值为a，再算出算出p值。
data m;
input gg gt tt;
n=gg+gt+tt;
pg=(2*gg+gt)/(2*n);
pt=(2*tt+gt)/(2*n);
agg=pg**2*n;
agt=2*pg*pt*n;
att=pt**2*n;
a=(gg-agg)**2/agg+(gt-agt)**2/agt+(tt-att)**2/att;
p=1-probchi(a,1);
cards;
26 24 6
108 52 21
41 57 7
43 53 30
84 109 54
23 23 2
59 45 13
61 80 26
;
proc print;
run; |
|
发表于 2011-2-28 11:24
| 