Random effects error for multigroup analysis

[CReemts]

I am running a four-model path analysis that includes lmer and glmmTMB models with random effects. When I try to do a multigroup analysis, I get this error: "No random effects terms specified in formula". Is the multigroup analysis not possible with mixed models?

Reproducible example (gives the same error whether one or both of the models include random effects):
data(meadows)
meadows$random<-as.factor(rep(1:6))

jutila <- psem(
lmer(rich ~ elev + mass + (1|random), data = meadows),
lmer(mass ~ elev + (1|random), data = meadows)
)

jutila.multigroup <- multigroup(jutila, group = "grazed")