[R] ANother s-plus to R problem

[Prof Brian Ripley]

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)

I can't help asking: If you have S-PLUS (sic) code, why not use S-PLUS?

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595