SASのPHREG手順で区分線形回帰モデルを実装する方法は?
たとえば、X=T でノットが 1 つの場合:
Y = β_10 + β_11 . X ≤ T の場合は X
Y = β_20 + β_21 . X >T の場合 X
連続性の制約があるモデルを考えると、次のようになります。
Y = β_10 + β_11 . X if X ≤ T
Y = β_10 + (β_11 - β_21) T + β_21 . X if X >T
すなわち：
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 に未知のベータ係数があることです。
ご協力いただきありがとうございます！