data(iris3)
tr <- sample(1:50, 25)
train <- rbind(iris3[tr,,1], iris3[tr,,2], iris3[tr,,3])
test <- rbind(iris3[-tr,,1], iris3[-tr,,2], iris3[-tr,,3])
cl <- factor(c(rep("s",25), rep("c",25), rep("v",25)))
z <- pcalda(train, cl)
pred <- predict(z, test)
## plot the projected data.
grpplot(pred$x, pred$class, main="PCALDA: Iris")
Run the code above in your browser using DataLab