Dear all, I have solve using proc genmod to perform modified poisson regression to generate Risk ratio,
here, how can I adjusted the site effect, would the following program appropriate?
proc genmod data = stroke;
class enroll_order group(desc) score(desc) IVT sites/PARAM=GLM;
model score = group ivt sites/ dist=poisson link=log covb corrb;
repeated subject = enroll_order/type = un;
estimate "beta" group 1 -1/exp;
lsmeans group/diff exp cl;
run;
Besides, how can I use mixed model to get risk ratio with sites as random effect?
I'm confuse to get the RR ratio, and the class group(desc) seems to wrong in proc glimmix,
How can I calculate RR with robust estimator as Zou.et(2004)
proc glimmix data=stroke EMPIRICAL; class site group; model Score= group/ dist=poisson link=log SOLUTION; random site/ subject= enroll_order; run;
... View more