states <- as.data.frame(state.x77[,c("Murder", "Population","Illiteracy", "Income", "Frost")])
 fit <- lm(Murder ~ Population + Illiteracy + Income + Frost, data=states)
 stepAIC(fit, direction="backward")