# fit model
lmod1 <- lm(murder ~ hs_grad + urban + poverty + single,
data = crime2009)
#fit without DC
lmod2 <- lm(murder ~ hs_grad + urban + poverty + single,
data = crime2009[-9, ])
#compare coefficients of models
coef_compare(lmod1, lmod2)
Run the code above in your browser using DataLab