Plot for an epplab Object

Description

The function offers three informative plots for an epplab object.

Usage

## S3 method for class 'epplab'
plot(x, type = "kernel", angles = "radiants",
  kernel = "biweight", which = 1:10, as.table = TRUE, ...)

Arguments

Details

The option which can restrict the output to certain simulation runs. In case of many simulations, this might improve the readability.

For type="kernel", the default, it plots a kernel density estimate for each of the chosen directions. In the case of type="histogram" the corresponding histograms. For type="angles" it plots the angles of the first chosen direction against all others. Whether the angles are given in degrees or radiants, depends on the value of angles.

Author(s)

Daniel Fischer, Klaus Nordhausen

See Also

xyplot, densityplot, histogram, density

Examples

library(tourr)
data(olive)
res <- EPPlab(olive[,3:10],PPalg="PSO",PPindex="KurtosisMin",n.simu=10, maxiter=20)

# Plot with kernel estimator
plot(res)

# Just the best 5 and then 8
plot(res,which=c(1:5,8))

# Plot as histogram
plot(res,type="histogram")

# Plot an angles plot
plot(res,type="angles")

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(REPPlab)
Loading required package: rJava
Loading required package: lattice
Loading required package: LDRTools
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/REPPlab/plot.epplab.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.epplab
> ### Title: Plot for an epplab Object
> ### Aliases: plot,epplab-method plot-method plot.epplab
> ### Keywords: hplot methods
> 
> ### ** Examples
> 
> 
> library(tourr)
> data(olive)
> res <- EPPlab(olive[,3:10],PPalg="PSO",PPindex="KurtosisMin",n.simu=10, maxiter=20)
Simulation 0... finished (I -5.009891687 in 0.644s)
Simulation 1... finished (I -22.431740100 in 0.181s)
Simulation 2... finished (I -29.306060952 in 0.040s)
Simulation 3... finished (I -2.600885740 in 0.023s)
Simulation 4... finished (I -9.894486466 in 0.167s)
Simulation 5... finished (I -25.094953519 in 0.157s)
Simulation 6... finished (I -5.762172036 in 0.027s)
Simulation 7... finished (I -9.840642653 in 0.026s)
Simulation 8... finished (I -46.866192671 in 0.053s)
Simulation 9... finished (I -21.358160792 in 0.026s)
Warning message:
In EpplabOutputConv(jepplab, maxiter) :
  There were 10 non-converged simulation runs!
> 
> # Plot with kernel estimator
> plot(res)
> 
> # Just the best 5 and then 8
> plot(res,which=c(1:5,8))
> 
> # Plot as histogram
> plot(res,type="histogram")
> 
> # Plot an angles plot
> plot(res,type="angles")
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>