代码之家 › 专栏 › 技术社区 › Marishka Usacheva

数据帧中出错tmp替换有x数据有y

lm regression r

Marishka Usacheva · 技术社区 · 7 年前

# Create variables for child's EA:

dat$cldeacdi <- rowMeans(dat[,c('cdcresp', 'cdcinv')],na.rm=T)
dat$cldeacu <- rowMeans(dat[,c('cucresp', 'cucinv')],na.rm=T)

# Create a residual score for child EA:

dat$cldearesid <- resid(lm(cldeacu ~ cldeacdi, data = dat))

我收到以下消息:

Error in `$<-.data.frame`(`*tmp*`, cldearesid, value = c(-0.18608488908881,  : 
  replacement has 366 rows, data has 367

我搜索了此错误,但找不到任何可以解决此问题的方法。此外,我为mom的EA创建了完全相同的代码,它很好地保存了剩余的代码,没有任何错误。如果有人能帮我解决这个问题,我将不胜感激。

2 回复 | 直到 7 年前

LyzandeR 7 年前

我有种感觉你有 NA 在您的数据中。看看这个例子:

#mtcars data set
test <- mtcars
#adding just one NA in the cyl column
test[2, 2] <- NA

#running linear model and adding the residuals to the data.frame
test$residuals <- resid(lm(mpg ~ cyl, test))
Error in `$<-.data.frame`(`*tmp*`, "residuals", value = c(0.382245430809409,  : 
  replacement has 31 rows, data has 32

作为验证:

length(resid(lm(mpg ~ cyl, test)))
#31
nrow(test)
#32

lm 将运行 na.omit 因此,如果您有任何带有NA的行,这些行将被消除,从而导致更少的结果。

如果你跑步 不适用。省略 在您的 dat 数据集(即。 dat <- na.omit(dat) 在代码的最开始,您的代码应该可以工作。

Dawson A. Verley 5 年前

这是一条古老的线索,但也许这可以帮助其他人面对同样的问题。对于LyzandeR的观点,检查NA作为第一道防线。此外,确保没有任何因素 x

推荐文章

Shawn Hemelstrand · 如何在ggplot中的许多回归线中为一条特定回归线上色以匹配注释?

2 年前

Arvind Sharma · 在R中编写一个用于回归的循环,替换自变量进行稳健性检查

2 年前

krassowski · 如何将'VGAM::cumulative'包装到助手函数中('object not found'问题)?

3 年前

Nazanin · CNN模型的核尺寸

3 年前

Anna Carolina de RoldÃ£o · R中的错误消息:错误`[.data.frame`(m,labs):选择了未定义的列

3 年前

D. Smel · 如何在线性模型的两侧生成对数函数

6 年前

MTT · 最适合散点图的回归

6 年前

santobedi · scikit学习中的多输出高斯过程回归

7 年前

Ben · FELM+Stargazer-将工具变量估计与OLS对齐

7 年前

Anx8 · minepy:缓冲区的维度数错误

7 年前

数据帧中出错*tmp*替换有x数据有y

数据帧中出错tmp替换有x数据有y