GPoM (version 1.4)

visuOutGP: visuOutGP : get a quick information of gPoMo output

Description

The algorithm aims to get a quick information about the outputs obtained with gPoMo.

Usage

visuOutGP(
  ogp,
  selecmod = NULL,
  id = 1,
  prioMinMax = "data",
  opt3D = "TRUE",
  maxPages = NULL,
  seeEq = 1
)

Value

A Matrix describing the terms composing each model by row. The first row corresponds to the model detection (1 unclarified, 2 diverging, 0 is fixed point, -n with n an integer, is period-n cycle' )

Arguments

ogp

The output list obtained from gPoMo.

selecmod

A vector of the selected model. Maximum 24 models can be presented at the same time.

id

The type of model to identify. id = 1 corresponds to the unidentified models, that is, potentialy chaotic models).

prioMinMax

Gives the priority for the plots among: "data", "model", "dataonly" and "modelonly".

opt3D

Provides a 3D plot (x,y,z) when opt = 'TRUE' (the rgl library is required).

maxPages

The maximum of pages to be displayed (4 by default, but this may be insufficient when too many models remain)

seeEq

Indicates if equations should be displayed (seeEq = 1, by default) or not (seeEq = 0).

Author

Sylvain Mangiarotti

Examples

# load data
data("Ross76")
# # time vector
tin <- Ross76[seq(1, 3000, by = 8), 1]
# single time series
data <- Ross76[seq(1, 3000, by = 8), 3]
dev.new()
plot(tin, data, type = 'l', main = 'Observed time series')
# global modelling
# results are put in list outputGPoM
outputGPoM <- gPoMo(data, tin=tin, dMax = 2, nS=c(3), show = 0,
                    nPmin = 9, nPmax = 12, method = 'rk4',
                    IstepMin = 200, IstepMax = 201)
visuOutGP(outputGPoM)