turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- SAS Procedures
- /
- Contrast statement in Phreg for linear trend

Topic Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

12-16-2013 02:41 PM

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.

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to 57078

12-18-2013 02:50 PM

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