创建含时变协变量的KM曲线与Cox PH模型，及数据编码正确性咨询

我将部分接受治疗的患者设为试验组，未接受治疗的设为对照组。为规避永生时间偏倚（immortal time bias），我把治疗作为时变协变量进行编码，想咨询这种编码方式是否正确。

举两个患者的情况说明：

• 患者P1：研究启动30天后接受治疗，接受治疗30天后死亡；
• 患者P2：未接受治疗，研究结束时（第80天）仍存活。

我对上述数据的编码如下：

ID  Treatment   Time    Death
P1  0           30      0
P1  1           30      1
P2  0           80      0

随后我使用以下SAS代码绘制KM生存曲线：

/*Survival curve*/
proc lifetest data=treatment;
time Time*Death(0); /*注：原代码中"Dearg"应为笔误，修正为"Death"*/
strata Treatment;
run;

同时运行Cox比例风险模型：

/*Cox PH Model*/
proc phreg data=treatment;
model Time*Death(0)=treatment;
run;

内容的提问来源于stack exchange，提问作者user1813615