ArchaeoPhases (version 1.4)

MultiPhaseTimeRange: Phase Time Range for multiple groups

Description

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

Usage

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

Arguments

data

dataframe containing the output of the MCMC algorithm. The MCMC samples should be in calendar year (BC/AD).

position_minimum

numeric vector containing the column number corresponding to the minimum of the dates included in each group

position_maximum

numeric vector containing the column number corresponding to the maximum of the dates included in each group. By default, position_maximum = position_minimum + 1.

level

probability corresponding to the desired level of confidence

Value

Returns 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 year (in format 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))
# }