MultiPhaseTimeRange {ArchaeoPhases}R Documentation

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, Anne.Philippe@univ-nantes.fr and

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

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


[Package ArchaeoPhases version 1.5 Index]