data(abr1)
cl <- factor(abr1$fact$class)
dat <- abr1$pos
model <- pcalda(dat,cl)
## Second component versus first
plot(model,dimen=c(1,2),main = "Training data",ep=2)
## Pairwise scatterplots of several components
plot(model,main = "Training data",ep=1)
## The first component
plot(model,dimen=c(1),main = "Training data")
Run the code above in your browser using DataLab