Question

我正在尝试使用 R 中的逐步 AIC 进行前向变量选择，但我认为我没有得到想要的结果。具体来说，该函数应该从没有变量开始，并不断添加变量并获取它们的 AIC 值。但是，当我运行它时，我只得到所有变量的 AIC 值。我哪里错了？这是我的代码-


model.full <- lm(distance ~ ., data = FAA_unique_without_speed_air)
model.null<-lm(distance ~ 1,  data = FAA_unique_without_speed_air)
modAIC <- MASS::stepAIC(model.full,direction='forward', scope=model.full, k = 2)

输出 -

Start:  AIC=9161.49
distance ~ aircraft + duration + no_pasg + speed_ground + height + 
    pitch

Answer 1

我认为最好明确使用 stepAIC 的参数，而不是使用默认值。试试：

当您想要进行前向选择时，提供空模型作为初始模型 object。
在 scope 中同时提供一个上下搜索公式。

例如，使用基础库 iris 中的 datasets 数据框：

library(MASS)

model.full <- lm(Sepal.Width ~ ., data = iris)
model.null <- lm(Sepal.Width ~ 1, data = iris)

MASS::stepAIC(model.null, direction = "forward", scope = list(lower = model.null,
                                                              upper = model.full))

或者如果你想使用默认值，那么你应该明确模型中包含的默认上层组件：

stepAIC(model.null, direction = "forward", scope = ~ Sepal.Length + Species + Petal.Length)

但是，正如@BenBolker 所提到的，您应该使用您的数据发布一个可重现的示例，以便我们进行确认。

R中使用前向选择的逐步AIC

1 个答案: