Density of Event and Accumulation Dates

Estimates the event and accumulation density.

Usage

density_event(object, ...)

density_accumulation(object, ...)

# S4 method for class 'EventDate'
density_event(object, dates = NULL, calendar = NULL, n = 500, ...)

# S4 method for class 'EventDate'
density_accumulation(
  object,
  dates = NULL,
  calendar = NULL,
  type = c("activity", "tempo"),
  n = 500,
  ...
)

Arguments

object

An EventDate object.

...

Currently not used.

dates

A numeric vector of dates expressed as calendar years or rata die (if calendar is NULL).

calendar

An aion::TimeScale object specifying the calendar of dates (see aion::calendar()). If NULL (the default), rata die are expected.

n

A length-one non-negative numeric vector giving the desired length of the vector of quantiles for density computation.

type

A character string indicating the type of plot. It must be one of "activity" (default) or "tempo" (see details). Any unambiguous substring can be given.

Value

An aion::TimeSeries object.

See also

Other event date tools: event(), model_event, predict_event()

Author

N. Frerebeau

Examples

## Data from Peeples and Schachner 2012
data("zuni", package = "folio")

## Assume that some assemblages are reliably dated (this is NOT a real example)
zuni_dates <- c(
  LZ0569 = 1097, LZ0279 = 1119, CS16 = 1328, LZ0066 = 1111,
  LZ0852 = 1216, LZ1209 = 1251, CS144 = 1262, LZ0563 = 1206,
  LZ0329 = 1076, LZ0005Q = 859, LZ0322 = 1109, LZ0067 = 863,
  LZ0578 = 1180, LZ0227 = 1104, LZ0610 = 1074
)

## Model the event and accumulation date for each assemblage
model <- event(zuni, zuni_dates, rank = 10)
plot(model, select = 1:10, event = TRUE, flip = TRUE)