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

The 2025 SAS Hackathon has begun!

It's finally time to hack! Remember to visit the SAS Hacker's Hub regularly for news and updates.

Latest Updates

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