Mixed model parameterization

Dear list members,

I am measuring an index in plots, in different hours, different days, and
different months in 3 different "habitats".

12 plots by "habitats".

I would like to estimate the mean index in the 3 environments.

Is this model appropriate to achieve my goal?

model <- lmer(index ~  habitat + (1 | plot/hour/day/month), data = mydata,
REML = FALSE)

Best,

Manuel

*Manuel Sp?nola, Ph.D.*
Instituto Internacional en Conservaci?n y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola at una.cr <mspinola at una.ac.cr>
mspinola10 at gmail.com
Tel?fono: (506) 8706 - 4662
Personal website: Lobito de r?o <https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

	[[alternative HTML version deleted]]