Chapter Contents |
Previous |
Next |
The CATMOD Procedure |
PROC CATMOD can compute three different types of logits with the use of keywords in the RESPONSE statement. Other types of response functions can be generated by specifying appropriate transformations in the RESPONSE statement.
proc catmod; model r=a|b; run;
proc catmod; direct x1 x2 x3; model r=x1 x2 x3; run;
Since the preceding statements do not include a RESPONSE statement, generalized logits are computed. See Example 22.3 for another example.
When the dependent variable has two responses, the parameter estimates from the CATMOD procedure are the same as those from a logistic regression program such as PROC LOGISTIC (see Chapter 39, "The LOGISTIC Procedure"). The chi-square statistics and the predicted values are also identical. In the two-response case, PROC CATMOD can be made to model the probability of the maximum value by either (1) organizing the input data so that the maximum value occurs first and specifying ORDER=DATA in the PROC CATMOD statement or (2) specifying cumulative logits (CLOGITS) in the RESPONSE statement.
Caution: Computational difficulties may occur if you use a continuous variable with a large number of unique values in a DIRECT statement. See the "Continuous Variables" section for more details.
proc catmod; response clogits; direct x; model r=a x; run;
The preceding statements correspond to a simple analysis that addresses the question of existence of an association between the independent variables and the ordinal dependent variable. However, there are some commonly used models for the analysis of ordinal data (Agresti 1984) that address the structure of association (in terms of odds ratios), as well as its existence.
If the independent variables are class variables, a typical analysis for such a model uses the following statements:
proc catmod; weight wt; response clogits; model r=_response_ a b; run;
On the other hand, if the independent variables are ordinally scaled, you might specify numeric scores in variables x1 and x2, and use the following statements:
proc catmod; weight wt; direct x1 x2; response clogits; model r=_response_ x1 x2; run;
Refer to Agresti (1984) for additional details of estimation, testing, and interpretation.
Chapter Contents |
Previous |
Next |
Top |
Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.