lm.full <- lm(abund ~ . - Present, data=Dickcissel) lm.step <- step(lm.full) summary(lm.step)