contrasts.mcemGLMM {mcemGLM}

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)