Import MCMC Samples
## Install the latest version install.packages("fasti", repos = "https://tesselle.r-universe.dev")
Two different files are generated by ChronoModel:
Chain_all_Events.csv that contains the MCMC samples of each event created in the modeling, and
Chain_all_Phases.csv that contains all the MCMC samples of the minimum and the maximum of each group of dates if at least one group is created.
## Read events from ChronoModel output_events <- system.file("chronomodel/ksarakil/Chain_all_Events.csv", package = "fasti") chrono_events <- read_chronomodel_events(output_events) ## Plot events plot(chrono_events) #> Picking joint bandwidth of 49.2
## Read phases from ChronoModel output_phases <- system.file("chronomodel/ksarakil/Chain_all_Phases.csv", package = "fasti") chrono_phases <- read_chronomodel_phases(output_phases) ## Plot phases plot(chrono_phases)
Oxcal generates a CSV file containing the MCMC samples of all parameters (dates, start and end of phases).
## Read OxCal MCMC samples output_oxcal <- system.file("oxcal/ksarakil/MCMC_Sample.csv", package = "fasti") oxcal_mcmc <- read_oxcal(output_oxcal)
The phase boundaries cannot be extracted automatically from Oxcal output. Use
as_phases() to get the phase boundaries:
BCal generates a CSV file containing the MCMC samples of all parameters (dates, start and end of groups).
## Read BCal MCMC samples output_bcal <- system.file("bcal/output/rawmcmc.csv", package = "fasti") bcal_mcmc <- read_bcal(output_bcal)
The group boundaries cannot be extracted automatically from BCal output. Use
as_phases() to get the group boundaries: