Hi there,
I was trying to fit an ordinal mixed effects model using PROC GLIMMIX procedure. The code is as below. The outcome Y has 3 ordered value (1,2,3). Y_baseline is the outcome measured at baseline time. The model includes random intercept and random-effect of time, and the baseline outcome as covariate. The data includes 53 subjects with 181 observations. However the model did not converge, even if I combined the three values into 2 values and run regular logistic mixed effects mode, it still did not converge. Is it due to that the sample size is too small? Is there any solution? Thank you in advance.
PROC GLIMMIX data=final1;
CLASS pt Y_baseline(ref=first);
MODEL Y(desc order=internal)= Y_baseline days/DIST=multinomial LINK=cumlogit SOLUTION CL;
RANDOM int days/subject=pt;
RUN;
Kane
