1
$\begingroup$

Im building a Cox Regression in SAS. My code looks like this:

proc phreg data = mydata; 
   class sex hospital;
   model time*status(0) = sex hospital age bmi;
   strata treatment_variable;
run;

Now this code does not calculate a hazard ratio for treatment A vs treatment B. If I remove the strata statement and instead include the treatment_variable inside the class and model statement I get a hazard ratio for the variable. The coefficients for the other variables in the model changes slightly when I do this.

My questions are:

  1. What purpose does the strata statement serve? Why not include the treatment variable in the model? Most of the time you are interested in the effect of treatment A vs treatment B when controlling for some other variables.
  2. Why does the coefficients of the other variables change?

Hope my question is clear. Thanks!

$\endgroup$

1 Answer 1

3
$\begingroup$

Stratification fits separate baseline hazards to each stratum and combines their likelihoods. This can address violation of proportional hazards across the stratification variable, but the drawbacks are that you may lose a bit of estimation efficiency and, as you discovered, you cannot actually draw inference on the stratification variable.

If proportional hazards are violated for your parameter of interest (treatment) then it would make sense to look into time-dependent covariates.

$\endgroup$
3
  • $\begingroup$ So I should stratify on treatment if the proportional hazards assumption for the treatment variable is violated? If that is not the case and if I would like to measure the effect of treatment A vs treatment B while controlling for my covariates, I guess I should not stratify but rather include the treatment variable in the model? Thus I can say that "patients with treatment A is 50% more likely to experience the event compared to patients with treatment B while controlling for sex, hospital, age and bmi"? $\endgroup$
    – Jam.Wil
    Commented Oct 26, 2023 at 7:22
  • $\begingroup$ Also, if I stratify on treatment. What does the HR of my covariates tell me? $\endgroup$
    – Jam.Wil
    Commented Oct 26, 2023 at 7:23
  • $\begingroup$ The same as they would in the non-stratified model, only now no longer assuming the same baseline hazard across treatment groups. $\endgroup$
    – PBulls
    Commented Oct 26, 2023 at 14:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.