SAS如何计算泊松分布的置信区间和p值

1316

收藏 2024-10-20

我有一个队列，恶性肿瘤预期发病0.1323人，实际发病4人，SIR(标准化发病率比)=4/0.1323=30.2。假定恶性肿瘤发病满足泊松分布，我想用SAS软件进行零假设为SIR=1，α=0.05的假设检验，计算p值和95%置信区间，请问用SAS应该如何计算

我看别人的文献里，计算出来的置信区间是（8.2，77.4），应该不是用的近似法，应该是用精确法（exact method）计算的，但我怎么都复现不出别人的p值和置信区间来，请指教，在SAS里应该怎么计算泊松分布的置信区间

扫码加我拉你入群

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

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

全部回复

extension

2024-10-21 13:18:46

data _null_;
n=4;/*发病例数*/;
eir=0.1323;/*预期发病率*/;
sir=30.2;/*标准化发病率比*/;
alpha=0.05;/*一类错误*/;
lcl_exact=round(cinv(alpha/2,2*n)/2/eir,0.1);/*精确法区间下限*/;
ucl_exact=round(cinv(1-alpha/2,2*(n+1))/2/eir,0.1);/*精确法区间上限*/;
put _all_;
run;

扫码加我拉你入群

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

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

drxx

2024-10-23 16:52:39

extension 发表于 2024-10-21 13:18
data _null_;
n=4;/*发病例数*/;
eir=0.1323;/*预期发病率*/;

谢谢大佬，确实复现出来了

太感谢了！！！跟R语言算出来的也是一样的。另外我还想请教一下，还是这个例子里面，p值用SAS该怎么计算呢

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群