I'm running PROC GENMOD to do a Poisson regression in SAS Software 9.4 (TS1M8) with SAS/STAT 15.3. The dependent variable is a 0-1 variable. The MODEL statement has options DIST=POI LINK=LOG.
One of my CLASS variables has four values, 0, 1, 2, 3, with ORDINAL (THERMOMETER) parametrization.
The Class Level Information for that variable (which I'm calling Ordinal_variable here) in the output of PROC GENMOD is
CLASS
Value
Design Variables
Ordinal_variable
0
0
0
0
1
1.000
0
0
2
1.000
1.000
0
3
1.000
1.000
1.000
Note that there are other variables (including continuous ones) in my MODEL statement, so I can't just use PROC FREQ to get risk ratios.
I'd like to use ESTIMATE statements (with options /E EXP) to get risk ratios for the comparisons between successive levels of the ordinal class variable (that is, risk ratios of the dependent variable for Ordinal_variable equaling 1 vs 0, 2 vs 1, and 3 vs 2).
My question: I can't figure out how what the arguments for the ESTIMATE statement should be.
(NOTE: I already know how to use reference parametrization and ESTIMATE statements to get risk ratios for 1 vs 0, 2 vs 0, and 3 vs 0 -- have the parametrization for that variable in the CLASS statement be ORDER=INTERNAL DESC PARAM=REF and then use three estimate statements, ESTIMATE that_variable 0 0 1 /E EXP; ESTIMATE that_variable 0 1 0 /E EXP; ESTIMATE that_variable 1 0 0 /E EXP;
What I'd like to do here, though, is take advantage of the ordinal nature of my predictor to get risk ratios for the successive steps.)
... View more
