phases_gap {ArchaeoPhases}

R Documentation

Gap or hiatus between two successive phases (for phases in temporal order constraint)

Description

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

Usage

phases_gap(a_chain, b_chain, level = 0.95)

Arguments

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.

Value

A list with the following components:

hiatus

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.

level

Probability corresponding to the confidence level of the interval.

call

The function call.

Author(s)

Anne Philippe, Anne.Philippe@univ-nantes.fr,

Marie-Anne Vibet, Marie-Anne.Vibet@univ-nantes.fr, and

Thomas S. Dye, tsd@tsdye.online

Examples

  data(Phases); attach(Phases)
  phases_gap(Phase.1.beta, Phase.2.alpha, 0.95)
  phases_gap(Phase.1.beta, Phase.2.alpha, 0.50)

---