plot.chaos01: Plot the additional results of 0-1 test for chaos.

Description

This function plot the Pc to Qc plot and Mc/Dc plot as described in Gottwald and Melbourne (2004).

Usage

# S3 method for chaos01
plot(x, plotvar = c("PQ", "MD", "BB", "COG"),
  mdcol = NULL, step = NULL, col2 = 2, col = 1, main = NULL,
  xlab = NULL, ylab = NULL, type = NULL, ylim = NULL, ...)

Arguments

the object of "chaos01" class, produced by testChaos01 function when parameter out = "TRUE". Subset the output of the function to get the results for the concrete c. See the example.

plotvar

list/vector define what should be plotted.

c("PQ", "MD", "BB", "COG") - plot Pc - Qc figure, and one of the Mc/Dc, bounding box or centre of gravity figure, depending on the appraoch used for the computation.
"PQ" - only Pc-Qc should be plotted.
"MD" - only Mc/Dc plot should be plotted.
"BB" - only bounding box plot should be plotted.
"COG" - only centre of gravity plot should be plotted.

Default is c("PQ", "MD", "BB" , "COG").

mdcol

vector of length 2 or NULL. If NULL colors in MD plot will be the same as in 'col' argument. If vector of length 2, first color stands for the Mc line and second color for the Dc line.

NULL - use color defined in 'col' argument.
c(4,3) - use blue for the Mc line and green for the Dc line.
c("firebrick", "cadetblue") - use of color names is also possible.

When used, it overrides 'col' argument. It is possible to set colors as numbers, or by the string name.

Default is NULL.

step

integer, set the step with which is plotted the bounding box or the cetre of gravity plot. Default is floor(length(x$pc)/20). This will plot 20 bounding boxes/centres of gravity.

col2

color for the second object in the plot if plotvar = "BB" (bounding boxes), or plotvar = "COG" (centre of gravity dots). Default is 2.

col

color for the lines in plots as defined in plot(). Default is 1.

main

string an overall title for the plot: see title

xlab

string a title for the x axis: see title

ylab

string a title for the y axis: see title

type

string what type of plot should be drawn: see plot

ylim

numeric vectors of length 2, giving the x and y coordinates ranges: see plot.window

...

arguments to be passed as graphical parameters.

Details

When plotvar = c("PQ", "MD"), or plotvar = c("MD", "PQ") the settings for main, xlab, ylab, ylim, would affect both plots, what does not make sense in most cases. To prevent this, setting of main, xlab, ylab and ylim only affects the first figure and second is set to default values for the given figure.

References

Gottwald G.A. and Melbourne I. (2004) On the implementation of the 0-1 Test for Chaos, SIAM J. Appl. Dyn. Syst., 8(1), 129<U+2013>145.

Martinovic T. (2019) Alternative approaches of evaluating the 0-1 test for chaos, Int. J. Comput. Math.

Examples

Run this code

# NOT RUN {
vec.x <- gen.logistic(mu = 3.55, iter = 2000)

# Output for each value of c
res2 <- testChaos01(vec.x, out = TRUE)

plot(res2[[1]], plotvar = c("PQ", "MD"), mdcol = c(4,3))
# }