数据子集的R回归的残差

Question

我想要获得回归残差，但仅从数据中获得残差子集：

我的R代码：

reg = lm(Y ~ X1+X2+.....+Xn,data=fic)
step_reg = step(reg, direction= "both")
summary(step_reg)

fic是一个有n列的数据帧，名为X1，X2，...Xn。

要获取所有残差：step_reg2$residuals

但我只想得到符合条件的行的残差，例如X1 = 'xxxx'，请问有什么解决方案？

Answer 1

您可以使用用于回归的数据来子集残差，如下所示：

reg <- lm(Sepal.Length ~ Sepal.Width + Petal.Length + Petal.Width, data=iris)
step_reg <- step(reg, direction= "both")
step_reg$residuals[iris$Species=="setosa"]

如果有缺失值：

x <- iris
x[1,2] <- NA
reg <- lm(Sepal.Length ~ Sepal.Width + Petal.Length + Petal.Width, data=x)
reg$residuals[iris[names(reg$residuals), "Species"] == "setosa"]