* Stata
stset time, failure(event)               // event=1 事件，0 删失
sts graph, by(group)                     // Kaplan-Meier 生存曲线
sts test group                           // log-rank 检验
stcox x1 x2                              // Cox 比例风险模型（HR）
estat phtest, detail                     // PH 假设检验

## R
library(survival)
fit <- survfit(Surv(time, event) ~ group, data = df)
plot(fit); survdiff(Surv(time, event) ~ group, data = df)
cfit <- coxph(Surv(time, event) ~ x1 + x2, data = df)
summary(cfit); cox.zph(cfit)             // PH 假设检验