plot_pbsEDM_Evec: Plot output from 'pbsEDM_Evec()' as six panel plot, with data...

View source: R/plotting.R

plot_pbsEDM_EvecR Documentation

Plot output from 'pbsEDM_Evec()' as six panel plot, with data shown the same as for 'plot.pbsEDM()' and sixth panel with predictions vs observation for each 'E'

Description

Plots time series of 'N_t' and 'Y_t', 2d phase plots of 'N_t' vs 'N_t-1' and 'Y_t' vs 'Y_t-1', and 3d phase plot of 'Y_t' vs 'Y_t-1' vs 'Y_t-2'. Sixth panel shows predictions vs observations for different values of 'E'. See vignette "Analyse a simple time series" for full details.

Usage

plot_pbsEDM_Evec(E_res, ...)

Arguments

E_res

List of 'pbsEDM' objects as output from 'pbsEDM_Evec()'

...

extra arguments to 'plot.pbsEDM()', such as label.cex

Value

Six-panel figure in current plot environment

Author(s)

Andrew Edwards

Examples


  aa_Evec <- pbsEDM_Evec(NY_lags_example$N_t)
  plot_pbsEDM_Evec(aa_Evec, portrait = FALSE)

# For manuscript figures:
E_results <- pbsEDM_Evec(NY_lags_example$N_t)

postscript("six_panels_all.eps",
           height = 5.36,
           width = 9,
           horizontal=FALSE,
           paper="special")
plot_pbsEDM_Evec(E_results,
    last.time.to.plot = NULL, # to automatically plot all
    portrait = FALSE)
dev.off()


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