put everything in one function
added some plots
added some lmer analysis