I have various datasets I need to analyse regarding soil properties, all in the same fashion, with one fixed effect (which is a position along a transect, indicating different land uses). Now my main level of replication is across different transects, which will obviously have some form of random variance associated with them, and so I want to account for this in my statistical analyses.
lme4 I specified a mixed model to this specification
model <- lmer(variable.of.interest ~ transect.position + (1|transect))
Now, when I analyse the above model against a model without the transect position term, I get exactly the results I was expecting, and then plugging the above model into
anova(), I get the F values, d.f. etc. that I need.
However, I can't figure out how to say, for my overall report, that the random effect of transect does not make any difference to the overall analyses (i.e. I can't get a p-value, F value, d.f. etc.).