data(lipset)
results <- qmethod(lipset, nfactors=3, nstat=33,
nqsorts=9, 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