R 语言经典实例（原书第 2 版）

线性回归和方差分析

｜

341

11.15 绘制回归残差

11.15.1 问题

需要直观地显示回归残差。

11.15.2 解决方案

你可以使用 broom 包将模型结果放入数据框中，然后使用 ggplot 绘图：

m <- lm(y ~ x1 + x2)

library(broom)

augmented_m <- augment(m)

ggplot(augmented_m, aes(x = .fitted, y = .resid)) +

geom_point()

使用先前方法中的线性模型 m，我们可以创建一个简单的残差图：

library(broom)

augmented_m <- augment(m)

ggplot(augmented_m, aes(x = .fitted, y = .resid)) +

geom_point()

输出如图 11-5 所示。

.fitted

.resid

30

20

10

0

-

10

-

20

-

100 100 2000

图 11-5：模型残差图

你还可以使用基础 R 的 plot 方法快速查看，但它将生成基础 R 图形输出，而不是

ggplot 图形：