I would like to evaluate the goodness-of-fit in an ordinal logistic regression model (proportional odds). My outcome has 4 levels (1-4 with 4 being a more severe outcome). My understanding is that in order to run the Hosmer-Lemeshow Goodness-of-Fit Test in SAS, your outcome has to be binary. When I request this test in SAS 9.4 I get the message "LACKFIT is ignored since there are more than 2 response levels". However, when I run the same model in SAS enterprise v7.15, I get the an output for the test. This is what my code looks like (with some additional class variables that I didn't write down: proc logistic data=cohort;
class gender(ref='1') specialist_visit(ref='1') gp_visit(ref='1') ...../param=ref;
model outcome (descending)= gender specialist_visit gp_visit .../lackfit;
run; This is the output from sas enterprise ( I just listed the first and last line): Group Total Observed outcome = 4 Observed outcome = 3 Observed outcome = 2 Observed outcome = 1 Expected outcome = 4 Expected outcome = 3 Expected outcome = 2 Expected outcome = 1 1 100 50 25 12 10 40.9 30.6 12.1 15.9 … … … … … … … … … … 10 100 10 15 20 55 13.2 22.5 18.0 46.7 Hosmer and Lemeshow Goodness-of-Fit Test Chi-Square DF Pr > ChiSq 36.3246 26 0.0859 My question is: how is sas enterprise producing this outcome as I couldn't find any documentation on it. Is there an update to the enterprise version to allow this test to run there but not in sas 9.14? Thank you
... View more