BookmarkSubscribeRSS Feed
RachelB_
Calcite | Level 5

Hello

I'm trying proc glimmix where the outcome is dichotomous and the model does not converge.  Are there any tips for setting up the proc that will make it more likely to converge?  The model does converge when the outcome is continuous.

Thanks,

Rachel

2 REPLIES 2
SteveDenham
Jade | Level 19

First, take a look at the paper by Kathleen Kiernan, Jill Tao and Phil Gibbs from SGF 2012:

http://support.sas.com/resources/papers/proceedings12/332-2012.pdf

Some of what I say will follow their recommendations.  First, when you say it doesn't converge, do you mean that the convergence history is "smooth" in the  objective function, or does it jump around?  If it is smooth, then it may be just a matter of increasing the number of iterations (NLOPTIONS statement), or controlling the outer loop convergence (pconv= option in the PROC GLIMMIX statement).  However, if it is jumpy, then it is time to look at section II in the Kiernan paper.  The first thing I think of is changing the technique to NRRIDG or NEWRAP.  These work better than the default QUANEW for Bernoulli data (I assume that the response for each individual is 0/1).

I encourage to share pertinent parts of the output and program with us, so that we might give a more complete answer.

Good luck.

Steve Denham

RachelB_
Calcite | Level 5

Thank you Steve for steering me to that paper.  The model converged when I increased the iterations, and changed the technique to NEWRAP.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

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
  • 2 replies
  • 1315 views
  • 0 likes
  • 2 in conversation