Question

我正在尝试在104个变量中找到数量有限的预测变量（最大值= 6）。因此，我正在使用逐步回归（对于每个变量，我有10532个值）。我尝试过MATLAB：

mdl = stepwiselm(Pr, obs,'PEnter', 0.06)

但是，它给了我大约70个变量。另外，我尝试使用R软件包leaps解决问题：

b <- leaps::regsubsets(obs ~ ., data=Pr, nbest=1, nvmax=6)

我收到以下错误：

“ jumps.exhaustive（a，true.big）中的错误：详尽搜索将为S L O W，必须指定true.big = T“

我知道这应该是解决此问题的简便方法，但是我似乎无法弄清楚正确的格式。预先谢谢你。

Answer 1

使用

leaps::regsubsets(obs ~ ., data=Pr, nbest=1, nvmax=6, really.big=T)

或者您可以尝试

library(MASS)
# Fit the full model 
full.model <- lm(obs ~ ., data=Pr)
# Stepwise regression model
step.model <- stepAIC(full.model, direction = "both", 
                      trace = FALSE)
summary(step.model)

使用逐步回归找到指定数量的预测变量

1 个答案: