dalmatian (version 0.3.0)

convergence.dalmatian: Convergence

Description

Compute convergence diagnostics for a dalmatian object.

Usage

# S3 method for dalmatian
convergence(object, pars = NULL,
  nstart = start(object$coda), nend = end(object$coda),
  nthin = coda::thin(object$coda), raftery = NULL, ...)

Arguments

object

Object of class dalmatian created by dalmatian().

pars

List of parameters to assess. If NULL (default) then diagnostics are computed for the fixed effects and random effects standard deviations in both the mean and variance models.

nstart

Start point for computing summary statistics (relative to true start of chain).

nend

End point for computing summary statistics (relative to true start of chain).

nthin

Thinning factor for computing summary statsitics (relative to full chain and not previously thinned output).

raftery

List of arguments to be passed to raftery.diag(). Any values not provided will be set to their defaults (see help(raftery.diag()) for details).

...

Ignored

Value

List containing Gelman-Rubin and Raftery convergence diagnostics and effective sampel sizes for the selected parameters.

Examples

Run this code
# NOT RUN {
## Load output from previously run model
load(system.file("Pied_Flycatchers_1","pfresults.RData",package="dalmatian"))

## Compute convergence diagnostics
pfconvergence <- convergence(pfresults)

# }

Run the code above in your browser using DataCamp Workspace