older: Bayesian Test for Anteriority/Posteriority

In ArchaeoStat/ArchaeoPhases: Post-Processing of Markov Chain Monte Carlo Simulations for Chronological Modelling

Bayesian Test for Anteriority/Posteriority

Description

A Bayesian test for checking the following assumption: "event x is older than event y".

Usage

older(x, y, ...)

## S4 method for signature 'numeric,numeric'
older(x, y)

## S4 method for signature 'EventsMCMC,missing'
older(x, y)

Arguments

x	A numeric vector giving the output of the MCMC algorithm for the first parameter, or an EventsMCMC object.
y	A numeric vector giving the output of the MCMC algorithm for the second parameter.
...	Currently not used.

Details

For a given output of MCMC algorithm, this function estimates the posterior probability of the event x < y by the relative frequency of the event "the value of event x is less than the value of event y" in the simulated Markov chain.

Methods (by class)

• older(x = numeric, y = numeric): Returns a length-one numeric vector (the posterior probability of the assumption: "event x is older than event y").

• older(x = EventsMCMC, y = missing): Returns a numeric matrix of posterior probabilities.

Author(s)

A. Philippe, M.-A. Vibet, N. Frerebeau

Examples

## Coerce to MCMC
eve <- as_events(mcmc_events, calendar = CE(), iteration = 1)
eve <- eve[1:10000, ]

## Test for anteriority
older(eve)

## Test for hiatus
hia <- hiatus(eve)
as.data.frame(hia)

ArchaeoStat/ArchaeoPhases documentation built on Jan. 19, 2025, 2:33 p.m.