BookmarkSubscribeRSS Feed
Tiny_Kane
Obsidian | Level 7

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

3 REPLIES 3
Tiny_Kane
Obsidian | Level 7

Proc glimmix does not have OUTTEST= option.

And the ConvergenceStatus=CS did not output any information, only one row data, reason=Did not converge., status=1, pdG=1. That's it. 

 

SAS Innovate 2025: Save the Date

 SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!

Save the date!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 3 replies
  • 928 views
  • 0 likes
  • 2 in conversation