如何获得分段的R包中所有回归线的截距的95%CI限制?

时间:2018-10-08 10:23:19

标签: r

目前,我正在使用“分段”的R包。

我有简单的数据,其中y通过与一个断点的线性关系依赖于x。 X的范围是20到100,y的范围是8到16。断点在62。

我可以得到(或从标准误差中估计)斜率和第一条回归线(在坐标系统中更接近于0,在断点左侧)的截距的置信区间。

但是我只能得到第二个回归线的斜率及其置信区间和截距本身在R之外(那个,即在我的断点右边,因此远离0),而不是拦截的CI。

如何获取两条回归线的截距+置信区间的上限和下限?

    out.lm<-lm(b~a)
    library(segmented)
    segreg<-segmented(out.lm, seg.Z=~a)
    slope(segreg)
    confint(segreg)
    intercept(segreg)
    print.segmented(segreg)
    predict(segreg)
    summary(segreg)

这是b的样子,下面是a:

 [1] 14.100000 15.300000 16.833333 10.333333 12.533333 12.866667 14.433333 16.000000
  [9] 14.533333 12.233333 15.533333 11.633333 15.733333 12.766667 14.200000 15.200000
 [17] 13.300000 15.500000 12.700000 13.800000 13.866667 13.866667 13.600000 15.566667
 [25] 15.533333 13.900000 14.200000 14.733333 13.800000 14.233333 13.833333  7.733333
 [33] 12.633333 12.833333 13.700000 13.266667 11.900000 14.400000 10.900000 11.300000
 [41] 12.833333 14.200000 12.266667 13.100000 14.633333 12.200000 12.266667  8.166667
 [49] 11.466667  9.833333 13.133333 13.766667 12.466667 12.633333 12.333333 12.333333
 [57] 13.066667 11.033333 11.500000 12.133333 14.166667 12.333333 10.966667 11.433333
 [65] 13.333333 12.033333  9.866667 11.966667  9.800000 10.300000 10.766667  7.833333
 [73] 13.766667 11.433333  8.466667  9.933333  6.466667 10.233333  7.300000 10.533333
 [81] 10.433333  9.400000 10.300000  9.133333  9.633333  9.600000  8.466667  7.733333
 [89]  9.233333  8.600000 11.400000  7.233333  8.800000  7.700000  8.533333 10.000000
 [97]  6.933333 10.200000  4.700000  9.100000 10.033333  8.066667  9.266667 10.666667
[105]  9.700000  9.333333  7.900000  9.166667
> print(a)
  [1] 33 33 32 31 31 30 30 30 29 29 29 29 29 29 29 28 28 28 27 27 27 27 27 27 27 27 27 27 27
 [30] 27 27 27 26 26 26 26 26 26 26 26 26 26 26 25 25 25 25 25 25 24 24 24 24 24 24 24 24 24
 [59] 24 24 23 23 23 23 23 23 22 22 22 22 22 22 21 21 21 21 21 21 21 21 21 21 21 21 20 20 20
 [88] 20 19 19 19 19 19 19 19 18 18 18 18 18 18 18 17 17 16 16 16 16

对不起,如果我在这里做错了什么。我从未在任何地方发布过有关编程的任何文章,这实际上不关我的事。我了解基本的R,但仅此而已。但是,尽管对我来说这似乎很简单,但是我在互联网上的任何地方都找不到答案。但是,我在任何地方都找不到值。这非常令人沮丧:-/

0 个答案:

没有答案