diagnostic: Calculation of a proposal for burn-in phase and thinning rate

Description Usage Arguments Value

Description

The proposed burn-in is calculated by dividing the Markov chains into m blocks and calculate the 95% credibility intervals and the respective mean. Starting in the first one, the block is taken as burn-in as long as the mean of the current block is not in the credibility interval of the following block or vice versa. The thinning rate is proposed by the first lag which leads to a chain autocorrelation less than dependence. It is not easy to automate these choices, so it is highly recommended to verify the chains manually.

Usage

1
diagnostic(chain, dependence = 0.8, m = 10)

Arguments

chain

vector of Markov chain samples

dependence

allowed dependence for the chain

m

number of blocks

Value

vector of burn-in and thinning


BaPreStoPro documentation built on May 2, 2019, 3:34 p.m.