topic Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX in Statistical Procedures

How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

Diels_O — Tue, 06 Sep 2022 13:28:33 GMT

Dear all, I met a problem in the PROC GLIMMIX,

How can I estimate the RR ratio using PROC GLIMMIX with dist = poisson and link = log?

proc glimmix data=stroke EMPIRICAL;
class group enroll_order;
model Score = group / dist=poisson link=log SOLUTION;
random site/ subject= enroll_order;
run;

Also, when adding the EMPIRICAL, my computer memory was not sufficient,

I will be apropriate if I can get you help! Sincerely, Chagnwei.

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

SteveDenham — Tue, 06 Sep 2022 13:40:54 GMT

In this case, you can probably get estimates if you specify 'site' in the CLASS statement. The more I look at this, I really need an explanation of enroll_order and how it is coded. I don't think your RANDOM statement correctly models the design in your data, as I don't see how to include 'enroll_order' as a subject for each site. Without more information, I suspect that there is a different enroll_order for each site, which would mean that the two variables are confounded. Please include more information regarding the data structure, and also include any messages in the log that may be generated.

SteveDenham

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

Diels_O — Tue, 06 Sep 2022 16:07:26 GMT

proc glimmix data=stroke EMPIRICAL;
   class group site enroll_order;
   model Score = group baseline_NIHSS baseline_ASPECTS onset_to_randomization / dist=binomia link=log SOLUTION;
   random site/ subject= enroll_order;
run;

the group refers to trt;

the enroll_order refers to ID of each individuals;

Also I have change the dist from dist = poisson to dist = binomial;

How can I get the RR ratio?

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

Diels_O — Wed, 07 Sep 2022 09:06:53 GMT

Dear Steve, the enroll_order was not related to the sites;

beside, the Score was a binary variable

, and the reasult

proc glimmix data=stroke EMPIRICAL;
   class site trt   enroll_order;
   model score = trt / dist=poisson link=log SOLUTION;
   store glim;
run;
proc plm restore = glim;
	lsmeans trt / ilink diff exp cl;
run;

It works.

But when try to add random effect, it shows insufficient memory

proc glimmix data=stroke EMPIRICAL;
   class site trt   enroll_order;
   model score = trt / dist=poisson link=log SOLUTION;
   random site/subject = enroll_order;
run;

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

SteveDenham — Fri, 09 Sep 2022 13:05:18 GMT

If you apply the %NLmeans macro to the results from using a STORE statement (see this note https://support.sas.com/kb/62/362.html ), you should get the RR.

Also of interest is this note https://support.sas.com/kb/23/003.html , which includes the Zou et al. method you mentioned earlier in one of these threads.

SteveDenham

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

Diels_O — Sun, 11 Sep 2022 06:07:15 GMT

Many thanks to you! Thank you Steve!

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

SteveDenham — Tue, 13 Sep 2022 12:30:32 GMT

I'll make two suggestions, but won't guarantee either will help.

First, since enroll_order is sequential, and there are 948 subjects, you might make things easier if you remove enroll_order from the CLASS statement.

Second, you have a multilevel design, where enroll_order characterizes the residual error, and site is a random effect. To assist convergence, you might try:

proc glimmix data=stroke EMPIRICAL;
   class site trt   ;
   model score = trt / dist=poisson link=log SOLUTION;
   random intercept/subject = site;
run;

Hope this helps.

SteveDenham

Re: How to calculate RR with quasi-Poisson distribution using PROC GLIMMIX

Diels_O — Wed, 14 Sep 2022 02:12:49 GMT

Thank you again Steve!

By reading other books, I write the following code:

ods trace on;
ods output diffs = out1;
proc glimmix data=stroke;
   class site trt   enroll_order;
   model score = trt/ dist=binomia link=log SOLUTION;
random intercept trt/subject = site;
lsmeans trt/diff CL;
run;

ods trace off;
proc print data = out1;run;

data exp1;
set out1;
rr = exp(estimate);
lowcl = exp(lower);
uppercl = exp(upper);
run;

Just one question, is the "random intercept trt/ subject = site" equal to "random intercept/subject = site"? This two run out the same result.

Sincerely, Diels