mc_anova_II: ANOVA type II table for mcglm objects via Wald test.

View source: R/mc_anova_II.R

mc_anova_IIR Documentation

ANOVA type II table for mcglm objects via Wald test.

Description

Performs Wald tests to generate type-II analysis-of- variance tables per response for model objects produced by mcglm.

Usage

mc_anova_II(object, verbose = TRUE)

Arguments

object

An object of mcglm class.

verbose

a logical if TRUE print some information about the tests performed. Default verbose = TRUE.

Value

Type II ANOVA table for mcglm objects.

Author(s)

Lineu Alberto Cavazani de Freitas, lineuacf@gmail.com

See Also

mc_anova_I, mc_anova_III and mc_anova_disp.

Examples


library(mcglm)
library(Matrix)
library(htmcglm)

data("soya", package = "mcglm")

form.grain <- grain ~ water * pot
form.seed <- seeds ~ water * pot

soya$viablepeasP <- soya$viablepeas / soya$totalpeas
form.peas <- viablepeasP ~ water * pot

Z0 <- mc_id(soya)
Z1 <- mc_mixed(~0 + factor(block), data = soya)

fit_joint <- mcglm(linear_pred = c(form.grain,
                                   form.seed,
                                   form.peas),
                   matrix_pred = list(c(Z0, Z1),
                                      c(Z0, Z1),
                                      c(Z0, Z1)),
                   link = c("identity",
                            "log",
                            "logit"),
                   variance = c("constant",
                                "tweedie",
                                "binomialP"),
                   Ntrial = list(NULL,
                                 NULL,
                                 soya$totalpeas),
                   power_fixed = c(TRUE,TRUE,TRUE),
                   data = soya)

mc_anova_II(fit_joint)


htmcglm documentation built on July 21, 2022, 5:10 p.m.