请选择 进入手机版 | 继续访问电脑版
楼主: TFBHG24865
2199 2

怎么用sas作蒙特卡洛求列联表的精确p值 [推广有奖]

  • 6关注
  • 0粉丝

本科生

20%

还不是VIP/贵宾

-

威望
0
论坛币
9 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
773 点
帖子
52
精华
0
在线时间
114 小时
注册时间
2013-10-24
最后登录
2018-5-5

TFBHG24865 发表于 2014-6-22 20:23:04 |显示全部楼层 |坛友微信交流群
30论坛币

最佳答案

jingju11 查看完整内容

不断的重复。我猜想可能如下所示。 京剧 http://blog.sina.com.cn/s/blog_a3a926360101kq9r.html
关键词:蒙特卡洛 列联表 蒙特卡 蒙特卡洛
jingju11 发表于 2014-6-22 20:23:05 |显示全部楼层 |坛友微信交流群
不断的重复。我猜想可能如下所示。
京剧

http://blog.sina.com.cn/s/blog_a3a926360101kq9r.html

使用道具

cos-kaka 发表于 2014-6-24 18:48:28 |显示全部楼层 |坛友微信交流群
data Y;
input a b c d;
XPL_FISH = cdf('hyper', a, a+b+c+d, a+b, a+c);
XPR_FISH = 1-cdf('hyper', a-1, a+b+c+d, a+b, a+c);
if a/(a+b)<c/(c+d) then do;
        p = pdf('hyper', a, a+b+c+d, a+b, a+c);
        do e=c to 0 by -1;
                p1 = pdf('hyper', e, a+b+c+d, c+d, a+c);
                if p1<p then leave;
        end;
        XP2_FISH = XPL_FISH + cdf('hyper', e, a+b+c+d, c+d, a+c);
end;
else do;
        p = pdf('hyper', c, a+b+c+d, c+d, a+c);
        do e=a to 0 by -1;
                p1 = pdf('hyper', e, a+b+c+d, a+b, a+c);
                if p1<p then leave;
        end;
        XP2_FISH = XPR_FISH + cdf('hyper', e, a+b+c+d, a+b, a+c);
end;
cards;
53 430 15 176
;
run;

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 20:08