Get variance components for binomial model with binary response.
var_comps_binary(mod, expct)
mod
expct
Fixed, random and residual variance
Useful links