Read OxCal Output

Reads MCMC output.

Usage

read_oxcal(file, ...)

# S4 method for character
read_oxcal(file, BP = FALSE)

Arguments

file

the name of the file which the data are to be read from. Each row of the table appears as one line of the file. If it does not contain an absolute path, the file name is relative to the current working directory, getwd(). Tilde-expansion is performed where supported. This can be a compressed file (see file).

Alternatively, file can be a readable text-mode connection (which will be opened for reading if necessary, and if so closed (and hence destroyed) at the end of the function call). (If stdin() is used, the prompts for lines may be somewhat confusing. Terminate input with a blank line or an EOF signal, Ctrl-D on Unix and Ctrl-Z on Windows. Any pushback on stdin() will be cleared before return.)

file can also be a complete URL. (For the supported URL schemes, see the ‘URLs’ section of the help for url.)

...

Further arguments to be passed to read.table.

BP

A logical scalar: should the data be converted from BP to CE? This should not be TRUE unless you change the default settings in 'OxCal' or 'ChronoModel'.

Value

An EventsMCMC object.

References

Bronk Ramsey, C. (2009). Bayesian Analysis of Radiocarbon Dates. Radiocarbon, 51(1), 337-360. doi:10.1017/S0033822200033865 .

Author

T. S. Dye, N. Frerebeau

Examples

if (FALSE) {
## Import OxCal Output
path_output <- system.file("oxcal/ksarakil/MCMC_Sample.csv", package = "fasti")
(oxcal <- read_oxcal(path_output))
}