[R] a question on stepAIC

From: Minya Pu <minya.pu_at_gmail.com>
Date: Thu, 6 Dec 2007 10:47:30 -0800


 Hi,
 I am trying to use stepAIC to do forward stepwise selection for a logistic model.
I did

 hi <- glm(melanoma ~ ., family = binomial, data=data)

lo <- glm(melanoma ~ 1, family = binomial, data=data)

 stepAIC(lo, scope=list(upper = hi , lower = lo), direction = "forward" )

But the full model hi did not converge since I have many many candidate covariates:

Warning message:

algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart,

 Question: Can I do the forward selection without specifying a maximal model hi?

Thanks,

Minya Pu

        [[alternative HTML version deleted]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 06 Dec 2007 - 18:50:55 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 06 Dec 2007 - 21:30:17 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.