rmaModel<-lmer(DIVERSITY~YEAR+(1|CIRCLE),data=data)

anova(rmaModel)

Type III Analysis of Variance Table with Satterthwaite's method

Sum Sq Mean Sq NumDF DenDF F value Pr(>F)

YEAR 4.4674 4.4674 1 672.72 31.328 3.174e-08 ***

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

I'm just not sure if this is what I should be doing, if there might be a better way and how to interpret my results correctly. Also, what is the best way to visualize this in R (aka graphs)?