# NOT RUN {
data(bassriver)
# Using the "pooled annual" calibration model:
sst <- predict_seatemp(bassriver$d18o, d18osw=0.0,
prior_mean=30.0, prior_std=20.0)
predictplot(x=bassriver$depth, y=sst, ylim=c(20, 40),
ylab="SST (<U+00B0>C)", xlab="Depth (m)")
# }
Run the code above in your browser using DataLab