Bclim (version 2.2)

BclimInterp: Interpolation function for producing palaeoclimate histories on a grid

Description

The function takes output from previous stages of Bclim and produces interpolated values of palaeoclimate (in 3 dimensions) and palaeoclimate volatility.

Usage

BclimInterp(Bclim.data, Bclim.res, tgrid = seq(0, 14, by = 0.1))

Arguments

Bclim.data

Output from the function BclimMixSer or BclimMixPar

Bclim.res

Output from the function BclimMCMC

tgrid

A grid (not necessarily regular) on which to interpolate climate

Value

• Returns a list containing the following elements:
• clim.interp
  Interpolated values of 3D climate
• vol.interp
  Interpolated values of 3D climate volatility
• time.grid
  The time grid used for interpolation

Details

A Bclim run can be separated into 4 stages: 1 - Turning each individual pollen layer into climate marginal data posteriors, 2 - Approximating the marginal data posteriors as mixtures of Gaussians, 3 - Running an MCMC stage to estimate climates and climate volatilities, 4 - Interpolating the climates and volatilities onto a regular grid. This function runs stage 4.

References

See Arxiv paper at http://arxiv.org/abs/1206.5009.

See Also

The main Bclim function is BclimRun. See also the other 3 stages: BclimLayer, BclimMixSer (or BclimMixPar), and BclimMCMC,

Examples

## See help for BclimMCMC and webpage at http://mathsci.ucd.ie/~parnell_a/Bclim.html