BookmarkSubscribeRSS Feed
57078
Calcite | Level 5

Hi,

I have a 6 level categorical variable and want a linear trend for the hazard ratios.

Should the contrast statement have 5 levels ignoring the reference group?

Code:

proc phreg data=trcs ;

class rumol5A  /order=internal ;

model age*censor1(1) =rumol5A /rl;

contrast ' LINEAR EFFECT'   rumol5A -2 -1 0 1 2  /e;

run;

I've found online that these are the numbers to use for a linear effect for 5 coefficients but I can't find an easy explanation of why. Can someone confirm the above and provide an explanation?

Thanks in advance.

1 REPLY 1
SteveDenham
Jade | Level 19

I think under the default parameterization (REF, with the reference level being the last level), your code looks right.  It excludes the last level.  However, that may not be what you want or need.  What values does rumol5A take on?  You may need to change the parameterization to get what you want.  I know if you used:

class rumol5A/order=internal param=GLM;

then the correct contrast would be:

contrast 'Linear effect' rumol5A -5 -3 -1 1 3 5 /e;

as the GLM parameterization requires all levels to be included.

As far as the coefficients, they provide a linear only contrast that has coefficients that sum to zero.  To derive them, take the ordinal values of the levels, find their mean, subtract the mean from the ordinal value, and convert all to integers by multiplying each by an appropriate constant.

Steve Denham

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 Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 1 reply
  • 2546 views
  • 0 likes
  • 2 in conversation