[R] step and glm : keep all models

[Prof Brian Ripley]

On Thu, 24 Jan 2008, Sylvain Bertrand wrote:

> Hi everyone,
>
> I'm running the following command:
>
>> step(mydata.glm,directions="both",trace=T)
>
> which returns the model with the lowest AIC.
> Now I'd like to get the list of all the models (or at least the formulas)
> that were used in between.
>
> I've noticed the "keep" option, but couldnt find any help on how to use it.
>
> Would anyone here be able to help me with this?

Did you try the example on the help page?  The final display shows you how 
to see the list of models.  Or

step(lm1, keep=function(model, aic) formula(model))$keep[1,]

(It's a bit peculiar: a formula is of length 3 but prints like length 1.)

-- 
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