Given a model estimated using vglmer, this function
performs marginally augmented variational Bayes (MAVB) to improve the
approximation quality.
MAVB(object, samples, verbose = FALSE, var_px = Inf)This function returns a matrix with samples rows and columns
for each fixed and random effect.
Model fit using vglmer.
Number of samples to draw.
Show progress in drawing samples.
Variance of working prior for marginal augmentation. Default
(Inf) is a flat, improper, prior.
This function returns the improved estimates of the parameters. To use MAVB when generating predictions, one should use predict_MAVB. At present, MAVB is only enabled for binomial models.
Goplerud, Max. 2022a. "Fast and Accurate Estimation of Non-Nested Binomial Hierarchical Models Using Variational Inference." Bayesian Analysis. 17(2): 623-650.