1

私は 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 である必要があると思いますが、それについては間違っている可能性があり、他の切片と勾配を計算する方法はまだわかりません。どんな助けでも大歓迎です。

4

0 に答える 0