Я работаю с 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, ниже следует:
[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, но это так.Тем не менее, я не мог найти ответ на свой вопрос в Интернете, хотя мне он кажется очень простым.Однако нигде не могу найти значения.Это очень расстраивает: - /