r - 将残差绑定(bind)到具有缺失值的输入数据集

Question

这个问题在这里已经有了答案:





How do I deal with NAs in residuals in a regression in R?

(3 个回答)


4年前关闭。




我正在寻找一种方法来绑定(bind)lm输入数据集的残差。方法必须加上NA对于缺失的残差，残差应该对应于正确的行。

样本数据:

N <- 100 
Nrep <- 5 
X <- runif(N, 0, 10) 
Y <- 6 + 2*X + rnorm(N, 0, 1) 
X[ sample(which(Y < 15), Nrep) ] <- NA
df <- data.frame(X,Y)

residuals(lm(Y ~ X,data=df,na.action=na.omit))

残差应该绑定(bind)到df。

Answer 1

只需更改 na.action至na.exclude :

residuals(lm(Y ~ X, data = df, na.action = na.exclude))

na.omit和 na.exclude两者都对预测变量和标准进行逐案删除。它们的区别仅在于提取器功能如 residuals()或 fitted()对于 na.exclude 省略的情况，将用 NA 填充它们的输出，因此具有与输入变量相同长度的输出。

(这是找到的最佳解决方案 here )