R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Screeplot for an ICS Object

screeplot.ics

R Documentation

Screeplot for an ICS Object

Description

Plots the kurtosis measures of an ics object against its index number. Two versions of this screeplot are available.

Usage

## S3 method for class 'ics'
screeplot(x, index = NULL, type = "barplot", 
          main = deparse(substitute(x)), ylab = "generalized kurtosis", 
          xlab = "component", ...)

Arguments

`x`	object of class `ics`.
`index`	index of the components to be plottes. If NULL all components are used.
`type`	`barplot` if a barplot or `lines` if a line plot is preferred.
`main`	main title of the plot.
`ylab`	y-axis label.
`xlab`	x-axis label.
`...`	other arguments for the plotting functions.

Author(s)

Klaus Nordhausen

Examples

set.seed(654321)
A <- matrix(c(3,2,1,2,4,-0.5,1,-0.5,2),ncol=3)
eigen.A <- eigen(A)
sqrt.A <- eigen.A$vectors %*% (diag(eigen.A$values))^0.5 %*% t(eigen.A$vectors)
normal.ic <- cbind(rnorm(800), rnorm(800), rnorm(800))
mix.ic <- cbind(rt(800,4), rnorm(800), runif(800,-2,2))

data.normal <- normal.ic %*% t(sqrt.A)
data.mix <- mix.ic %*% t(sqrt.A)

par(mfrow=c(1,2))
screeplot(ics(data.normal))
screeplot(ics(data.mix), type="lines")
par(mfrow=c(1,1))
rm(.Random.seed)

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(ICS)
Loading required package: mvtnorm
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ICS/screeplot.ics.Rd_%03d_medium.png", width=480, height=480)
> ### Name: screeplot.ics
> ### Title: Screeplot for an ICS Object
> ### Aliases: screeplot.ics
> ### Keywords: hplot
> 
> ### ** Examples
> 
> set.seed(654321)
> A <- matrix(c(3,2,1,2,4,-0.5,1,-0.5,2),ncol=3)
> eigen.A <- eigen(A)
> sqrt.A <- eigen.A$vectors %*% (diag(eigen.A$values))^0.5 %*% t(eigen.A$vectors)
> normal.ic <- cbind(rnorm(800), rnorm(800), rnorm(800))
> mix.ic <- cbind(rt(800,4), rnorm(800), runif(800,-2,2))
> 
> data.normal <- normal.ic %*% t(sqrt.A)
> data.mix <- mix.ic %*% t(sqrt.A)
> 
> par(mfrow=c(1,2))
> screeplot(ics(data.normal))
> screeplot(ics(data.mix), type="lines")
> par(mfrow=c(1,1))
> rm(.Random.seed)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>