ggs_ppsd {ggmcmc}

Posterior predictive plot comparing the outcome standard deviation vs the distribution of the predicted posterior standard deviations.

Description

Histogram with the distribution of the predicted posterior standard deviations, compared with the standard deviations of the observed outcome.

Usage

ggs_ppsd(D, outcome, family = NA, bins = 30)

Arguments

D

Data frame whith the simulations. Notice that only the posterior outcomes are needed, and so either the ggs() call limits the parameters to the outcomes or the user provides a family of parameters to limit it.

outcome

vector (or matrix or array) containing the observed outcome variable. Currently only a vector is supported.

family

Name of the family of parameters to plot, as given by a character vector or a regular expression. A family of parameters is considered to be any group of parameters with the same name but different numerical value between square brackets (as beta[1], beta[2], etc).

bins

integer indicating the total number of bins in which to divide the histogram. Defaults to 30, which is the same as geom_histogram()

Value

A ggplot object.

References

Fernández-i-Marín, Xavier (2016) ggmcmc: Analysis of MCMC Samples and Bayesian Inference. Journal of Statistical Software, 70(9), 1-20. doi:10.18637/jss.v070.i09

Examples

data(linear)
ggs_ppsd(ggs(s.y.rep), outcome=y)