目前,我正在使用“分段”的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,但仅此而已。但是,尽管对我来说这似乎很简单,但是我在互联网上的任何地方都找不到答案。但是,我在任何地方都找不到值。这非常令人沮丧:-/