topic Re: Change in behavior of the group= option in the RANDOM statement in PROC GLIMMIX in Statistical Procedures

Change in behavior of the group= option in the RANDOM statement in PROC GLIMMIX

SteveDenham — Fri, 15 Nov 2013 16:26:38 GMT

As a watch out, be aware that the output changes in 9.4/12.3 from 9.3/12.1 for the following (buried inside a macro

proc glimmix data=for_Stats1 ic=pq;

by param ;

nloptions tech=nrridg;

class grp_no studyday anml_nbr ;

model value = grp_no

studyday

grp_no*studyday

/ ddfm=kr2;

random studyday /residual type=&covtype subject= anml_nbr(grp_no) group=grp_no;

%if "&covtype" = "AR(1)" OR "&covtype" = "ARH(1)" %then %do;

random intercept/subject=anml_nbr(grp_no) group=grp_no;

%end;

covtest homogeneity;

ods output fitstatistics = &outdata._a convergencestatus = &outdata.status ;

run;

If &covtype is AR(1) or ARH(1), then under 9.3/12.1, a separate intercept is reported for each grp_no. However, under 9.4/12.3, I do not get separate estimates, only a pooled estimate, almost as if the group= option was ignored.

Steve Denham

Re: Change in behavior of the group= option in the RANDOM statement in PROC GLIMMIX

lvm — Fri, 15 Nov 2013 17:10:22 GMT

I think you should report this to technical support.

Re: Change in behavior of the group= option in the RANDOM statement in PROC GLIMMIX

SteveDenham — Fri, 15 Nov 2013 17:26:44 GMT

- That was my first step after the great gnashing of teeth and rending of garments this morning. Figured I should get the news out here as well.

Steve Denham

Re: Change in behavior of the group= option in the RANDOM statement in PROC GLIMMIX

SteveDenham — Fri, 15 Nov 2013 19:37:07 GMT

I would like to thank the good folks at technical support (and in particular, Jill) for taking the time to work on this, when the problem is with my code and NOT an inconsistency between versions. So, I apologize to everyone, and remind us all: If you do something more than once in a program using a PROC, make sure the syntax used is IDENTICAL, as even small differences can make you insane. And in particular, if you apply the group= option in one RANDOM statement, make sure you apply it in all of the RANDOM statements if you are fitting an AR + random effect model.

Steve Denham