R前向选择强制变量保持在等式中

Question

我正在使用755个观测值和16个变量进行逻辑回归。我正在使用glm函数进行变量选择。 glm找到了8个变量的最佳模型。我希望这些变量被强制留在并使用glm和step找到下一个最好的9变量模型（见下文）。我想这样做，直到我为9-16变量的模型进行了前向选择 （选择了所有16个变量）。

我的代码看起来像

饱和模型

full=glm(PREVAP ~ SEX + TOTCHOL + AGE + SYSBP + DIABP + as.factor(CURSMOKE)    
+ CIGPDAY + BMI + as.factor(DIABETES) + as.factor(BPMEDS) + HEARTRTE + 
  GLUCOSE + as.factor(EDUC) + TIME + HDLC + LDLC, data=training, 
  family=binomial(link="logit"))
summary(full)
anova(full,test="Chisq")
full.forward <- step(null,     
scope=list(lower=null,upper=full),direction="forward", 
      family=binomial(link="logit"))

这给了我一个8个因子的模型 我需要在下一个模型中强制使用这些因子，并使用前向选择找到9个因子的模型。怎么办？

我被告知bestglm和glmnet也允许这样但我不知道这些包。

你能帮忙吗？这些包有很多选择。

此致 Mary A. Marion