# NOT RUN {
# data generation from '?prcomp', slightly modified
C <- chol(S <- toeplitz(.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