For binary response data, PROC GENMOD fits the following model by default:
where p is the probability of the response level identified as the first level in the "Response Profile" table in the output and is the normal cumulative distribution function. By default, the covariate vector x contains an intercept term. This is sometimes called Abbot’s formula.
Because of the symmetry of the normal (and logistic) distribution, the effect of reversing the order of the two response values is to change the signs of in the preceding equation.
By default, response levels appear in ascending, sorted order (that is, the lowest level appears first, and then the next lowest, and so on). There are a number of ways that you can control the sort order of the response categories and, therefore, which level is assigned the first ordered level. One of the most common sets of response levels is {0,1}, with 1 representing the event with the probability that is to be modeled.
Consider the example where Y
takes the values 1 and 0 for event and nonevent, respectively, and EXPOSURE is the explanatory variable. By default, PROC
GENMOD assigns the first ordered level to response level 0, causing the probability of the nonevent to be modeled. There are
several ways to change this.
Besides recoding the variable Y
, you can do the following:
Explicitly state which response level is to be modeled by using the response variable option EVENT= in the MODEL statement:
model Y(event='1') = Exposure;
Specify the nonevent category for the response variable in the response variable option REF= in the MODEL statement:
model Y(ref='0') = Exposure;
Specify the response variable option DESCENDING
in the MODEL
statement to assign the lowest ordered value to Y
=1:
model Y(descending)=Exposure;
Assign a format to Y
such that the first formatted value (when the formatted values are put in sorted order) corresponds to the event. For the
following example, Y
=0 could be assigned formatted value ‘nonevent’ and Y
=1 could be assigned formatted value ‘event.’ Since ORDER=FORMATTED by default, Y
=1 becomes the first ordered level.
proc format; value disease 1='event' 0='nonevent'; run; proc genmod; model y=exposure / dist=bin; format y disease.; run;
Arrange the input data set so that Y
=1 appears first and use the ORDER=DATA option in the PROC GENMOD statement. Because ORDER=DATA sorts levels in order of their
appearance in the data set, Y
=1 becomes the first ordered level. Note that this option causes classification variables to be sorted by their order of appearance
in the data set, also.