Weźmy następujący przykład:
set.seed(342)
x1 <- runif(100)
x2 <- runif(100)
y <- x1+x2 + 2*x1*x2 + rnorm(100)
fit <- lm(y~x1*x2)
Tworzy to model y oparty na x1 i x2 przy użyciu regresji OLS. Jeśli chcemy przewidzieć y dla danego x_vec, moglibyśmy po prostu użyć wzoru, który otrzymujemy z summary(fit)
.
Co jednak, jeśli chcemy przewidzieć dolną i górną prognozę y? (dla danego poziomu ufności).
Jak więc zbudowalibyśmy formułę?