将残差合并到缺少obs na.exclude的数据中不起作用

Question

我有以下命令来估计线性回归的残差，然后将其合并到数据中。数据在感兴趣的变量中有缺失值，但正如我在许多变量上做的那样，我希望 NA 值回归的残差只是 NA。

我在 text 上看到 na.action=na.exclude 应该可以很好地完成工作，但它似乎什么也没做。你知道为什么 na.action 不起作用吗？

d <- lm(log(br_prix) ~ sectoryear, data=test,na.action = na.exclude)

P1 <- as.data.frame(d$residuals)

P1 <- P1 %>% dplyr::rename(lnprice = "d$residuals")

test <- cbind(test,P1)