dplyr model creation
do(mod = lm(num_purchases ~ amount, data=., na.action = na.exclude)) %>% tidy(mod) %>%