I have a multilevel model that has countries nested in geopolitical groupings. Let’s say the model is:

y ~ x + (1 | continent / country).

I have three set of intercepts,

A global intercept.

An intercept for ‘continent’

And and intercept for ‘continent:country’

I can calculate each separately using ranef().

Do I use coef(model)$‘continent:country’ to get the final intercept for each country (which does not seem to be the case),

or do I calculate use ranef() for each of the above and sum them up?

coef(model)$continent seems to give me global intercept + ranef(model)$continent, as explained by @paul.buerkner here.

but coef(model)$‘continent:country’ gives me something else.