thin.twDEMC: thin twDEMC

Description Usage Arguments Author(s) See Also Examples

Description

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

Usage

1
2
3
## 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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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 31, 2017, 3:44 a.m.