plot_phase_3d: Plot 3d phase plot of 'Y_t' vs 'Y_t-1' vs 'Y_t-2'

In luke-a-rogers/pbsedm: An R package to implement some of the methods of empirical dynamic modelling

Plot 3d phase plot of 'Y_t' vs 'Y_t-1' vs 'Y_t-2'

Description

Shows the first-differenced values in three-dimensional phase space, as points showing lags of 0, 1 and 2. Cobwebbing shows how one point iterates to the next point in the time series. Can also highlight a $t^*$ value and the points that should be excluded from the library of nearest neighbours.

Usage

plot_phase_3d(
  obj,
  par.mgp.3d = c(3, 10, 0),
  par.mai.3d = c(0.1, 0.1, 0.1, 0.1),
  par.mar.3d = c(3, 0, 0, 0),
  x.lab = expression(italic(Y)[t - 2]),
  y.lab = expression(italic(Y)[t - 1]),
  z.lab = expression(italic(Y)[t]),
  last.time.to.plot = NULL,
  axis.range = NA,
  late.num = 3,
  pt.type = "p",
  late.col = "red",
  early.col = "black",
  early.col.lines = "lightgrey",
  axes.col = "darkblue",
  par.mar.phase = c(3, 0, 1, 0),
  par.mgp = c(1.5, 0.5, 0),
  tstar = NA,
  angle_view = NA,
  label = NULL,
  label.cex = 0.7
)

Arguments

obj	'pbsEDM' object
par.mgp.3d	'par(mgp)' values
par.mai.3d	'par(mai)' values
par.mar.3d	'par(mar)' values
x.lab	x axis label
y.lab	y axis label
z.lab	z axis label
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])
axis.range	range of axes, if NA then calculated from values; all three axes have the same range
late.num	final number of 'N_t' time steps to plot in a different colour; not showing 'N_t' here but keeping consistency with other plots
pt.type	'type' value for 'points()'
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
axes.col	colour of orthogonal origin axes that go through (0, 0, 0)
par.mar.phase	'par(mar)' to reset to after plotting 3d figure
par.mgp	'par(mgp)' to reset to after plotting 3d figure
tstar	focal point to highlight (for 'pbsSmap' vignette)
angle_view	manually specify the angle for viewing (to rotate plot)
label	label to annotate plot, such as '(a)' etc. for six-panel figure
label.cex	size of label annotation

Value

Plots figure to current device

Author(s)

Andrew Edwards

Examples

 aa <- pbsEDM(NY_lags_example,
              lags = list(N_t = 0:2),
              first_difference = TRUE)
 plot_phase_3d(aa)

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