## Imbrie and Kipp example
## load the example data
data(ImbrieKipp)
data(SumSST)
data(V12.122)
## merge training and test set on columns
dat <- join(ImbrieKipp, V12.122, verbose = TRUE)
## extract the merged data sets and convert to proportions
ImbrieKipp <- dat[[1]] / 100
V12.122 <- dat[[2]] / 100
## Fit a MAT model
(ik.mat <- mat(ImbrieKipp, SumSST, method = "chord"))
## Reconstruct pH for the RLGH core
v12.pH <- predict(ik.mat, V12.122)
## draw the reconstruction
reconPlot(v12.pH, use.labels = TRUE, display.error = "bars",
xlab = "Depth", ylab = "Summer Seas-surface Temperature")
Run the code above in your browser using DataLab