私は R で実行しているスリーピースの線形回帰モデルを持っており、大規模な母集団で年齢を超えた体重をモデル化しています。私のデータセットは hdata と呼ばれます。ブレーク ポイントのすべての組み合わせを実行する反復手順を通じて、最小の残差二乗誤差モデルに関連付けられた 2 つのブレークポイントを見つけました。2 つのブレークポイントが指定された区分回帰のコードは次のとおりです。
piecewise=lm(hdata$weight ~ hdata$age*(hdata$age < 0.7) + hdata$age*(hdata$age>=0.7 & hdata$age<2) + hdata$age*(hdata$age>2))
私が見たとき:
summary(piecewise)
次の出力が得られます。
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1723.012 210.042 8.203 3.67e-16 ***
hdata$age 18.494 1.089 16.975 < 2e-16 ***
hdata$age < 0.7TRUE -1690.051 210.476 -8.030 1.48e-15 ***
hdata$age >= 0.7 & hdata$age < 2TRUE -721.200 213.882 -3.372 0.000758 ***
hdata$age > 2TRUE 478.094 210.194 2.275 0.023016 *
hdata$age:hdata$age < 0.7TRUE 2022.896 45.477 44.481 < 2e-16 ***
hdata$age:hdata$age >= 0.7 & hdata$age < 2TRUE 603.453 30.532 19.764 < 2e-16 ***
hdata$age:hdata$age > 2TRUE NA NA NA NA
これらの推定値から、モデルに関連付けられた 3 つの切片と 3 つの勾配を計算したいのですが、その方法がわかりません。簡単にするために、(Intercept) Est1 に関連付けられた推定値、hdata$age に関連付けられた推定値を Est2、というように Est7 まで呼び出します。最初の切片は Est1 + Est3 で、最初の勾配は Est2 + Est6 である必要があると思いますが、それについては間違っている可能性があり、他の切片と勾配を計算する方法はまだわかりません。どんな助けでも大歓迎です。