clam - an introduction

clam¹ is an age-depth modelling tool which calibrates any radiocarbon dates, and then repeatedly chooses random age estimates from the (radiocarbon or other) dates and draws curves through them in order to estimate ages and their uncertainties.

clam works through R (or Rstudio). The first time you are using clam on your computer, you will have to install its R package, by typing within the terminal of R:

install.packages('clam')

The above command will also download the companion packages rice and rintcal, which deal with radiocarbon calibration aspects. To use clam, first the package has to be loaded:

library(clam)

## Loading required package: rice

## Loading required package: rintcal

The package comes with a pre-loaded dataset, called Example. Let’s start by running this default core, using all default settings. Note R’s comments about where the files are going to be placed, as we will need this information later. Press y and/or Enter to accept any suggestions about where it will place files, or alternatively provide a folder location using the option coredir.

clam()

## The run's files will be put in this folder: /private/var/folders/w_/ktch6mh953z13lmhs4_9bj_80000gn/T/RtmpFCaMlf/clam_runs/Example
## The run's files will be put in this folder: /private/var/folders/w_/ktch6mh953z13lmhs4_9bj_80000gn/T/RtmpFCaMlf/clam_runs/Example

## 
##  Warning, some dates lie partly outside the calibration curve!

##  Calibrating dates...

##  Interpolating, sampling.....

## 
## Example's 95% confidence ranges span from 20 to 531 yr (average 234 yr)

##   Fit (-log, lower is better): 6.78

clam will calibrate any C-14 dates, sample the age estimates, draw curves, calculate age ranges, produce graphs and provide information about the confidence ranges. The output graph, above, contains: the calibrated distributions in blue (and a green/blue cal BP for the core’s surface), the age-depth model’s 95% ranges (grey envelope) and its mean (black curve). Besides the graph, clam also produces files containing age estimates (95% ranges and the mean) for each depth (by default every cm from the top to the bottom core depth). This file can be found in the core’s folder, together with files giving the calibrated ranges for the dated depths.

Please don’t over-interpret the Fit value, which is also reported. Sometimes age-depth models with low fit values, so with a good match of the dates to the model, will still be unrealistic, e.g. with very abrupt changes or even age reversals. Note also that the uncertainty estimates provided by classical modelling are likely to be underestimates of the true uncertainties², and that Bayesian age-models such as OxCal’s P Sequence or R packages rbacon or BChron are more robust alternatives which provide more realistic error estimates.

There are many options to produce different types of age-depth models. The default, type=1, is linear interpolation between the dated depths. Alternatives are linear or higher polynomial regression (type=2, the degree is set with the smooth parameter), cubic spline (type=3), smooth spline (type=4) or locally weighted spline (type=5). Hiatuses can be set at specific depths using hiatus, and bits where sedimentation was abrupt can also be modelled (slump, e.g., tephra layers). If some dates appear outlying they can be labelled such and excluded from the age-model (e.g. if the sixth and ninth dates counting from the top one appear outlying, type outliers=c(6,9)). If you prefer your ages as BC/AD, use BCAD=TRUE. Note that such decisions can have large impacts on the resulting age-depth model - please be critical about your assumptions.

By default, ages are calculated for every cm from the top to the bottom dated depth. This can be adapted by specifying dmin, dmax, and every (default 1). A file with depths can also be provided. This must be in the core’s folder (see later), start with the core’s name, and end in _depths.txt. Then use depths.file=TRUE as option.

Here is an example using some of the above options:

clam(ask=FALSE, type=4, outliers=6, hiatus=470, dmax=800, slump=c(220, 250), coredir=clam_dir)

## The run's files will be put in this folder: /private/var/folders/w_/ktch6mh953z13lmhs4_9bj_80000gn/T/RtmpFCaMlf/clam_runs/Example

## 
##  Warning, some dates lie partly outside the calibration curve!

##  Calibrating dates...

## 
##  section 1,

##  Using smoothing spline (smoothing 0.3), sampling

## .....

## 
##  section 2,

##  extrapolating beyond dated levels, dangerous!
## 

##  Using smoothing spline (smoothing 0.3), sampling

## .....

## 
##  Warning, some dates lie partly outside the calibration curve!

## 
## Example's 95% confidence ranges span from 25 to 1140 yr (average 364 yr)

##   Fit (-log, lower is better): 7.45

When clam runs it will tell you where the files will be stored. The default is to use an ‘umbrella’ folder called clam_runs or Cores, and store the runs as folders within it. After running the default core, a folder ‘Example’ will appear in the umbrella folder, and this will contain the input and output files of the clam run.

To run your own cores, make a new folder within the umbrella folder, and save it under a name, e.g., MyCore. Then place a file called MyCore.csv in that folder. This file will contain the dates and should have the same formatting as the Example.csv file in the Example folder. It has headers and six fields (or seven if you want to include a variable thickness which will draw a rectangle showing the thickness of the dated depth):

Note that depths with radiocarbon dates should have the dates in the second/radiocarbon column, and no entry in the cal age column. Dates already on the calendar scale should be empty in the second column and have their entry in the third/cal age column. Leave the offset column empty if no offset is assumed. All dates must have an error estimate - this cannot be left empty. If you produce your file in a spreadsheet program such as MS-Excel or LibreOffice, please also check its formatting in a plain-text editor to see if all looks OK (e.g., no orphan quotation marks, no empty lines, not too many commas).

Then run clam("MyCore").

For more information and options, run ?clam. If you have any questions or ideas for improvements, please contact me at maarten.blaauw@qub.ac.uk.

If you use clam in your publications, please cite Blaauw 2010³, the version used (e.g., 2.6.1), any non-default settings applied, and for radiocarbon dates their calibration curves (e.g., IntCal20⁴).