Расчет прогнозируемого интервала

9

У меня есть следующие данные, расположенные здесь . Я пытаюсь рассчитать 95% доверительный интервал для средней чистоты, когда процент углеводородов равен 1,0. В R я ввожу следующее.

> predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95)
   fit      lwr      upr
1 89.66431 87.51017 91.81845

Тем не менее, как я могу получить этот результат сам? Я попытался использовать следующее уравнение.

s_{N е вес} знак равно \sqrt{s^{2} (1 + \frac{1}{N} + \frac{({Икс}_{N е вес} - \bar{Икс})^{2}}{Σ ({Икс}_{я} - \bar{Икс})^{2}})}

$s_{new}=\sqrt{s^2\left(1+\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$

И я ввожу следующее в R.

> SSE_line = sum((purity - (77.863 + 11.801*hydro))^2)
> MSE = SSE_line/18
> t.quantiles <- qt(c(.025, .975), 18)
> prediction = B0 + B1*1
> SE_predict = sqrt(MSE)*sqrt(1+1/20+(mean(hydro)-1)^2/sum((hydro - mean(hydro))^2))
> prediction + SE_predict*t.quantiles
[1] 81.80716 97.52146

Мои результаты отличаются от функции предсказания R. Что я неправильно понимаю в интервалах прогнозирования?

r regression confidence-interval prediction-interval

— idealistikz
источник

Как вы рассчитываете MSE в своем коде?

Я добавил расчет в пост.

— idealistikz

1

как MMJ предложил вам попробоватьpredict(purity.lm, newdata=list(hydro=1.0), interval="prediction", level=.95)

— Vinux

16

Ваш predict.lmкод рассчитывает доверительные интервалы для установленных значений. Ваша рука вычисляет интервалы прогнозирования для новых данных. Если вы хотите получить тот же результат, predict.lmчто и при расчете руки, измените interval="confidence"на interval="prediction"

— MMJ
источник

1

Хороший ответ от dpel. Я бы добавил, что разницу между доверительным интервалом и интервалом прогнозирования можно указать, как показано ниже:

s_{N е вес} знак равно \sqrt{s^{2} (\frac{1}{N} + \frac{({Икс}_{N е вес} - \bar{Икс})^{2}}{Σ ({Икс}_{я} - \bar{Икс})^{2}})}

$s_{new}=\sqrt{s^2\left(\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$

s_{N е вес} знак равно \sqrt{s^{2} (1 + \frac{1}{N} + \frac{({Икс}_{N е вес} - \bar{Икс})^{2}}{Σ ({Икс}_{я} - \bar{Икс})^{2}})}

$s_{new}=\sqrt{s^2\left(1+\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$

Источник См слайд страницы 5/17 и 11/17

— lklklk
источник