BookmarkSubscribeRSS Feed
EJ0603
Calcite | Level 5

I'm creating a model using proc logistic with a code of the following form:

 

proc logistic data=work.MODEL_DATASET;

class cat_var1 cat_var2 cat_var3 cat_var4 cat_var4 cat_var5 cat_var6 cat_var7 /order=freq desc param = glm;

model outcome(event='1')= con_var1 convar2 cat_var1 cat_var2 cat_var3 cat_var4 cat_var4 cat_var5 cat_var6 cat_var4*cat_var6;

 

I changed the code slightly so that it used reference parameterization i.e. the only change I made was to replace param = glm with param = ref.  I was expecting to see almost identical results, but for some of the variables I'm seeing different parameter estimates and odds ratios.  Can anyone explain why this is?  I'm using SAS 9.4.

4 REPLIES 4
Reeza
Super User

I'm assuming the variables that are different are the class variables?

 

if so, that makes sense as you've changed your parameterization method and the odds ratio is in comparison to a level. If you've changed that level your estimate would change.  Your p-values should be the same. 

EJ0603
Calcite | Level 5

Thanks for the quick response - the variables that are different are some (but not all) of the class variables.  However, as far as I can see the same reference level has been used in both parameterizations and so I don't understand why the coefficients have changed.  Also, the value of the intercept has changed. 

EJ0603
Calcite | Level 5

Also, the odds ratios differ between the two parametrizations.  With the param=GLM, there is a convergence warning message which states "WARNING: The information matrix is singular and thus the convergence is questionable.  Try specifying a larger SINGULAR= value."  There isn't a convergence warning message when param=REF is used. 

 

Could this be the reason for the different odds ratios and parameter estimates between the two options?  When would there be a convergence issue with one parameterization and not another?

 

Reeza
Super User

I wouldn't expect convergence issues. 

Ideally you could provide sample data so we could replicate the issue. 

 

If this isn't possible you should contact SAS Tech Support for help. 

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
  • 4 replies
  • 2277 views
  • 2 likes
  • 2 in conversation