DatesHiatus: Test for the existence of a hiatus between two parameters

In RChronoModel: Post-Processing of the Markov Chain Simulated by ChronoModel or Oxcal

Description

Finds if it exists a gap between two dates that is the longest interval that satisfies : P(a_chain < IntervalInf < IntervalSup < b_chain | M) = level

Usage

DatesHiatus(a_chain, b_chain, level=0.95)

Arguments

a_chain	numeric vector containing the output of the MCMC algorithm for the first one-parameter (date) a
b_chain	numeric vector containing the output of the same MCMC algorithm for the second one-parameter (date) b
level	probability corresponding to the level of confidence used for the credible interval and the highest density region

Value

Returns the endpoints of the longest hiatus between two parameters

Author(s)

Anne Philippe <[e-mail]> and

Marie-Anne Vibet <[e-mail]>

Examples

  data(Events); attach(Events)
  DatesHiatus(Event.1, Event.12)
  DatesHiatus(Event.1, Event.12, level = 0.5)

Example output

            level HiatusIntervalInf HiatusIntervalSup 
           "0.95"              "NA"              "NA" 
            level HiatusIntervalInf HiatusIntervalSup 
              0.5           -1244.0           -1235.0 

RChronoModel documentation built on May 1, 2019, 10:06 p.m.