BookmarkSubscribeRSS Feed
abened
Calcite | Level 5

I am trying to estimate a pooled random effects proportion. I get an error message I have never encountered before:

NOTE: At the starting values the G matrix is the zero matrix. Adaptively determining the number of quadrature points is not meaningful. The QMIN=1 value is chosen. You can provide starting values for the covariance parameters with

      the PARMS statement.

It also says: (which I understand better)

NOTE: Convergence criterion (ABSGCONV=0.00001) satisfied.

NOTE: Estimated G matrix is not positive definite.

 

The confidence interval obtained is really, really, really wide. 

 

The code I used is:

PROC GLIMMIX DATA=IPDALL2 METHOD=QUAD;

CLASS COUNTRY;  

MODEL NUM/DENOM= /DIST=BIN SOLUTION CL;

RANDOM INT /SUBJECT=COUNTRY;

 

RUN;

 
 
 

 

2 REPLIES 2
ballardw
Super User

That is not an error message but a statement about how the procedure behaved. It is likely indicating that method=quad was not appropriate for the specific data used in this case.

SteveDenham
Jade | Level 19

@ballardw's response is right on.  Method=quad isn't appropriate--I would guess that the number of countries being fit or at least with very different proportions is small.  You could try METHOD=LAPLACE (which is what QUAD with one point is doing anyway), but the whole thing looks as if there is insufficient variability due to country for that to work well (that's the meaning of NOTE: Estimated G matrix is not positive definite.), and it is why the confidence interval is so large

 

If that doesn't help, you'll probably have to go to a marginal model, using the default pseudo-likelihood method.

 

Steve Denham

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
  • 2000 views
  • 0 likes
  • 3 in conversation