基于SAS PHREG的分段线性回归

如何在SAS的PHREG程序中实现分段线性回归模型

例如，在X=T处有一个结：

Y=β10+β11。如果X≤ T

Y=β20+β21。如果X>T，则为X

给定具有连续性约束的模型：

Y = β_10 + β_11 . X   if X ≤ T

Y = β_10 + (β_11 - β_21) T + β_21 . X   if X >T

i、 e:

在哪里

最后，我想将其纳入考克斯模型：

Proc PHREG
     Model  time * cas (censure) = X  S_1  ;
Run ;            

但问题是S_1的β系数未知

---

谢谢你的帮助

X是一个时间还是一个分类变量？我忘了提它了。谢谢你的注意。X应该是一个连续变量。

 S_1 =  ( β_11 - β_21 ) T  if X >T  and 0 otherwise.

Proc PHREG
     Model  time * cas (censure) = X  S_1  ;
Run ;