library(filinreg)
set.seed(666)
x <- rnorm(15)
y <- x + rnorm(15)
new <- data.frame(x = c(-3,3))
predict(lm(y ~ x), new, interval = "prediction")
fi <- filinreg(y ~ x, L = 100L, distr = "normal", lucky = TRUE)
fpred <- filinregPredictive(fi, new)
fiSummary(fpred)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.