Plots the kurtosis measures of an ics
object against its index number. Two versions of this screeplot are available.
1 2 3 4 
x 
object of class 
index 
index of the components to be plottes. If NULL all components are used. 
type 

main 
main title of the plot. 
ylab 
yaxis label. 
xlab 
xaxis label. 
names.arg 

labels 

... 
other arguments for the plotting functions. 
Klaus Nordhausen
plot.ics
, icsclass
and ics
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17  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)
screeplot(ics(data.normal), names.arg=paste("IC", 1:ncol(A), sep=""), xlab="")

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.