R-如何用程序在分段函数中获得断点



使用mtcars数据集,我试图确定mpg的折线回归拟合,作为hp和wt的函数,并带有断点。这是代码:

mpg = mtcars$mpg
wt = mtcars$wt
hp = mtcars$hp
reg = lm (mpg ~ hp +wt)
hp_thresh = 150
wt_thresh = 3
library(segmented)
seg.o = segmented (reg,seg.Z=~hp+wt, psi=list (hp=hp_thresh, wt=wt_thresh))
summary (seg.o)

如何以编程方式获取断点?

弄清楚了这一点,下面的技巧:

summary.segmented(seg.o)$psi [1,2]
summary.segmented(seg.o)$psi [2,2]

最新更新