Plots the first derivative of the tempo plot Bayesian estimate.

activity(object, ...)

# S4 method for EventsMCMC
activity(
  object,
  from = min(object),
  to = max(object),
  grid = getOption("ArchaeoPhases.grid")
)

# S4 method for CumulativeEvents
activity(object)

# S4 method for ActivityEvents,missing
plot(
  x,
  calendar = getOption("ArchaeoPhases.calendar"),
  main = NULL,
  sub = NULL,
  ann = graphics::par("ann"),
  axes = TRUE,
  frame.plot = axes,
  panel.first = NULL,
  panel.last = NULL,
  ...
)

Arguments

object

An EventsMCMC or a CumulativeEvents object.

...

Other graphical parameters may also be passed as arguments to this function, particularly, border, col, lwd or lty.

from

A length-one numeric vector giving the earliest date to estimate for (in years).

to

A length-one numeric vector giving the latest date to estimate for (in years).

grid

A length-one numeric vector specifying the number of equally spaced points of the temporal grid.

x

An ActivityEvents object.

calendar

A TimeScale object specifying the target calendar (see calendar()).

main

A character string giving a main title for the plot.

sub

A character string giving a subtitle for the plot.

ann

A logical scalar: should the default annotation (title and x and y axis labels) appear on the plot?

axes

A logical scalar: should axes be drawn on the plot?

frame.plot

A logical scalar: should a box be drawn around the plot?

panel.first

An an expression to be evaluated after the plot axes are set up but before any plotting takes place. This can be useful for drawing background grids.

panel.last

An expression to be evaluated after plotting has taken place but before the axes, title and box are added.

Value

  • activity() returns an ActivityEvents object.

  • plot() is called it for its side-effects: it results in a graphic being displayed (invisibly returns x).

References

Dye, T. S. (2016). Long-term rhythms in the development of Hawaiian social stratification. Journal of Archaeological Science, 71: 1-9. doi:10.1016/j.jas.2016.05.006 .

See also

Other event tools: elapse(), occurrence(), tempo()

Author

A. Philippe, M.-A. Vibet, T. S. Dye, N. Frerebeau

Examples

## Coerce to MCMC
eve <- as_events(mcmc_events, calendar = CE(), iteration = 1)
eve <- eve[1:10000, ]

## Tempo plot
tmp <- tempo(eve)
plot(tmp)

plot(tmp, interval = "credible", panel.first = grid())

plot(tmp, interval = "gauss", panel.first = grid())


## Activity plot
act <- activity(tmp)
plot(act, panel.first = grid())