点估计95%置信区间

shiyiming

事件发生数服从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

[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]