全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
6184 5
2013-10-18
悬赏 888 个论坛币 已解决
data a;
input a b c d;
cards;
12 132 9 78
;
run;
/*a为实验组的阳性数,b为实验组的阴性数
c为对照组的阳性数,D为对照组的阴性数*
如何生成两个新的变量p1 p2
p1为卡方检验的p值
p2为fisher精确检验的p值
不通过proc freq过程
怎么利用数学运算以及SAS函数完成?*/

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-10-18 09:23:03
复制代码
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 10:39:48
按照公式,卡方检验的P值应该可以这样求:
data a;
input a b c d;

*** Chi Square;
sum=sum(a,b,c,d);
chi_square=((a*d-c*b)**2*sum)/((a+b)*(c+d)*(a+c)*(b+d));
p_chi=1-probchi(chi_square,1);
put chi_square= p_chi=;

cards;
12 132 9 78
;
run;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 10:54:02
pobel 发表于 2013-10-18 10:39
按照公式,卡方检验的P值应该可以这样求:
data a;
input a b c d;
非常感谢,关键是fisher精确检验的p值不知道怎么算
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 10:56:04
Imasasor 发表于 2013-10-18 10:54
非常感谢,关键是fisher精确检验的p值不知道怎么算
Fisher的我也不知道。版主可以问一下京剧大侠。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 13:45:29
按照fisher检验最原始的公式,我也贴出一个。因为楼主的数据太大,不符合fisher检验的条件,然后在sas运行由于阶乘太大报错,所以我就改小了,使其总和<40. 顺便说一句楼上大神啊。
复制代码
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群