summary.epimcmc: Summary method for 'epimcmc' objects

In EpiILM: Spatial and Network Based Individual Level Models for Epidemics

Description

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

Usage

## S3 method for class '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

## 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)

EpiILM documentation built on Jan. 13, 2021, 1:07 p.m.