data(rhone)
pca1 <- dudi.pca(rhone$tab, scan = FALSE, nf = 3)
iv1 <- pcaiv(pca1, rhone$disch, scan = FALSE)
iv1
# iner inercum inerC inercumC ratio R2 lambda
# 6.27 6.27 5.52 5.52 0.879 0.671 3.7
# 4.14 10.4 4.74 10.3 0.984 0.747 3.54
plot(iv1)
Run the code above in your browser using DataLab