Reads MCMC output.
Usage
read_chronomodel_events(file, ...)
read_chronomodel_phases(file, ...)
read_chronomodel_model(file, ...)
# S4 method for character
read_chronomodel_events(file, BP = FALSE, sep = ",", dec = ".")
# S4 method for character
read_chronomodel_phases(file, BP = FALSE, sep = ",", dec = ".")
# S4 method for character
read_chronomodel_model(file)
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 (seefile
).Alternatively,
file
can be a readable text-mode connection (which will be opened for reading if necessary, and if soclose
d (and hence destroyed) at the end of the function call). (Ifstdin()
is used, the prompts for lines may be somewhat confusing. Terminate input with a blank line or an EOF signal,Ctrl-D
on Unix andCtrl-Z
on Windows. Any pushback onstdin()
will be cleared before return.)file
can also be a complete URL. (For the supported URL schemes, see the ‘URLs’ section of the help forurl
.)- ...
Further arguments to be passed to
read.table
.- BP
A
logical
scalar: should the data be converted from BP to CE? This should not beTRUE
unless you change the default settings in 'OxCal' or 'ChronoModel'.- sep
the field separator character. Values on each line of the file are separated by this character. If
sep = ""
(the default forread.table
) the separator is ‘white space’, that is one or more spaces, tabs, newlines or carriage returns.- dec
the character used in the file for decimal points.
Value
An EventsMCMC
or a PhasesMCMC
object.
References
Lanos, Ph., Philippe, A. & Dufresne, Ph. (2015). Chronomodel: Chronological Modeling of Archaeological Data using Bayesian Statistics. URL: https://www.chronomodel.fr.
See also
Other read methods:
check
,
coda
,
coerce()
,
read_bcal()
,
read_oxcal()
Examples
if (FALSE) {
## Import ChronoModel Output
path <- "chronomodel/ksarakil"
## Events
path_events <- system.file(path, "Chain_all_Events.csv", package = "fasti")
(chrono_events <- read_chronomodel_events(path_events))
## Phases
path_phases <- system.file(path, "Chain_all_Phases.csv", package = "fasti")
(chrono_phases <- read_chronomodel_phases(path_phases))
}