MultiPhasesGap {RChronoModel}

Gap/Hiatus between a succession of phases (for phases in temporal order constraint)

Description

This function finds, if it exists, the gap between two successive phases. This gap or hiatus is the longest interval [IntervalInf, IntervalSup] that satisfies : P(Phase1Max < IntervalInf < IntervalSup < Phase2Min | M) = level for each successive phase.

Usage

MultiPhasesGap(data, position_minimum, position_maximum = position_minimum+1,
 level = 0.95, max_decimal = 0)

Arguments

data

dataframe containing the output of the MCMC algorithm

position_minimum

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

position_maximum

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

level

probability corresponding to the level of confidence

max_decimal

maximum number of decimal

Details

For each i, MultiPhasesGap computes the gap 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

Returns a matrix of values containing the level of confidence and the endpoints of the gap for each pair of successive phases

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 phases 
  names(Phases)
  
  # Stipulating position_maximum
  MultiPhasesGap(Phases, c(4,2), c(5,3))
  
  # In this case, equivalent to 
  MultiPhasesGap(Phases, c(4,2))