Plots the first derivative of the tempo
plot Bayesian estimate.
activity(object, ...)
# S4 method for class 'EventsMCMC'
activity(
object,
from = min(object),
to = max(object),
grid = getOption("ArchaeoPhases.grid")
)
# S4 method for class 'CumulativeEvents'
activity(object)
# S4 method for class '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,
...
)
An EventsMCMC
or a CumulativeEvents
object.
Other graphical parameters may also be passed as
arguments to this function, particularly, border
, col
, lwd
or lty
.
A length-one numeric
vector giving the earliest date to
estimate for (in years).
A length-one numeric
vector giving the latest date to estimate
for (in years).
A length-one numeric
vector specifying the number of equally
spaced points of the temporal grid.
An ActivityEvents
object.
A aion::TimeScale
object specifying the target
calendar (see calendar()
).
A character
string giving a main title for the plot.
A character
string giving a subtitle for the plot.
A logical
scalar: should the default annotation (title and x
and y axis labels) appear on the plot?
A logical
scalar: should axes be drawn on the plot?
A logical
scalar: should a box be drawn around the
plot?
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.
An expression
to be evaluated after plotting has taken
place but before the axes, title and box are added.
activity()
returns an ActivityEvents
object.
plot()
is called it for its side-effects: it results in a graphic being
displayed (invisibly returns x
).
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 .
Other event tools:
elapse()
,
occurrence()
,
tempo()
## 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())