lm.mult <- lm(abund ~ clTma + NDVI + grass, data=Dickcissel) summary(lm.mult)