values <- data.frame(x_1 = 0, x_2 = 1) predict(regmodel, values, interval = "confidence")
values <- data.frame(x_1 = 0, x_2 = 1) predict(regmodel, values, interval = "prediction")