# data generation from '?prcomp', slightly modified
C <- chol(S <- toeplitz(0.9^(0:15)))
set.seed(17)
X <- matrix(rnorm(1600), 100, 16)
Z <- X %*% C
if (require("parameters") && require("psych")) {
pca <- principal_components(as.data.frame(Z), rotation = "varimax", n = 3)
pca
check_itemscale(pca)
}
Run the code above in your browser using DataLab