y=f=13 ; n=46021 # input data alpha=0.05 # to get a 95% interval se=sqrt(f) # estimate se # give intervals to y (CL=qnorm(alpha/2,y,se)) # lower CL (CU=qnorm(1-alpha/2,y,se)) # upper CL # give as rates per 100,000 CL*100000/n # lower CL CU*100000/n # upper CL

># give intervals to y > (CL=qnorm(alpha/2,y,se)) # lower CL [1] 5.933249 > (CU=qnorm(1-alpha/2,y,se)) # upper CL [1] 20.06675 > > # give as rates per 100,000 > CL*100000/n # lower CL [1] 12.89248 > CU*100000/n # upper CL [1] 43.60347