topic Re: Deciphering a Note in the SAS Log When Running PROC GLIMMIX with poisson distribution in Statistical Procedures

Deciphering a Note in the SAS Log When Running PROC GLIMMIX with poisson distribution

fdehkord — Wed, 13 Apr 2022 23:27:16 GMT

I have three observations per patients and the response variable is the number of days in ICU. There are 41 patients. The explanatory variables are age (continuous) gender (F, M), and varA (continuous). When I run the following codes. There will be no convergence

proc glimmix data= mydata;

class id Gender ;

model Length_of_ICU_stay= gender Age varA / dist=poi link=log

solution;

random int / subject=id;

random residual;

output out=gmxout pearson=pearson;

run;

When I run this code, I get the message “An R-side variance component is confounded with the profiled variance.

proc glimmix data= mydata;

class id Gender ;

model Length_of_ICU_stay= gender Age varA / dist=poi link=log

solution;

random int / subject=id residual;

output out=gmxout pearson=pearson;

run;

I need help to understand the difference between the two models fitted through the above codes and make sense or the message.

Thanks.

Re: Deciphering a Note in the SAS Log When Running PROC GLIMMIX with poisson distribution

JosvanderVelden — Thu, 14 Apr 2022 10:09:15 GMT

Have seen these threads: https://communities.sas.com/t5/Statistical-Procedures/Need-Help-Deciphering-a-Mysterious-Note-in-the-SAS-Log-When/td-p/568622 & https://communities.sas.com/t5/Statistical-Procedures/GEE-Variance-Component-An-R-side-variance-component-is/m-p/801878

Best regards, Jos

Re: Deciphering a Note in the SAS Log When Running PROC GLIMMIX with poisson distribution

SteveDenham — Thu, 14 Apr 2022 12:12:19 GMT

Recall that for a Poisson distribution the mean equals the variance. By specifying RANDOM residual without indicating a subject you have confounded the R side variance with the residual variance. With 3 visits per subject, perhaps this code will help (no guarantee, though):

proc glimmix data=  mydata;
   class id  Gender visit_no ;
nloptions maxiter=2000 tech=nrridg;
   model Length_of_ICU_stay= gender Age varA visit_no/ dist=poi link=log
         solution;
     random visit_no/type=cs residual subject=id;
   output out=gmxout pearson=pearson;

run

/*  Or this, for a G side approach  */
proc glimmix data=  mydata method=laplace;
   class id  Gender visit_no ;
  nloptions maxiter=2000 tech=nrridg;

 model Length_of_ICU_stay= gender Age varA / dist=poi link=log
         solution;
     random visit_no/type=cs  subject=id;
   output out=gmxout pearson=pearson;

run;;

I chose a compound symmetry variance structure as I doubt that the visits are equally spaced in time, even if indexed as 1, 2 and 3. A heterogeneous compound symmetry (csh) may give a superior fit, but simpler is usually better and has a higher chance of converging. I also added an NLOPTIONS statement that could result in convergence, if you are hitting up against the default 20 iteration limit.

SteveDenham.

Re: Deciphering a Note in the SAS Log When Running PROC GLIMMIX with poisson distribution

fdehkord — Fri, 15 Apr 2022 22:19:57 GMT

Thank you very much Steve. I appreciate your detailed explanation. It was very helpful.