SAS中文论坛

标题: 点估计95%置信区间 [打印本页]

作者: shiyiming 时间: 2010-12-13 14:37
标题: 点估计95%置信区间
事件发生数服从Poisson分布，设 R=事件累积发生密度，I=事件发生数，t=总人时（人×天），根据Poisson分布的正态近似原理有：
研究中止时事件累积发生密度点估计为 R=I/t*100%
相应95%置信区间为[(I-Z(1-a/2)√I)/t , (I+Z(1-a/2)√I)/t ] 。

现在I R t都可以算出来,都是一个数值,请问这种点估计的95%置信区间应该怎么编程?

作者: shiyiming 时间: 2010-12-13 22:20
标题: Re: 点估计95%置信区间
[code:3jcx39vz]data _null_;
a =0.05;
i =10;
t =1000;
if (i >0 and t >0) then do;
q =quantile('normal', 1-a/2)/sqrt(i);
lci =(1-q)/t; uci =(1+q)/t;
end;
put lci=uci=;
run;[/code:3jcx39vz]

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)