RandomEffectsExcess: Rnadom effects in excess

In oswaldosantos/INLAOutputs: Process selected outputs from the 'INLA' package

Description

Probability of random effects being higher than a cutoff value.

Usage

RandomEffectsExcess(..., cutoff = 0, rnd = 3)

Arguments

...	objects of class inla with random marginals.
cutoff	value above which random effects would be in excess. The cutoff is used for untransformed marginal posteriors, so cutoff = 0 (default) is used to calculate the probability of relative risk (odds ratio) > 1.
rnd	integer indicating the number of decimal places to be used.

Value

vector if only one model with one random effect is evaluated; data.frame or list otherwise.

References

Blangiardo, Marta, and Michela Cameletti. Spatial and Spatio-temporal Bayesian Models with R-INLA. John Wiley & Sons, 2015.

Examples

data(sp)

mod1 <- inla(aan ~ f(id, model = 'bym', graph = sp.adj),
             E = eaan,
             family = 'poisson', data = spn)

mod2 <- inla(aan ~ shvn + f(id, model = 'bym', graph = sp.adj),
             E = eaan,
             family = 'poisson', data = spn)

see <- RandomEffectsExcess(mod1, mod2)
summary(see)

oswaldosantos/INLAOutputs documentation built on May 24, 2019, 5 p.m.