data(lipset)
results <- qmethod(lipset[[1]], nfactors = 3, rotation = "varimax")
title <- "Q method z-scores, lipset dataset"
subtitle <- paste0("Three factors, PCA, varimax. Printed on ",
Sys.Date())
plot(results, main = title, sub = subtitle)
Run the code above in your browser using DataLab