A thin wrapper of rcarbon::calibrate() that returns calibrated dates as a single list rather than a `CalDates`` object. Useful, for example, if you want to add a column of calibrated dates to an existing table with dplyr::mutate()

c14_calibrate(cra, error, ..., engine = c("rcarbon", "OxCal", "Bchron"))

Arguments

cra

A vector of uncalibrated radiocarbon ages.

error

A vector of standard errors associated with cra

...

Optional arguments passed to calibration function.

engine

Package to use for calibration, i.e. rcarbon::calibrate() ("rcarbon"), oxcAAR::oxcalCalibrate() ("OxCal"), or Bchron::BchronCalibrate() ("Bchron"). Default: "rcarbon".

Value

A list of cal objects.

See also

Other tidy radiocarbon functions: c14_sum(), cal()

Examples

data("shub1_radiocarbon") shub1_radiocarbon %>% dplyr::mutate(cal = c14_calibrate(cra, error))
#> [1] "Calibrating radiocarbon ages..." #> | | | 0% | |=== | 4% | |===== | 7% | |======== | 11% | |========== | 15% | |============= | 19% | |================ | 22% | |================== | 26% | |===================== | 30% | |======================= | 33% | |========================== | 37% | |============================= | 41% | |=============================== | 44% | |================================== | 48% | |==================================== | 52% | |======================================= | 56% | |========================================= | 59% | |============================================ | 63% | |=============================================== | 67% | |================================================= | 70% | |==================================================== | 74% | |====================================================== | 78% | |========================================================= | 81% | |============================================================ | 85% | |============================================================== | 89% | |================================================================= | 93% | |=================================================================== | 96% | |======================================================================| 100% #> [1] "Done."
#> Registered S3 method overwritten by 'cli': #> method from #> print.boxx spatstat
#> # A tibble: 27 x 9 #> lab_id context phase sample material cra error outlier cal #> <chr> <int> <chr> <chr> <chr> <int> <int> <lgl> <list> #> 1 RTD-7951 23 Phase… Context 126… Bolboschoe… 12166 55 FALSE <cal[,2… #> 2 Beta-11… 24 Phase… SHUB1/105 gazelle du… 12310 60 FALSE <cal[,2… #> 3 RTD-7317 26 Phase… Context 83,… Bolboschoe… 12289 46 FALSE <cal[,2… #> 4 RTD-7318 27 Phase… Context 86,… Zilla sp. … 12332 46 FALSE <cal[,2… #> 5 RTD-7948 24 Phase… Context 120… Bolboschoe… 12478 38 FALSE <cal[,2… #> 6 RTD-7947 22 Phase… Context 166… Bolboschoe… 12322 38 FALSE <cal[,2… #> 7 RTD-7313 22 Phase… Context 67,… Bolboschoe… 12346 46 FALSE <cal[,2… #> 8 RTD-7311 22 Phase… Context 67,… Vitex sp. … 12367 65 FALSE <cal[,2… #> 9 RTD-7312 22 Phase… Context 67,… Zilla sp. … 12405 50 FALSE <cal[,2… #> 10 RTD-7314 22 Phase… Context 67,… Bolboschoe… 12273 48 FALSE <cal[,2… #> # … with 17 more rows