如何在 R 中创建线性模型的预测区间并可视化

如何在 R 中创建线性模型的预测区间并可视化

线性回归是一种常用的统计方法，用于建立自变量（X）与因变量（Y）之间的线性关系模型。在进行线性回归分析时，除了对因变量的预测值感兴趣外，我们还希望能够了解预测的不确定性范围。通过创建线性模型的预测区间，我们可以获得对预测值的置信区间，从而评估模型的准确性和稳定性。本文将介绍如何在 R 语言中创建线性模型的预测区间，并通过可视化的方式展示结果。

首先，我们需要准备一些数据以进行线性回归分析。以下是一个示例数据集，其中包含自变量 X 和因变量 Y 的观测值：

# 创建示例数据
X <- c(1, 2, 3, 4, 5)
Y <- c(3, 5, 7, 9, 11)

# 将数据合并为数据框
data <- data.frame(X, Y)

接下来，我们可以使用 lm() 函数来拟合线性回归模型，并使用 predict() 函数进行预测。同时，我们可以使用 predict() 函数的 interval 参数来计算预测区间。以下是相应的代码：

# 拟合线性回归模型
model <- lm(Y ~ X, data = data)

# 进行预测，并计算预测区间
predictions <- predict(model, interval = "prediction")

在上述代码中，lm() 函数用于拟合线性回归模型，其中 Y ~ X 表示因变量 Y 与自变量 X 之间的线