标题:
[学习资料]
handy-weinberg平衡性检验的卡方值计算方法 分享~~
[打印本页]
作者:
wwy
时间:
2011-2-28 11:24
标题:
handy-weinberg平衡性检验的卡方值计算方法 分享~~
请教：
handy-weinberg平衡性检验的卡方值怎么计算的吗？其中的B项怎么计算的。三个数字算的，有公式吗？
下载
(64.81 KB)
2011-2-28 11:23
方法一、有强者厉害者已经做出了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;
图片附件:
未命名.jpg
(2011-2-28 11:23, 64.81 KB) / 下载次数 1
http://cdcman.cn/attachment.php?aid=4629&k=6013a59707614217e70577ba2fb11207&t=1719390421&sid=eGxs9G
作者:
wwy
时间:
2011-2-28 11:25
不好意思哈。附件只能上传PDF等格式，excell和word上传不了。
有意者留邮箱吧。
作者:
ajackieb
时间:
2011-9-10 16:39
需要EXCEL格式模板,邮箱:Jackiemail2000@163.com.谢谢!
欢迎光临 疾控家园 (http://cdcman.cn/)
Powered by Discuz! 7.2
图片附件: 未命名.jpg (2011-2-28 11:23, 64.81 KB) / 下载次数 1