weights argument in
mixed_model(). These are simple multipliers on the log-likelihood contributions of each group/cluster, i.e., we presume that there are multiple replicates of each group/cluster denoted by the weights.
The internal implementation of the
negative.binomial() family is now considerably faster.
More stable and accurate calculations using the matrixStats package.
Mixed models for ordinal clustered outcomes using the continuation ratio model, with forward and backward formulation.
The new function
VIF() calculates variance inflation factors for mixed models fitted in the package.
The CRAN version of the package has now only three basic vignettes to conform with the time restrictions of CRAN checks. The whole list of available vignettes can be found in the website of the package https://drizopoulos.github.io/GLMMadaptive/.
Corrected a typo in
anova.MixMod() that was reporting the AIC as the BIC of the fitted model.
The offset is now passed in the calculation of initial values.
Added support for the DHARMa package.
The new vignette Goodness of Fit for MixMod Objects describes how to check the fit of mixed models fitted by GLMMadaptive.
Added support for the effects package.
There is a new section in the vignette Methods for MixMod Objects illustrating the use of the effects package.
marginal_coefs() has a faster implementation. Compared to the previous implementation the results will be slightly different.
The new vignette Optimization and Numerical Integration in GLMMadaptive describes how to control the optimization and numerical integration procedures in the package.
predict() method now works for zero-inflated and hurdle models.
Hurdle Beta mixed effects models are now available using the
hurdle.beta.fam family object.
The new function
scoring_rules() calculates proper scoring rule for subject-specific predictions from mixed models for categorical data.
Added support for the emmeans package.
Hurdle Poisson and negative binomial models are now implemented using the family objects
A new vignette illustrating multiple comparisons with the multcomp package.
simulate() have gained the logical argument
sandwich to invoke the use of robust/sandwich standard errors in the calculations.
Zero-inflated Poisson and negative binomial models are now implemented using the family objects
zi.negative.binomial(), respectively. In addition, taking into advantage of the fact that users can specify their own log density functions for the outcome, two-part / hurdle model can also be implemented.
A new vignette illustrates how the zero-inflated models can be fitted.
“mean subject”: only the fixed effects part corresponding to predictions for the average subject (but not population averaged predictions in case of nonlinear link functions).
“marginal”: predictions using the marginalized coefficients that correspond to population averaged predictions.
“subject specific”: predictions at the subject level. These can be also calculated for subjects not originally in the dataset (i.e., estimates of the random effects are internally obtained).
simulate() method is available to simulate data from fitted mixed models. This can be used for instance to perform replication / posterior predictive checks.