I have a model such as this model <- lmer(minutes_after_sunrise ~ forest * forest_strata * habitat + (1 | forest_strata), data = data)