I would like to use MixedModels.jl to fit regression models with random effects.
Looking at the official docs and some examples on the Internet I’ve found different syntaxes.
For example if I want to fit a Logit model (my outcome is binary) I can do it…:
fit(GeneralizedLinearMixedModel, @formula(y ~ 1 + x + (1 | id)), df, Bernoulli()) fit!(glmm(@formula(y ~ 1 + x + (1 | id)), df, Binomial(), LogitLink()))
What’s the difference between both syntaxes?