我有一个三件式线性回归模型,我在R中运行,以模拟大量人口中的体重 . 我的数据集称为hdata . 通过贯穿断点的所有组合的迭代过程,我发现了与最低残差平方误差模型相关的两个断点 . 我的分段回归的代码,指定了两个断点:

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

根据这些估计,我想计算三个截距,以及与模型相关的三个斜率,但我不知道如何做到这一点 . 为简单起见,我调用与(拦截)Est1相关的估计,与hdata $ age Est2相关的估计,依此类推......直到Est7 . 我认为第一个拦截应该是Est1 Est3,第一个斜率应该是Est2 Est6,但我可能错了,仍然不知道如何计算其他拦截和斜率 . 任何帮助,将不胜感激 .