transformSPD function

Apply taphonomic corrections or other transformations to an SPD.

transformSPD(
  x,
  correction = expression(PrDens/(5.726442 * 10^6 * (CalBP + 2176.4)^-1.3925309))
)

Arguments

• x: An object of class CalSPD, compositeKDE or stackCalSPD.
• correction: An expression for transforming the SPD. Available input terms include: CalBP, the vector of calBP year within the time range; and PrDens, a matching vector of summed probability. The default expression is the taphonomic correction formula proposed by Surovell et al 2009.

Returns

An object of the same class as x

Examples

## Not run:

data(emedyd)
region1 = subset(emedyd,Region==1)
x = calibrate(x=region1$CRA, errors=region1$Error,normalised=FALSE)
bins = binPrep(sites=region1$SiteName, ages=region1$CRA,h=50)
region1.spd = spd(x=x,bins=bins,timeRange=c(16000,8000))
region1.spd.corrected = transformSPD(region1.spd)
## End(Not run)

References

Surovell, T.A., Finley, J.B., Smith, G.M., Brantingham, P.J., Kelly, R., 2009. Correcting temporal frequency distributions for taphonomic bias. Journal of Archaeological Science 36, 1715–1724.