diag_mc_stderr: Batch Means Standard Errors

In njtierney/dsmcmc: tidy mcmc.list using data.table

Description

Calculate Monte Carlo Standard Errors Using Batch Means

Usage

diag_mc_stderr(x, b_size = "sqroot", warn = TRUE)

Arguments

x	'mcmc.list' or 'data.table' object
b_size	character, method for determining size of batch (see Details)
warn	logical, give a warning of there are too few samples in the MCMC output

Details

For the batch size, the options are 'sqroot', for using the using the square root of the MCMC output length, or 'cuberoot', for using the cube root of the output length. The default is 'sqroot'.

References

Galin L. Jones, Murali Haran, Brian S. Caffo, and Ronald Neath (2006). "Fixed-Width Output Analysis for Markov Chain Monte Carlo," Journal of the American Statistical Association, 101, 1537–1547

Examples

library(coda)
data(line)
r <- diag_mc_stderr(line)
r

njtierney/dsmcmc documentation built on Oct. 4, 2021, 10:17 p.m.