如何在 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 系数。
谢谢你的帮助!