calcBaseTemp: calcBaseTemp

In twDEMC: parallel DEMC

Description

calculate Base temperature from given temperatures of data streams and corresponding number of observations.

Usage

calcBaseTemp(temp, nObs, TFix = rep(NA_real_, length(temp)), 
    iNonFixTemp = which(!is.finite(TFix)))

Arguments

temp	numeric vector (nResultComp): temperature, i.e. variance inflation factor
nObs	integer vector (nResultComp): number of observations
TFix	numeric vector (nResultComp): fixed temperature for components, non-finite for those with varying temperature , alternatively a named vector listing only the components with fixed temperatures (temp must have names then too)
iNonFixTemp	integer vector: index of result components, which Temperature is fixed

Details

The Temperature, i.e. variance inflation factor, scales with the number of observations by T = 1 + (T0 - 1) * n. A common base temperature is the maximium across calculated stream base temperatures

Value

numeric scalar: The temperature, i.e. variance inflation factor, for a single observation.

Author(s)

Thomas Wutzler

See Also

calcStreamTemp, twDEMC

calcTemperatedLogDen.default

Examples

data(twdemcEx1)   
.nObs <- c(parms=getNParms(twdemcEx1), obs=length(twdemcEx1$dInfos[[1]]$argsFLogDen$obs) )
.T <- getCurrentTemp(twdemcEx1)
calcBaseTemp( .T, .nObs[names(.T)], TFix=c(parms=1) )

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