I'm running a 3- level multilevel regression model using proc glimmix and one of the classes is HPSU so when i run the following code the estimate for HPSU random effect is 9.13E-11 and the standard error is zero.
Is there something i'm missing?
Here is the code I'm using:
PROC GLIMMIX DATA=DIRT METHOD=LAPLACE ;
CLASS HPSU HNUMBER;
MODEL HEPBRESF (EVENT=LAST) = / CL DIST=BINARY LINK=LOGIT ALPHA=0.1 SOLUTION ODDSRATIO (DIFF=FIRST
LABEL) DDFM=BW ;
WEIGHT HPWGT ;
RANDOM INTERCEPT / SUBJECT=HPSU TYPE=VC SOLUTION CL ;
RANDOM INTERCEPT / SUBJECT=HNUMBER (HPSU) TYPE=VC SOLUTION CL ;
COVTEST / WALD;
RUN;