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