Fits generalized linear mixed models for a single grouping factor under
    maximum likelihood approximating the integrals over the random effects with an 
    adaptive Gaussian quadrature rule; Jose C. Pinheiro and Douglas M. Bates (1995) 
    <doi:10.1080/10618600.1995.10474663>.

Allows users to fit a cosinor model using the 'glmmTMB' framework. 
    This extends on existing cosinor modelling packages, including 'cosinor' 
    and 'circacompare', by including a wide range of available link functions 
    and the capability to fit mixed models. The cosinor model is described by 
    Cornelissen (2014) <doi:10.1186/1742-4682-11-16>.

Approximate frequentist inference for generalized linear mixed model analysis with expectation propagation used to circumvent the need for multivariate integration. In this version, the random effects can be any reasonable dimension. However, only probit mixed models with one level of nesting are supported. The methodology is described in Hall, Johnstone, Ormerod, Wand and Yu (2018) <arXiv:1805.08423v1>.

Implements Bayesian spatial and spatiotemporal
    models that optionally allow for extreme spatial deviations through
    time. 'glmmfields' uses a predictive process approach with random
    fields implemented through a multivariate-t distribution instead of
    the usual multivariate normal.  Sampling is conducted with 'Stan'.
    References: Anderson and Ward (2019) <doi:10.1002/ecy.2403>.

Fit Generalized Linear Models to continuous and count outcomes, as well as estimate the prevalence of misrepresentation of an important binary predictor. Misrepresentation typically arises when there is an incentive for the binary factor to be misclassified in one direction (e.g., in insurance settings where policy holders may purposely deny a risk status in order to lower the insurance premium). This is accomplished by treating a subset of the response variable as resulting from a mixture distribution. Model parameters are estimated via the Expectation Maximization algorithm and standard errors of the estimates are obtained from closed forms of the Observed Fisher Information. For an introduction to the models and the misrepresentation framework, see Xia et. al., (2023) <https://variancejournal.org/article/73151-maximum-likelihood-approaches-to-misrepresentation-models-in-glm-ratemaking-model-comparisons>.