I am trying to model the data using a two-level model. Student nested in school. My dependent variable "pass" is dichotomously coded (1, 0). My independent variable includes both dichotomously coded variable first and status (1, 0) and continuous variable "practice score" as below
PROC GLIMMIX
DATA=pred METHOD=LAPLACE NOCLPRINT;
CLASS SCHOOLNUM;
MODEL PASS (EVENT="1")=first practice_score status/CL DIST=BINARY LINK=LOGIT SOLUTION;
RANDOM INTERCEPT/ SUBJECT=SCHOOLNUM S CL TYPE=VC;
COVTEST /WALD;
run;
Based on the fixed effects estimates as in the picture below, the probability of pass for first=1, status=1, is exp(-5.7755+0.7758+0.00872461.1849)/(1+exp(-5.7755+0.7758+0.008724+1.1849))=0.02175, which is too low from expected.
Where was wrong with my SAS code above? Any help would be much appreciated. Thanks!
... View more