5

我正在尝试将 lm 的残差放回原始 data.frame 中:

fit <- lm(y ~ x, data = mydata, weight = ind)
mydata$resid <- fit$resid

如果残差的长度与 mydata 的行数相同,则第二行通常会起作用。但是,就我而言, ind 的一些元素是NA. 因此,剩余长度通常小于行数。也是fit$resid一个“数字”向量,所以我没有标签可以将它与 mydata data.frame 合并回来。有没有一种优雅的方式来实现这一目标?

4

1 回答 1

7

ind我认为如果只是一个向量 应该很容易。

sel <- which(!is.na(ind))
mydata$resid <- NA
mydata$resid[sel] <- fit$resid
于 2013-06-12T17:03:36.010 回答