getWeights {cAIC4}R Documentation

Optimize weights for model averaging.

Description

Function to constructed an optimal vector of weights for model averaging of Linear Mixed Models based on the proposal of Zhang et al. (2014) of using Stein's Formular to derive a suitable criterion based on the conditional Akaike Information Criterion as proposed by Greven and Kneib. The underlying optimization used is a customized version of the Augmented Lagrangian Method.

Usage

getWeights(models)

Arguments

models

An list object containing all considered candidate models fitted by lmer of the lme4-package or of class lme.

Value

An updated object containing a vector of weights for the underlying candidate models, value of the object given said weights as well as the time needed.

WARNINGS

For models called via gamm4 or gamm no weight determination via this function is currently possible.

Author(s)

Benjamin Saefken & Rene-Marcel Kruse

References

Greven, S. and Kneib T. (2010) On the behaviour of marginal and conditional AIC in linear mixed models. Biometrika 97(4), 773-789.

Zhang, X., Zou, G., & Liang, H. (2014). Model averaging and weight choice in linear mixed-effects models. Biometrika, 101(1), 205-218.

Nocedal, J., & Wright, S. (2006). Numerical optimization. Springer Science & Business Media.

See Also

lme4-package, lmer, getME

Examples

data(Orthodont, package = "nlme")
models <- list(
    model1 <- lmer(formula = distance ~ age + Sex + (1 | Subject) + age:Sex,
               data = Orthodont),
    model2 <- lmer(formula = distance ~ age + Sex + (1 | Subject),
               data = Orthodont),
    model3 <- lmer(formula = distance ~ age + (1 | Subject),
                 data = Orthodont),
    model4 <- lmer(formula = distance ~ Sex + (1 | Subject),
                data = Orthodont))

foo <- getWeights(models = models)
foo



[Package cAIC4 version 0.9 Index]