palaeoSig (version 2.1-3)

predict.agelme: Predicts the Calibrated age

Description

This function uses the output from 'agelme' to predict the Calibrated ages for specified depths.

Usage

# S3 method for agelme
predict(object, v = 1, depth,...)

Value

A list with three items

• v Whether constant variance or mu variance used.

• fit A data.frame of the predictions including 95% confidence interval.

  • DepthThe depths for the predicted ages

  • EstagePredicted age

  • LowlimLower 95% confidence interval

  • UpplimUpper 95% confidence interval

  • TsdTotal standard deviation

• data A data.frame containing the age and depth information of the radiocarbon dates.

Arguments

object

An 'agelme' model

v

Using constant (1) or mu (2) variance

depth

A vector of the depths to be predicted

...

Other arguments, currently unused.

Author

Einar Heegaard <[e-mail]

References

Heegaard, E., Birks, HJB. & Telford, RJ. 2005. Relationships between calibrated ages and depth in stratigraphical sequences: an estimation procedure by mixed-effect regression. The Holocene 15: 612-618

Examples

data(STOR)

fit.mod <- with(STOR,agelme(depthup,depthdo,cageup,cagedo))

#Predicting using the constant variance model,
#for each cm between 70 and 400 cm.
fit.pre <- predict(fit.mod,1,70:400)
plot(fit.pre)