ArchaeoPhases (version 1.8)

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

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).

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.

Examples

# NOT RUN {
  # 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))

# }