Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- WARNING: Obtaining minimum variance quadratic unbiased estimates as st...

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 05-14-2017 06:27 AM
(8994 views)

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;

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

- Tags:
- glimmix

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

9 REPLIES 9

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Why you include subject term in the independent variables ?

Try option type=chol.

random intercept/ subject=fid type=chol ;

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

I am glad that worked and pround of your effort .

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

tissue5454

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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.

Are you ready for the spotlight? We're accepting content ideas for **SAS Innovate 2025** to be held May 6-9 in Orlando, FL. The call is **open **until September 16. Read more here about **why** you should contribute and **what is in it** for you!

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.