Hi SASers,
I'm observing an irritating problem, using stdrate: I receive a missing value for LowerCL.
The problem occurs only, when I ...:
- calculate DIRECTLY standardized rates
- using the GAMMA option to calculate confidence intervals
- if there are a lot of observed events (from different tests I would say, the problem starts to occur, when there are mor than 110'000 obserbved events)
Because the problem occurs only with a lot of events, I'm not able to provide sample data. The code looks like this:
ods exclude all;
proc stdrate
data=work.ms_akut
refdata=work.resp
method=direct
cl=gamma
stat=rate(mult=1000)
;
population
event=obs_type
total=pop
;
reference
total=pop_resp
;
strata sex age_class ;
by canton
ods output StdRate=output_stdrate;
run;
ods exclude none;
Does anybody know this problem? Any explanation?
Any suggestions are very much appreciated....
Have a nice day,
Reto