Tests whether hazard function of different groups are equal, controlling for confounding covariates by partitioning into homogeneous strata.
Definition
Let:
- : number of strata defined by a set of covariates
- : number of groups being compared within each stratum
- : log-rank statistic for group in stratum
- : variance-covariance in stratum
Test Hypotheses
Procedure
For each stratum , compute and as in the k-sample log-rank test.
Global Test
Pool across strata:
For (two groups):
Interpretation
The stratified test controls for confounding covariates by testing within homogeneous strata. For example, comparing treatments while controlling for disease stage: compute the log-rank statistic separately within each stage, then pool the results. This produces a treatment effect estimate adjusted for stage.
Example: Controlling for Sex in Larynx Cancer
Goal: Compare survival across larynx cancer stages while controlling for sex.
# Stratified log-rank: stage effect adjusted for sex
fit_strat <- survdiff(Surv(time, status) ~ as.factor(stage) + strata(sex), data = dat)
fit_stratProcedure:
- Partition data into strata (male, female)
- Within each stratum, compute and as in the k-sample test
- Pool across strata: ,
- Global test:
Interpretation: A significant result means that stage differences in survival persist after controlling for sex. The stratified test answers: “Is the stage effect still there, even when accounting for sex differences?”
Interpretation
The stratified test controls for confounding by testing within homogeneous strata, then pooling results globally. It’s a bridge between non-parametric methods (stratification) and regression (Cox PH).