View source: R/contrasts.mcemGLMM.R
contrasts.mcemGLMM | R Documentation |
Contrast testing for a model fitted with mcemGLMM
.
contrasts.mcemGLMM(object, ctr.mat)
object |
a model fitted with the mcemGLMM function. |
ctr.mat |
contrast matrix. Each row corresponds to a linear combination of the fixed effect coefficients. |
A matrix with each row corresponding to a contrast and with the following columns:
contrasts's point estimator.
standard error for each fitted contrast.
Wald statistic for each fitted contrast.
p-value adjusted for multiple comparison (Bonferroni.)
Felipe Acosta Archila <acosta@umn.edu>
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.