Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

batchSE

Batch Standard Error


Description

Effective standard deviation of population to produce the correct standard errors.

Usage

batchSE(x, batchSize=100)

Arguments

x

An mcmc or mcmc.list object.

batchSize

Number of observations to include in each batch.

Details

Because of the autocorrelation, the usual method of taking var(x)/n overstates the precision of the estimate. This method works around the problem by looking at the means of batches of the parameter. If the batch size is large enough, the batch means should be approximately uncorrelated and the normal formula for computing the standard error should work.

The batch standard error procedure is usually thought to be not as accurate as the time series methods used in summary and effectiveSize. It is included here for completeness.

Value

A vector giving the standard error for each column of x.

Author(s)

Russell Almond

References

Roberts, GO (1996) Markov chain concepts related to sampling algorithms, in Gilks, WR, Richardson, S and Spiegelhalter, DJ, Markov Chain Monte Carlo in Practice, Chapman and Hall, 45-58.

See Also


coda

Output Analysis and Diagnostics for MCMC

v0.19-4
GPL (>= 2)
Authors
Martyn Plummer [aut, cre, trl], Nicky Best [aut], Kate Cowles [aut], Karen Vines [aut], Deepayan Sarkar [aut], Douglas Bates [aut], Russell Almond [aut], Arni Magnusson [aut]
Initial release
2020-09-30

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.