spDates (version 1.1)

modelDates: Perform regression of archaeological dates on great circle distances from a hypothetical origin. Dates can be filtered to retain only the earliest dates per distance bins (Hamilton and Buchanan 2007). Bootstrap is executed to account for uncertainty in calibrated dates. Regression can be either reduced major axis or ordinary least squares. If using ordinary least squares, regression is performed both on time-versus-distance and on distance-versus-time.

Description

Perform regression of archaeological dates on great circle distances from a hypothetical origin. Dates can be filtered to retain only the earliest dates per distance bins (Hamilton and Buchanan 2007). Bootstrap is executed to account for uncertainty in calibrated dates. Regression can be either reduced major axis or ordinary least squares. If using ordinary least squares, regression is performed both on time-versus-distance and on distance-versus-time.

Usage

modelDates(
  ftrSites,
  c14bp,
  origin,
  binWidth = 0,
  nsim = 999,
  method = "rma",
  ncores = 1
)

Value

a dateModel object.

Arguments

ftrSites

A SpatialPointsDataFrame object with associated earliest C14 dates per site and respective calibrated distributions (CalDates objects) in a field named "cal". Result of applying filterDates() and calAll().

c14bp

A string. Name of the field with the radiocarbon ages in C14 BP format.

origin

A SpatialPointsDataFrame object. The site considered as hypothetical origin of expansion.

binWidth

A number. Width of the spatial bins in km, calculated as distance intervals from the hypothetical origin. Default is 0 (no bins).

nsim

A number. Number of simulations to be run during the bootstrapping procedure. Default is 999.

method

A string. Method to be used in the regression. One of "rma" or "ols". Default is "rma".

ncores

A number. Number of cores used for parallel processing. Default is 1.

Examples

# \donttest{
data(neof)
data(centers)
jericho <- centers[centers$Site=="Jericho",]
model <- modelDates(neof, "C14Age", jericho, method="ols")
# }
# \dontshow{
data(neof)
data(centers)
jericho <- centers[centers$Site=="Jericho",]
model <- modelDates(neof[1:2,], "C14Age", jericho, method="ols",
nsim=1, ncores=2)
# }