怎么不利用proc freq过程在data步计算fisher检验的p值

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函数完成？*/

1. data a;
   
2. input a b c d;
   
3. XPL_FISH = cdf('hyper', a, a+b+c+d, a+b, a+c);
   
4. XPR_FISH = 1-cdf('hyper', a-1, a+b+c+d, a+b, a+c);
   
5. if a/(a+b)<c/(c+d) then do;
   
6.         p = pdf('hyper', a, a+b+c+d, a+b, a+c);
   
7.         do e=c to 0 by -1;
   
8.                 p1 = pdf('hyper', e, a+b+c+d, c+d, a+c);
   
9.                 if p1<p then leave;
   
10.         end;
    
11.         XP2_FISH = XPL_FISH + cdf('hyper', e, a+b+c+d, c+d, a+c);
    
12. end;
    
13. else do;
    
14.         p = pdf('hyper', c, a+b+c+d, c+d, a+c);
    
15.         do e=a to 0 by -1;
    
16.                 p1 = pdf('hyper', e, a+b+c+d, a+b, a+c);
    
17.                 if p1<p then leave;
    
18.         end;
    
19.         XP2_FISH = XPR_FISH + cdf('hyper', e, a+b+c+d, a+b, a+c);
    
20. end;
    
21. cards;
    
22. 12 132 9 78
    
23. 11 4 2 6
    
24. ;
    
25. run;

复制代码

按照公式，卡方检验的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;

pobel 发表于 2013-10-18 10:39
按照公式，卡方检验的P值应该可以这样求：
data a;
input a b c d;

非常感谢，关键是fisher精确检验的p值不知道怎么算

Imasasor 发表于 2013-10-18 10:54
非常感谢，关键是fisher精确检验的p值不知道怎么算

Fisher的我也不知道。版主可以问一下京剧大侠。

按照fisher检验最原始的公式，我也贴出一个。因为楼主的数据太大，不符合fisher检验的条件，然后在sas运行由于阶乘太大报错，所以我就改小了，使其总和<40. 顺便说一句楼上大神啊。

1. data a;
   
2. input a b c d;
   
3. n=sum(a,b,c,d);
   
4. p_ = (fact(a+b)*fact(a+c)*fact(b+d)*fact(c+d))/(fact(a)*fact(b)*fact(c)*fact(d)*fact(n));
   
5. array p[*] p1 - p15;
   
6. do i = 1 to 15;
   
7.    a1 = 22-i;
   
8.    b1 = i+1;
   
9.    c1 = i-1;
   
10.    d1 = 15-i;
    
11.    p[i] = (fact(a1+b1)*fact(a1+c1)*fact(b1+d1)*fact(c1+d1))/(fact(a1)*fact(b1)*fact(c1)*fact(d1)*fact(n));
    
12.    if p[i] > p_ then p[i] = 0;
    
13.    end;
    
14.    p_fisher = sum(of p1 - p15);
    
15. cards;
    
16. 12 11 9 5
    
17. ;
    
18. run;

复制代码