## Continue example from ?join
example(join)
## Fit a MAT model
swap.mat <- mat(swapdiat, swappH, method = "SQchord")
swap.mat
## Reconstruct pH for the RLGH core
rlgh.pH <- predict(swap.mat, rlgh, bootstrap = TRUE,
n.boot = 100)
## draw the reconstruction
reconPlot(rlgh.pH, use.labels = TRUE, display.error = "bars",
xlab = "Depth", ylab = "pH")
Run the code above in your browser using DataLab