MultiPhaseTimeRange: Phase time range for multiple groups

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

Phase time range for multiple groups

Description

Computes the shortest interval that satisfies P(PhaseMin < IntervalInf < IntervalSup < PhaseMax | M) = level for each phase

Usage

MultiPhaseTimeRange(
  data,
  position_minimum,
  position_maximum = position_minimum + 1,
  level = 0.95
)

Arguments

data	Data frame containing the output of the MCMC algorithm.
position_minimum	Numeric vector containing the column number corresponding to the minimum of the events included in each phase.
position_maximum	Numeric vector containing the column number corresponding to the maximum of the phases set in the same order as in position_minimum.
level	Probability corresponding to the desired level of confidence.

Details

For each i, MultiPhaseTimeRange() computes the time range interval for the phase defined by its minimum position_minimum[i] and its maximum position_maximum[i]. The default value of position_maximum corresponds to CSV files exported from ChronoModel software.

Value

A matrix of values containing the level of confidence and the endpoints of the shortest time range associated with the desired level. The result is given in calendar years (BC/AD).

Author(s)

Anne Philippe, [e-mail] and

Marie-Anne Vibet, [e-mail]

Examples

  # Data extracted from ChronoModel software
  data(Phases)
  # List of the name of the groups
  names(Phases)
  # Stipulating position_maximum
  MultiPhaseTimeRange(Phases, position_minimum = c(4, 2), position_maximum = c(5, 3))
  # In this case, equivalent to
  MultiPhaseTimeRange(Phases, position_minimum = c(4, 2))

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