plot_time_series: Plot the observed time series as either 'N_t' or 'Y_t'

plot_time_seriesR Documentation

Plot the observed time series as either 'N_t' or 'Y_t'

Description

First value must be 't=1'. For non-differenced values 'N_t', shows the time series with the final values in a different colour, and a title showing the final time step. For first-differenced values 'Y_t', shows the time series of those, plus a one-dimensional phase plot.

Usage

plot_time_series(
  values,
  X.or.N,
  par.mar.ts = c(3, 3, 1, 1),
  max_time = NULL,
  t.axis.range = NULL,
  last.time.to.plot = NULL,
  late.num = 3,
  late.col = "red",
  early.col = "black",
  early.col.lines = "lightgrey",
  start = 1,
  pt.type = "p",
  par.mgp = c(1.5, 0.5, 0),
  add.legend = TRUE,
  label = NULL,
  label.cex = 0.7
)

Arguments

values

vector of values to be plotted

X.or.N

"N" if raw non-differenced data, "X" for differenced data

par.mar.ts

'par(mar)' values

max_time

maximum time value for the time axis

t.axis.range

range of time axis

last.time.to.plot

last time value of N[t] to use when plotting, so final Y[t] used will be Y[t-1] (since Y[t] uses N[t+1])

late.num

final number of 'N[t]' time steps to plot in a different colour

late.col

colour in which to plot final 'late.num' time steps

early.col

colour in which to plot earlier time step points

early.col.lines

colour in which to plot earlier time step points

start

first time step (must be 1)

pt.type

'type' value for 'points()'

par.mgp

'par("mgp")' values

add.legend

logical, whether to add legend to the 'N[t]' time series

label

label to annotate plot, such as '(a)' etc. for six-panel figure

label.cex

size of label annotation

Author(s)

Andrew Edwards

Examples


  plot_time_series(NY_lags_example$N_t, X.or.N = "N")
  plot_time_series(NY_lags_example$Y_t, X.or.N = "X")


luke-a-rogers/pbsedm documentation built on June 3, 2024, 5:20 a.m.