dates_hiatus {ArchaeoPhases} | R Documentation |
Determines whether there is a hiatus between two MCMC chains and returns the longest interval that satisfies: P(a_chain < IntervalInf < IntervalSup < b_chain | M) = level
dates_hiatus(a_chain, b_chain, level = 0.95)
a_chain |
: Numeric vector containing the output of the MCMC algorithm for the first parameter. |
b_chain |
: Numeric vector containing the output of the MCMC algorithm for the second parameter. |
level |
Probability corresponding to the confidence level of the interval. |
A list with the following components:
A named vector where inf
is the lower endpoint of the
hiatus as a calendar year (AD/BC) or NA
if there is no hiatus at
level
, and sup
is the upper endpoint of the gap as a calendar
year (AD/BC), or NA
if there is no hiatus at level
.
The duration of the hiatus at level
.
Probability corresponding to the confidence level of the interval.
The function call.
Anne Philippe, Anne.Philippe@univ-nantes.fr,
Marie-Anne Vibet, Marie-Anne.Vibet@univ-nantes.fr, and
Thomas S. Dye, tsd@tsdye.online
data(Events); attach(Events) dates_hiatus(Event.1, Event.12) dates_hiatus(Event.1, Event.12, level = 0.5)