n <- 50 ; p <- 10
Xtrain <- matrix(rnorm(n * p), ncol = p)
Xtest <- Xtrain[1:5, ] + .4
fm <- pcaeigen(Xtrain, nlv = 5)
Ttrain <- fm$T
Ttest <- transform(fm, Xtest)
T <- rbind(Ttrain, Ttest)
group <- c(rep("Training", nrow(Ttrain)), rep("Test", nrow(Ttest)))
i <- 1
plotxy(T[, i:(i+1)], group = group,
pch = 16, zeroes = TRUE,
main = "PCA")
plotxy(T[, i:(i+1)], group = group,
pch = 16, zeroes = TRUE, asp = 1,
main = "PCA")
Run the code above in your browser using DataLab