我正在尝试在R中设置更改分数分析。到目前为止，我的代码看起来像这样：
modChange <- lm(water81~income + water80, data = concord)
summary(modChange)
coef(modChange)
x2.y <- modChange$coef[3]-modChange$coef[1]
score <- lm(x2.y~ income + water80, data = concord)
summary(score)
我理解变化得分是从第一次回归得到的（x2-y），然后在术语x1和x2上回归。一旦进入第二次回归，我的代码就是抛出
“model.frame.default中的错误（公式= x2.y~收入+ water80，数据 = concord，：变量长度不同（找到'收入'）“
任何有关我是否设置错误或采取新方向的见解都将不胜感激。
R版.99.48，在Windows上运行。
谢谢！