twDEMCBlock.twDEMCPops: twDEMCBlock twDEMCPops
In twDEMC: parallel DEMC

Description Usage Arguments Details Author(s) Examples

initialize twDEMCBlockInt by former run and append results to former run

## S3 method for class 'twDEMCPops'
twDEMCBlock(x,
 ..., TEnd = numeric(0), TEnd0, TStart = numeric(0), 
    TStart0, TFix, nObs, doRecordProposals = FALSE, extendRun = TRUE)

`x`	list of class twDEMCPops, result of `twDEMCBlockInt`
`...`	further arguments to `twDEMCBlockInt`
`TEnd`	numeric vector (nResComp) of end temperatures for each result component, if not given stays at current temperature, if in addition TSpec is specified, it has no effect
`TEnd0`	numeric scalar: alternative way of specifying end tempertaure: by baseTemperature for all data streams
`TStart`	numeric vector (nResComp): specifing starting temperature, if neigher TSpec, TStart, or TStart0 are given stays at the current temperature
`TStart0`	numeric scalar: specifing starting temperature by baseTemperature for all data streams
`TFix`	see `twDEMCSACont`, argument `ctrlT$TFix`. Must be given if specifying TEnd0 or TStart0
`nObs`	see `twDEMCSACont`, argument must be given if specifying TEnd0 or TStart0
`doRecordProposals`	if TRUE then an array of each proposal together with the results of fLogDen are recorded and returned in component Y
`extendRun`	if set to FALSE then only the new samples are returned

pops, dInfos, and blocks are reused from x or overwritten by arguments

TSpec: If TSpec is not explicitly given as an array for twDEMCBlockInt it is constructed as follows: , If TEnd is given use it to explicitely define end Temperature for all data streams. , If TEnd0 is given calculate TEnd for all data streams by scaling for number of observations (see calcStreamTemp). , If TStart0 is given calcualte TCurr for all data streams , else use current Temperature of streams in x as starting temperature

Thomas Wutzler

data(twdemcEx1) 		# previous run of twDEMCBlock
class(twdemcEx1)
twdemcEx1$thin			# thinning interval
(nGen0 <- getNGen(twdemcEx1))		# number of generations

# extend by 16 generations
nGen <- 16
#mtrace(twDEMCBlock.twDEMCPops)
res <- twDEMCBlock( twdemcEx1, nGen=nGen )

identical( nGen0+nGen, getNGen(res) )
plot( as.mcmc.list(res), smooth=FALSE )