data(octane)
octane <- octane[, -1] # remove the dependent variable y
pca=PcaHubert(octane, k=10)
screeplot(pca, type="lines")
pca2 <- PcaHubert(octane, k=2)
plot(pca2, id.n.sd=6)
pca7 <- PcaHubert(octane, k=7)
plot(pca7, id.n.sd=6)
Run the code above in your browser using DataLab