I was able to run your model without errors, but there are some results that make me wonder if your model could be over specified. Here is the code I used:
PROC GLIMMIX DATA = test1/*(where=(stratclus<999))*/ METHOD=QUAD(FASTQUAD QPOINTS=20) plots=all;
CLASS RAHHIDPN STRATCLUS GENDER(REF='1') RACE_ETH(ref='1') MARITAL(ref='1');
PARMS (0.01) (1.72) (0.5);
MODEL MOBILITY = BSAGE_C BSAGE_C*BSAGE_C GENDER GENDER*BSAGE_C RACE_ETH RACE_ETH*BSAGE_C RACE_ETH*GENDER
TIME TIME*TIME TIME*BSAGE_C TIME*GENDER TIME*RACE_ETH/solution obsweight=HCNSWGTR ;/**LEVEL 1 VARIABLES**/
RANDOM INTERCEPT /subject=STRATCLUS type=VC ;
RANDOM INTERCEPT /subject=RAHHIDPN(STRATCLUS) type=VC ;
NLOPTIONS GCONV=0.00001 tech=congra;
COVTEST/WALD;
ODS OUTPUT parameterestimates=FIXED1 covparms=RAND1;
RUN;
The two changes were to reduce the number of quadrature points, and to use the PARMS statement to provide some better starting values for the random effects. Regarding my thought of over-specification, the F values for five of the terms in the model went to infinity. These were GENDER, RACE_ETH*GENDER, TIME*TIME, TIME*RSAGE_C and TIME*GENDER. You may want to look to see if there is something unusual with the strata and the number of observations available in the various cells.
SteveDenham
... View more