EpiILM (version 1.5.2)

summary.epimcmc: Summary method for epimcmc objects

Description

Summarize a epimcmc object and return an object of class summary.epimcmc.

Usage

# S3 method for epimcmc
summary(object, …)

Arguments

object

an S3 object of class epimcmc (i.e. the output of the epimcmc function).

…

potential further arguments (require by generic).

See Also

epimcmc, plot.epimcmc.

Examples

# NOT RUN {

## Example:  spatial SI model
# generate 100 individuals

x <- runif(100, 0, 10)

y <- runif(100, 0, 10)

covariate <- runif(100, 0, 2)

out1 <- epidata(type = "SI", n = 100, Sformula = ~covariate, tmax = 15,
               sus.par = c(0.1, 0.3), beta = 5.0, x = x, y = y)

alphapar1 <- matrix(c(1, 1, 1, 1), ncol = 2, nrow = 2)

betapar1 <- c(1, 1)

epi <- epimcmc(object = out1, tmin = 1, tmax = 15,
               niter = 1000, sus.par.ini = c(1, 1), beta.ini = 1,
               Sformula = ~covariate, pro.sus.var = c(0.5, 0.3), pro.beta.var = 0.1,
               prior.sus.dist = c("gamma", "gamma"), prior.beta.dist = "gamma",
               prior.sus.par = alphapar1, prior.beta.par = betapar1,
               adapt = TRUE, acc.rate = 0.5)

# summary of mcmc output
summary(epi)

# }