BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
Tissue5454
Obsidian | Level 7

Hi all, 

I am trying to do logistic regression with random effiect for binary data by GLIMMIX.

However, I get warinig message as follows;

 

WARNING: Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed

 

Anyone know how I can overcome this problem?

Thanks,

 

--------------------------------------------------------------------------------

proc glimmix data=logitbank method=quad;
class fid;
model ydat2(event='1') = D1013 fid*D1013
/dist = binary link=logit ddfm=none solution;
random intercept/ subject=fid;
output out =glmout pred =xbeta pred(ilink) =predprob;
covtest;
run;

---------------------------------------------------------------------------------

 

1 ACCEPTED SOLUTION

Accepted Solutions
Ksharp
Super User
Or try PARAM statement . and I quoted here.
"The PARMS statement specifies initial values for the covariance or scale parameters, or it requests a grid search over several values of these parameters in generalized linear mixed models."

View solution in original post

9 REPLIES 9
Ksharp
Super User

Why you include subject term in the independent variables ?

Try option type=chol.

 

random intercept/ subject=fid type=chol ;

Tissue5454
Obsidian | Level 7

Hi Ksharp!

 

Thank you for reply!

 

>Why you include subject term in the independent variables ?

 

Because I want to include not only random intercept but also random slope.

 

>Try option type=chol.

 

I received same warning message: WARNING: Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed...

Tissue5454

 

 

Ksharp
Super User
Or try PARAM statement . and I quoted here.
"The PARMS statement specifies initial values for the covariance or scale parameters, or it requests a grid search over several values of these parameters in generalized linear mixed models."

Tissue5454
Obsidian | Level 7

Hi Ksharp!

 

Thank you for reply!

 

I tried PARAM statement and succeed in getting result!

I think the warning message;

"WARNING: Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed" tend to come up when the model is too complicated.

In my case, independent variable "D1013" is nonsignificant .

 

Thank you again,

 

tissue5454

 

rplum64
Calcite | Level 5

Hey there, tissue5454

 

Do you mind sharing the syntax of your PARAM fix. I am having a similar code and error log, and would like to know how/where the PARAM fix was inserted. 

 

Many thanks, in advance.

 

-R

Tissue5454
Obsidian | Level 7

Sorry... I could not find the above file.

 

tissue5454

 

sld
Rhodochrosite | Level 12 sld
Rhodochrosite | Level 12

Lots of good information in this recent SGF paper by Kathleen Kiernan (2018) of SAS, including possible solutions to the warning message that you are getting

 

Insights into Using the GLIMMIX Procedure to Model Categorical Outcomes with Random Effects

 

This older paper by Kiernan, Tao, and Gibbs (2012) is good, too

 

Tips and Strategies for Mixed Modeling with SAS/STAT® Procedures

Tissue5454
Obsidian | Level 7

Thank you sld !


Your information is valuable for me and rplum64.

These papers containe some examples which use "PARMS" statement.

 

Thank you again for your help.

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
  • 9 replies
  • 8197 views
  • 6 likes
  • 4 in conversation