What's the difference between (~1 +....) and (1 | ...) and (0 | ...) etc.?
Say you have variable V1 predicted by categorical variable V2, which is treated as a random effect, and continuous variable V3, which is treated as a linear fixed effect. Using lmer syntax, simplest model (M1) is:
V1 ~ (1|V2) + V3
This model will estimate:
P1: A global intercept
P2: Random effect intercepts for V2 (i.e. for each level of V2, that level's intercept's deviation from the global intercept)
P3: A single global estimate for the effect (slope) of V3
The next most complex model (M2) is:
V1 ~ (1|V2) + V3 + (0+V3|V2)
This model estimates all the parameters from M1, but will additionally estimate:
P4: The effect of V3 within each level of V2 (more specifically, the degree to which the V3 effect within a given level deviates from the global effect of V3), while enforcing a zero correlation between the intercept deviations and V3 effect deviations across levels of V2.
This latter restriction is relaxed in a final most complex model (M3):
V1 ~ (1+V3|V2) + V3
In which all parameters from M2 are estimated while allowing correlation between the intercept deviations and V3 effect deviations within levels of V2. Thus, in M3, an additional parameter is estimated:
P5: The correlation between intercept deviations and V3 deviations across levels of V2
Usually model pairs like M2 and M3 are computed then compared to evaluate the evidence for correlations between fixed effects (including the global intercept).
Now consider adding another fixed effect predictor, V4. The model:
V1 ~ (1+V3*V4|V2) + V3*V4
would estimate:
P1: A global intercept
P2: A single global estimate for the effect of V3
P3: A single global estimate for the effect of V4
P4: A single global estimate for the interaction between V3 and V4
P5: Deviations of the intercept from P1 in each level of V2
P6: Deviations of the V3 effect from P2 in each level of V2
P7: Deviations of the V4 effect from P3 in each level of V2
P8: Deviations of the V3-by-V4 interaction from P4 in each level of V2
P9 Correlation between P5 and P6 across levels of V2
P10 Correlation between P5 and P7 across levels of V2
P11 Correlation between P5 and P8 across levels of V2
P12 Correlation between P6 and P7 across levels of V2
P13 Correlation between P6 and P8 across levels of V2
P14 Correlation between P7 and P8 across levels of V2
Phew, That's a lot of parameters! And I didn't even bother to list the variance parameters estimated by the model. What's more, if you have a categorical variable with more than 2 levels that you want to model as a fixed effect, instead of a single effect for that variable you will always be estimating k-1 effects (where k is the number of levels), thereby exploding the number of parameters to be estimated by the model even further.
I'm not 100% sure, so check carefully, and you'll probably want to put priors on the variances instead of the 1.0E-12, but maybe something like this?
model {
for( i in 1 : nData ) {
crime[i] ~ dbern( mu[i] )
logit(mu[i]) <- base + b0[counties[i], period[i]] + b1[counties[i], period[i]] * police[i]
}
base ~ dnorm( 0 , 1.0E-12)
for (j in 1 : nCounties){
for (k in 1 : nPeriod) {
b0[j, k] ~ dnorm(b0h[j], 1.0E-12)
b1[j, k] ~ dnorm(b1h[j], 1.0E-12)
}
}
for ( j in 1 : nCounties ) {
b0h[j] ~ dnorm( 0 , 1.0E-12 )
b1h[j] ~ dnorm( 0 , 1.0E-12 )
}
}
Best Answer
Not enough reputation to comment, so I'll post this as an answer. There are a number of questions like this already around. you might want to look at this message.
However,
(1|group1/group2)should work with all but very old versions of lme4, so if that gives you an error, there is probably something wrong with the way you set up your data. Note that once your data are correctly set up,(1|group1/group2)and(1|group1) + (1|group2)should give the same results.