spd: Summed probability distributions (SPD) of radiocarbon dates.

View source: R/aggregation.R

spdR Documentation

Summed probability distributions (SPD) of radiocarbon dates.

Description

The function generates Summed probability distributions (SPD) of radiocarbon dates, with optional binning routine for controlling inter-site or inter-phase variation in sampling intensity.

Usage

spd(
  x,
  timeRange,
  bins = NA,
  datenormalised = FALSE,
  spdnormalised = FALSE,
  runm = NA,
  verbose = TRUE,
  edgeSize = 500
)

Arguments

x

A CalDates class object containing the calibrated radiocarbon dates.

timeRange

A vector of length 2 indicating the start and end date of the analysis in cal BP.

bins

A vector containing the bin names associated with each radiocarbon date. If set to NA, binning is not carried out.

datenormalised

Controls for calibrated dates with probability mass outside the timerange of analysis. If set to TRUE the total probability mass within the time-span of analysis is normalised to sum to unity. Should be set to FALSE when the parameter normalised in calibrate is set to FALSE. Default is FALSE.

spdnormalised

A logical variable indicating whether the total probability mass of the SPD is normalised to sum to unity.

runm

A number indicating the window size of the moving average to smooth the SPD. If set to NA no moving average is applied. Default is NA

verbose

A logical variable indicating whether extra information on progress should be reported. Default is TRUE.

edgeSize

Extra margin in C14 Age time to handle edge effect when datenormalise is set to TRUE. Default is 500.

Details

The binning routine consists of computing summed probability distribution of all dates associated to a given bin, divided by the number of contributing dates. This controls for any striking differences in sampling intensity, and ensures that each site phase is equally contributing to the final SPD (see Timpson et al 2014 for details). Bins can be generated using the binPrep, whilst the sensitivity to parameter choice can be explored with the binsense function.

Value

An object of class CalSPD with the following elements

  • metadata A data.frame containing relevant information regarding the parameters used to create the SPD as well as sample size and number of bins

  • grid A CalGrid class object containing the summed probability associated to each calendar year between timeRange[1] and timeRange[2]

References

Timpson, A., et al, (2014). Reconstructing regional population fluctuations in the European Neolithic using radiocarbon dates: a new case-study using an improved method. Journal of Archaeological Science 52: 549-557. DOI:10.1016/j.jas.2014.08.011

See Also

calibrate for calibrating radiocarbon dates; binPrep for preparing bins.


rcarbon documentation built on Aug. 24, 2023, 5:11 p.m.