R语言model.frame.default报错：变量长度不一致问题求助

核心问题：变量名不匹配 #

你的报错信息里出现了变量x.bar，但你贴出的回归代码里用的是x2——这说明拟合模型时实际使用的变量名和你写的代码不一致。要么是你运行的代码和贴出来的版本存在手误（比如把x.bar写成了x2），要么是你的数据集里第二个自变量本来就叫x.bar，而非x2。

验证与修复步骤 #

1. 确认模型中的变量名
   运行以下代码查看拟合模型实际包含的变量：

   names(fit$model)
   

   如果输出里有x.bar而不是x2，需要统一变量名：

   • 回归代码改为：

     fit <- lm(y ~ x + x.bar + x:x.bar)
     

   • 预测代码对应改为：

     predict(fit, data.frame(x=-0.7, x.bar=1)) - predict(fit, data.frame(x=-0.5, x.bar=1))
     

2. 更高效的计算方式
   不用两次调用predict，可以直接用回归系数计算影响。你的模型形式为：
   [
   y = \beta_0 + \beta_1x + \beta_2x_2 + \beta_3x \cdot x_2
   ]
   当x2=1时，x从-0.5变到-0.7的影响可直接计算为：
   [
   (\beta_1 + \beta_3) \times (-0.7 - (-0.5))
   ]
   对应R代码：

   (coef(fit)["x"] + coef(fit)["x:x2"]) * (-0.7 - (-0.5))
   

   （如果变量名是x.bar，将代码里的x:x2替换为x:x.bar即可）

内容的提问来源于stack exchange，提问作者Sairah S