BookmarkSubscribeRSS Feed
deleted_user
Not applicable
Hello,

I have two categorical variables (A and B) that I want to add as Random factor in GLIMMIX. A has 4 categories and B has 44. Finally B is nested in A.
When I added as random factor, the model results said that "Estimate G matrix is not positive definite". How can I solve this problem? Did I add the variables in a correct way??

proc GLIMMIX data=Resumen;
class A B Cat;
model total = intra flux Cat length / dist=poisson link=log solution;
random B B(A);
run;
1 REPLY 1
SteveDenham
Jade | Level 19
While this may not solve the problem completely, note that the statement:

random B B(A);

specifies factor B twice. That may be the source of the error message.

Try

random A B(A);

This should specify A as a random effect, and B within A as a second effect.

Good luck.

hackathon24-white-horiz.png

2025 SAS Hackathon: There is still time!

Good news: We've extended SAS Hackathon registration until Sept. 12, so you still have time to be part of our biggest event yet – our five-year anniversary!

Register Now

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

SAS Training: Just a Click Away

 Ready to level-up your skills? Choose your own adventure.

Browse our catalog!

Discussion stats
  • 1 reply
  • 1647 views
  • 0 likes
  • 2 in conversation