contrasts.mcemGLMM: Contrast estimation for mcemGLMM objects

In mcemGLM: Maximum Likelihood Estimation for Generalized Linear Mixed Models

Contrast estimation for mcemGLMM objects

Description

Contrast testing for a model fitted with mcemGLMM.

Usage

contrasts.mcemGLMM(object, ctr.mat)

Arguments

object	a model fitted with the mcemGLMM function.
ctr.mat	contrast matrix. Each row corresponds to a linear combination of the fixed effect coefficients.

Value

A matrix with each row corresponding to a contrast and with the following columns:

Estimate

contrasts's point estimator.

Std. Err.

standard error for each fitted contrast.

Wald

Wald statistic for each fitted contrast.

Adj. p-value

p-value adjusted for multiple comparison (Bonferroni.)

Author(s)

Felipe Acosta Archila <[e-mail]>

Examples

set.seed(72327)
data(exdata)
fit <- mcemGLMM(obs ~ z2 + x, random = ~ 0 + z1, 
                data = exdata, 
                family = "bernoulli", vcDist = "normal", 
                controlEM = list(verb = FALSE, EMit = 15, MCit = 10000), 
                initial = c(-0.13, -0.15, -0.21, 1.59, 0.002))
mat <- rbind("1 - 2" = c(0, -1, 0, 0), "1 - 3" = c(0, 0, -1, 0), "2 - 3" = c(0, 1, -1, 0))
contrasts.mcemGLMM(fit, mat)

mcemGLM documentation built on April 3, 2023, 5:43 p.m.