I like to estimate a mixed model with two Random Effects, that are independent of each other and among themselves.

I use Panel data with a nested structure (counties $j$ nested within regions $i$). The error component should decompose in the following way:

$$E_{ijt} = \alpha_i + \mu_{ij} + e_{ijt}$$

$\alpha_i$ is an unobservable region specific time-invariant effect (normal distributed)

$\mu_{ij}$ is a nested effect of county $j$ within region $i$ (normal distributed)

$e_{ijt}$ is a remainder disturbance (normal distributed)

I read somewhere that "lmer" is capable to fit such specifications with idenpendent Random Effects.

lmer(y~.... + (1|region) + (1|county)

Is the lmer-model adequate for my task?


1 Answer 1


Yes, lme4 should be fine in here. The formula you are asking is:

y ~ 1 + (1|region) + (1|region:county)

because in model you described country is nested in region.

For an extended description and examples see here, here or here.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.