RChronoModel (version 0.1)

PhaseStatistics: Summary statistics for a phase

Description

Estimation of several summary statistics of the beginning, the end and the duration of a phase.

Usage

PhaseStatistics(PhaseBeginning_chain, PhaseEnd_chain, level = 0.95,
  max_decimal = 0)

Arguments

PhaseBeginning_chain

numeric vector containing the output of the MCMC algorithm for the beginning of the phase

PhaseEnd_chain

numeric vector containing the output of the MCMC algorithm for the end of the phase

level

probability corresponding to the level of confidence used for the credible interval and the highest density region

max_decimal

maximum number of decimal

Value

• Returns a list of values corresponding to the summary statistics:
• 1
  Statistics of the beginning of the phase
• 2
  Statistics of the end of the phase
• 3
  Statistics of the duration of the phase
• 4
  The endpoints of the phase time range

Details

The summary statistics are those given by MarginalStatistics function. The time range is given by PhaseTimeRange function. The duration of the phase is computed as follow duration = end - beginning at each iteration of the MCMC output.

Examples

data(Phases); attach(Phases)
  PhaseStatistics(Phase.1.alpha, Phase.1.beta, 0.95, 0)
  PhaseStatistics(Phase.2.alpha, Phase.2.beta, 0.95, 0)