phases_gap: Gap or hiatus between two successive phases (for phases in...

In ArchaeoPhases: Post-Processing of the Markov Chain Simulated by 'ChronoModel', 'Oxcal' or 'BCal'

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, [e-mail],

Marie-Anne Vibet, [e-mail], and

Thomas S. Dye, [e-mail]

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)

ArchaeoPhases documentation built on June 22, 2022, 1:05 a.m.