BookmarkSubscribeRSS Feed
☑ This topic is solved. Need further help from the community? Please sign in and ask a new question.
RobertWF1
Quartz | Level 8

I'm running a repeated measures multiple regression for a healthcare dataset. Each member has two records, one pre-intervention and one post-intervention outcome value (# of skilled nursing facility visits).

 

My code looks something like this:

ods output GEEEmpPEst=parms;
proc genmod data=health_data;
   class member_id;
   y = time*class time class x1 x2 x3 / dist=negbin ;
   repeated subject=member_id / type=ind;
run;
ods output close;

 

When I run the procedure the regression coefficients are calculated but the stnd errors, confidence intervals, and p-values are blank.

 

Any ideas what's going wrong?

 

                            Algorithm converged.                                                       


                                                       GEE Model Information

                                    Correlation Structure                           Independent
                                    Subject Effect                  hce_member_id (1199 levels)
                                    Number of Clusters                                     1199
                                    Correlation Matrix Dimension                              8
                                    Maximum Cluster Size                                      8
                                    Minimum Cluster Size                                      2


                            ERROR: Error in estimation routine.                                        


                                                         GEE Fit Criteria

                                                       QIC            0.0000
                                                       QICu           0.0000


                                                Analysis Of GEE Parameter Estimates
                                                 Empirical Standard Error Estimates
 
                                                           Standard   95% Confidence
                              Parameter           Estimate    Error       Limits            Z Pr > |Z|

                              Intercept           223254.9    .        .        .         .      .    
                              time*class          -48.2170    .        .        .         .      .    
                              time                -130.091    .        .        .         .      .    
                              class               150.1242    .        .        .         .      .    
                              x1                  -120.939    .        .        .         .      .    
                              x2                  -11.0064    .        .        .         .      .    
                              x3                   -5.7767    .        .        .         .      .     

 

1 ACCEPTED SOLUTION

Accepted Solutions
ballardw
Super User

Your first clue is likely the text "ERROR: Error in the estimation routine.

 

You haven't provided any example data but I suspect this is going to come down to something in, or not in, the data.

 

Have you tried reducing the data set to maybe 10 or 15 Member_ids and see if that provides more of what you expect? If not, you may want to share the complete log from the trial. That way someone that knows more about Genmod than I do may see something that provides clues.

View solution in original post

1 REPLY 1
ballardw
Super User

Your first clue is likely the text "ERROR: Error in the estimation routine.

 

You haven't provided any example data but I suspect this is going to come down to something in, or not in, the data.

 

Have you tried reducing the data set to maybe 10 or 15 Member_ids and see if that provides more of what you expect? If not, you may want to share the complete log from the trial. That way someone that knows more about Genmod than I do may see something that provides clues.

SAS Innovate 2025: Call for Content

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 25. Read more here about why you should contribute and what is in it for you!

Submit your idea!

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
  • 1 reply
  • 341 views
  • 1 like
  • 2 in conversation