stackspd: Stacked Summed Probability Distribution

View source: R/aggregation.R

stackspdR Documentation

Stacked Summed Probability Distribution

Description

Generates and combines multiple SPDs based on a user defined grouping.

Usage

stackspd(
  x,
  timeRange,
  bins = NA,
  group = NULL,
  datenormalised = 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.

group

A character or factor vector containing the grouping variable.

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.

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

Controls edge effect by expanding the fitted model beyond the range defined by timeRange.

Value

An object of class stackCalSPD

Examples

## Not run: 
data(emedyd)
x = calibrate(x=emedyd$CRA, errors=emedyd$Error,normalised=FALSE)
bins = binPrep(sites=emedyd$SiteName, ages=emedyd$CRA,h=50)
res = stackspd(x=x,timeRange=c(16000,8000),bins=bins,group=emedyd$Region)

## End(Not run)

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