For each sample, calculate the time it takes for the density to decrease by half from the peak (or by another use-provided ratio)

Source: R/bayesian_inference_functions.R

For each sample, calculate the time it takes for the density to decrease by half from the peak (or by another use-provided ratio)

calc_half_life_from_peak(
  bayesian_soln,
  density_model,
  rc_meas = list(),
  calib_df = list(),
  prop_change = 0.5,
  bayesian_summary = NA,
  peak_range = NA
)

Arguments

bayesian_soln

The solution, a list-like object of class bd_bayesian_soln (see sample_theta).
density_model

The density model (see set_density_model).
rc_meas

The radiocarbon measurements (see import_rc_data; optional: if not provided, bayesian_summary must be provided).
calib_df

The calibration data frame (see load_calib_curve; optional: if not provided, bayesian_summary must be provided).
prop_change

The relative decrease in density to use for the duration calculation (default: 0.5).
bayesian_summary

The result of a call to summarize_bayesian_inference. (optional; if not provided, it is calculated, which requires that rc_meas and calib_df be provided).
peak_range

A range over which to search for the peak of the density function (default: NA, which means that tau_min to tau_max is used for the range).

Value

A vector of "half-lives" (proportional change set by prop_change)

Details

For each sample, calculate the time it takes for the density to decrease by half from the peak. Optionally, a different proportion can be used than the default prop_change = 0.5. For example, with prop_change = 0.1 the time it takes for the density to decrease by 10% is used. If the relative density is not reached, the half life for the sample is set to NA. If there is no interior peak in the range peak_range, which is tau_min to tau_max by default, the half life is set to NA.

See also