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-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

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