Turn on suggestions

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

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- how to specify strata specific effects in PHREG?

Options

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

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 05-14-2018 02:07 PM
(1452 views)

How to i modify the following sample code to get the strata specific effects for a and b? thanks

proc phreg data = abc;

class a b;

model time*event(0) = a b;

strata c d;

run;

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hi Alexchien,

It is important to understand that the strata statement in PHREG means that the baseline hazard function can depend on the variables you put in the strata statement. The variables you put in the model statement will have same effect across the strata variables - unless ofcourse you have also putted the strata variables in the model line.

If you want the effect of A and B to depend of the level of C and D, then you can do as follows:

```
data simulation;
do i=1 to 1000;
A=rand('bernoulli',0.5);
B=rand('bernoulli',0.5);
C=rand('bernoulli',0.5);
D=rand('bernoulli',0.5);
t=rand('exponential',10);
output;
end;
run;
proc phreg data=simulation;
class A B C D/param=glm;
model t=A*C*D B*C*D;
hazardratio A/at(C=all D=all);
hazardratio B/at(C=all D=all);
run;
*or, if you want each strata to has its own baseline hazard function, you should also use the strata statement;
proc phreg data=simulation;
class A B C D/param=glm;
model t=A*C*D B*C*D;
strata C D;
hazardratio A/at(C=all D=all);
hazardratio B/at(C=all D=all);
run;
```

Notice that you need the "glm" parametrization to make the interactions effects that allows you to use C and D as effectmodifiers.

2 REPLIES 2

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hi Alexchien,

It is important to understand that the strata statement in PHREG means that the baseline hazard function can depend on the variables you put in the strata statement. The variables you put in the model statement will have same effect across the strata variables - unless ofcourse you have also putted the strata variables in the model line.

If you want the effect of A and B to depend of the level of C and D, then you can do as follows:

```
data simulation;
do i=1 to 1000;
A=rand('bernoulli',0.5);
B=rand('bernoulli',0.5);
C=rand('bernoulli',0.5);
D=rand('bernoulli',0.5);
t=rand('exponential',10);
output;
end;
run;
proc phreg data=simulation;
class A B C D/param=glm;
model t=A*C*D B*C*D;
hazardratio A/at(C=all D=all);
hazardratio B/at(C=all D=all);
run;
*or, if you want each strata to has its own baseline hazard function, you should also use the strata statement;
proc phreg data=simulation;
class A B C D/param=glm;
model t=A*C*D B*C*D;
strata C D;
hazardratio A/at(C=all D=all);
hazardratio B/at(C=all D=all);
run;
```

Notice that you need the "glm" parametrization to make the interactions effects that allows you to use C and D as effectmodifiers.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

thanks for the explanation.

📢

**ANNOUNCEMENT**

The early bird rate has been extended! Register by March 18 for just $695 - $100 off the standard rate.

Check out the agenda and get ready for a jam-packed event featuring workshops, super demos, breakout sessions, roundtables, inspiring keynotes and incredible networking events.** **

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.