如何解决R语言中使用geom_abline时出现的‘subscript out of bounds’错误？

核心问题：系数名称匹配错误 #

你遇到的subscript out of bounds错误，是因为提取线性模型系数时使用的名称和模型实际输出的系数名称不匹配：

1. 线性模型的截距项名称是**(Intercept)**（首字母大写），而你写的是(intercept)（全小写），导致找不到对应下标。
2. 自变量ln(MW)的系数名称是**ln(MW)**，但你用了带反引号的\ln(MW)``，同样无法匹配。

解决步骤 #

步骤1：确认模型系数的正确名称 #

先运行以下代码查看模型的系数名称：

names(coef(SEC.lm))

你会看到输出类似：

[1] "(Intercept)" "ln(MW)"

步骤2：修正geom_abline的参数 #

把代码里的系数提取部分改成正确的名称：

SEC.lm <- lm(`Elution time (sec)` ~ `ln(MW)`, data = calibration_data) # 建议加上data参数，避免环境变量问题

ggplot(data = calibration_data,
       aes(x = `ln(MW)`, 
           y = `Elution time (sec)`)) +
  geom_point(pch = 0) + 
  geom_abline(intercept = coef(SEC.lm)[["(Intercept)"]],
              slope = coef(SEC.lm)[["ln(MW)"]])

更简便的替代方案：用geom_smooth直接绘制回归线 #

不需要手动提取系数，ggplot的geom_smooth可以直接拟合并绘制线性回归线，代码更简洁且不易出错：

ggplot(data = calibration_data,
       aes(x = `ln(MW)`, 
           y = `Elution time (sec)`)) +
  geom_point(pch = 0) + 
  geom_smooth(method = "lm", se = FALSE) # se=FALSE去掉置信区间，需要的话可以保留

额外注意事项 #

• 确保calibration_data数据框中确实存在Elution time (sec)和ln(MW)这两个变量，且没有大量缺失值导致模型拟合失败。可以用summary(SEC.lm)查看模型拟合情况，确认系数是否正常生成。

内容的提问来源于stack exchange，提问作者Kayley Brownbill