ArchaeoPhases (version 1.8)

phases_gap: 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.

Examples

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

# }