Description Usage Arguments Details References See Also Examples
This function plot the Pc to Qc plot and Mc/Dc plot as described in Gottwald and Melbourne (2004).
1 2 3 4 |
x |
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.
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.
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 |
xlab |
string a title for the x axis: see |
ylab |
string a title for the y axis: see |
type |
string what type of plot should be drawn: see |
ylim |
numeric vectors of length 2, giving the x and y coordinates ranges: see |
... |
arguments to be passed as graphical parameters. |
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.
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–145.
Martinovic T. (2019) Alternative approaches of evaluating the 0-1 test for chaos, Int. J. Comput. Math.
1 2 3 4 5 6 | 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))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.