[R] ANother s-plus to R problem

[Thomas Lumley]

On Fri, 5 Dec 2003, Prof Brian Ripley wrote:

> On Fri, 5 Dec 2003, Ognen Duzlevski wrote:
>
> > Hi all, thank you for replying so quickly!
> >
> > I have another problem:
> >
> >  step.wise <- stepwise(data.table[,names(good.motifs[,1])], data.table[1],
> > f.crit=fval.cutoff)
> >
> > s-plus has the stepwise() formula but R has step() and stepAIC() from base
> > and MASS packages. I cannot seem to figure out how to convert the above
> > stepwise to either step() or stepAIC().
>
> You can't.  This is an old-fashioned approach, and the closest equivalent
> in R is probably that of regsubsets in package leaps (in one of its
> stepwise modes: it is a little short of detail)
>

The default method for stepwise() in S-PLUS seems to be "efroymson". This
is not provided by regsubsets(), but it is in the Fortran code, so it
could be added.

	-thomas