phases_gap {ArchaeoPhases} | R Documentation |
This function finds, if it exists, a gap or hiatus between two successive phases.
This gap or hiatus is the longest interval that satisfies
P(Phase1Max_chain < IntervalInf < IntervalSup < Phase2Min_chain | M) = level
phases_gap(a_chain, b_chain, level = 0.95)
a_chain |
Numeric vector containing the output of the MCMC algorithm for the upper boundary of the older phase. |
b_chain |
Numeric vector containing the output of the MCMC algorithm for the lower boundary of the younger phase. |
level |
Probability corresponding to the level of confidence. |
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
.
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(Phases); attach(Phases)
phases_gap(Phase.1.beta, Phase.2.alpha, 0.95)
phases_gap(Phase.1.beta, Phase.2.alpha, 0.50)