0

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

例如,在 X=T 处有一个结:

Y = β_10 + β_11 。X 如果 X ≤ T

Y = β_20 + β_21 。X 如果 X >T

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

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

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

IE :

Y= β_0 + β_1 . X + S_1  

在哪里

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

最后,我想将它包含在 Cox 模型中:

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

但问题是 S_1 中有未知的 beta 系数。

谢谢你的帮助!

4

0 回答 0