Sorry, I should have mentioned this, equatiomatic can’t express this model because of the lack of concordance in the intercept between population (fixed) and group (random) level effects. See here: lme4 models with random effects with suppressed intercepts, (0 + f | g), fail · Issue #153 · datalorax/equatiomatic · GitHub .