statnet.common (version 4.9.0)

mcmc-utilities: Utility operations for mcmc.list objects

Description

colMeans.mcmc.list is a "method" for (non-generic) colMeans applicable to mcmc.list objects.

sweep.mcmc.list is a "method" for (non-generic) sweep applicable to mcmc.list objects.

lapply.mcmc.list is a "method" for (non-generic) lapply applicable to mcmc.list objects.

Usage

colMeans.mcmc.list(x, ...)
sweep.mcmc.list(x, STATS, FUN = "-", check.margin = TRUE, ...)
lapply.mcmc.list(X, FUN, ...)

Value

colMeans.mcmc returns a vector with length equal to the number of mcmc chains in x with the mean value for each chain.

sweep.mcmc.list returns an appropriately modified version of x

lapply.mcmc.list returns an mcmc.list each of whose chains had been passed through FUN.

Arguments

x

a mcmc.list object.

...

additional arguments to colMeans or sweep.

STATS, FUN, check.margin

See help for sweep.

X

An mcmc.list object.

See Also

colMeans, mcmc.list

sweep

lapply

Examples

data(line, package="coda")
summary(line) # coda
colMeans.mcmc.list(line) # "Method"
# \dontshow{
stopifnot(isTRUE(all.equal(summary(line)$statistics[,"Mean"],colMeans.mcmc.list(line))))
# }
data(line, package="coda")
colMeans.mcmc.list(line)-1:3
colMeans.mcmc.list(sweep.mcmc.list(line, 1:3))
# \dontshow{
stopifnot(isTRUE(all.equal(colMeans.mcmc.list(sweep.mcmc.list(line, 1:3)), colMeans.mcmc.list(line)-1:3)))
# }
data(line, package="coda")
colMeans.mcmc.list(line)[c(2,3,1)]
colMeans.mcmc.list(lapply.mcmc.list(line, `[`,,c(2,3,1)))
# \dontshow{
stopifnot(isTRUE(all.equal(colMeans.mcmc.list(line)[c(2,3,1)],colMeans.mcmc.list(lapply.mcmc.list(line, `[`,,c(2,3,1))))))
# }