bhpm.summary.stats {bhpm}

Summary Statistics for the Posterior Distributions in the model.

Description

Returns the Summary Statistics for the Posterior Distributions in the model.

Usage

	bhpm.summary.stats(raw, prob)

Arguments

raw

The output from a model simulation (e.g. bhpm.cluster.BB.hier3).

prob

The probability level for HPI. Default is 0.95.

Details

The function reports the mean, upper and lower bounds of the HPI (highest probabily interval), the standard deviation and MCMC standard error.

Value

Returns a list of the summary statistics for each sampled variable. Each element of the list is a data.frame containing at least the columns mean, median, hpi_lower, hpi_upper, SD and SE. Columns which may be used to indentify the individual variables are Outcome.Grp, Outcome, and Cluster.

Note

The MCMC error is found using the 'coda' summary function.

Author(s)

R. Carragher

Examples

data(bhpm.cluster.data1)
data <- subset(bhpm.cluster.data1, Cluster == '0.0-180.0')
raw = bhpm.npm(data, burnin = 100, iter = 200)
summ = bhpm.summary.stats(raw)


data(bhpm.cluster.data1)
raw = bhpm.npm(bhpm.cluster.data1)
summ = bhpm.summary.stats(raw)