thin.twDEMC: thin twDEMC

In twDEMC: parallel DEMC

Description

Reduces the rows of an twDEMC object (list returned by twDEMCBlockInt) to correspond to a thinning of newThin.

Usage

## S3 method for class 'twDEMC'
thin(x,
 newThin = x$thin, start = 0, end = NA, ..., doKeepBatchCall = FALSE)

Arguments

x	the twDEMC list to thin
newThin	finite numeric scalar: the target thinning factor, must be positive multiple of x$thin
start	numeric scalar: the start time of the chain. Note that time starts from zero. If a vector or matrix is supplied (e.g. nGenBurnin) then the maximum is used
end	numeric scalar: the maximum end time of the chains. Note that time starts from zero. If a vector or matrix is supplied (e.g. nGenBurnin) then the maximum is used
...
doKeepBatchCall	wheter to retain the batch call attribute of x

Author(s)

Thomas Wutzler

See Also

subChains.twDEMC

Examples

data(twdemcEx1)
x <- twdemcEx1
c( nGen=getNGen(twdemcEx1), thin=twdemcEx1$thin, nSample=getNSamples(twdemcEx1) )

    .nGenBurnin <- max(getNGen(twdemcEx1))-70
thinned <- thin(twdemcEx1, start=.nGenBurnin)	# removing burnin period
c( nGen=getNGen(thinned), thin=thinned$thin, nSample=getNSamples(thinned), nGenBurnin=thinned$nGenBurnin )	#18 sample describing 68 generations

thinned <- thin(twdemcEx1, start=.nGenBurnin, newThin=8)	
c( nGen=getNGen(thinned), thin=thinned$thin, nSample=getNSamples(thinned), nGenBurnin=thinned$nGenBurnin )	#9 samples describing 64 generations

twDEMC documentation built on May 2, 2019, 5:38 p.m.