I'd love a check if anyone is willing!
I am trying to see if there is a statistical difference in female size between sites. Over the years females were repeatedly sampled within sites. I have sampled females opportunistically. Meaning that females were sampled a different number of times between and within sites.
My formula is:
> lmerfit1<-lmer(size ~ (1|FEMALE), data=Data) > lmerfit2<-lmer(size ~ SITE+(1|FEMALE), data=Data) > anova(lmerfit1, lmerfit2) Data: Data Models: lmerfit1: size ~ (1 | FEMALE) lmerfit2: size ~ SITE + (1 | FEMALE) Df AIC BIC logLik Chisq Chi Df Pr(>Chisq) lmerfit1 3 2167.8 2179.6 -1080.9 lmerfit2 4 2169.8 2185.5 -1080.9 0 1 **1**
A p value of 1 leaves me concerned. The other female traits I ran thru this same formula made sense.