标题:
[学习资料]
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) / 下载次数 153
http://cdcman.s1.288idc.com/attachment.php?aid=4629&k=265208d5c7233979b76b2b283a10b3d6&t=1736586198&sid=7l9Y5j
作者:
wwy
时间:
2011-2-28 11:25
不好意思哈。附件只能上传PDF等格式,excell和word上传不了。
有意者留邮箱吧。
作者:
ajackieb
时间:
2011-9-10 16:39
需要EXCEL格式模板,邮箱:Jackiemail2000@163.com.谢谢!
欢迎光临 疾控家园 (http://cdcman.s1.288idc.com/)
Powered by Discuz! 7.2