data(napping)
nappeplot(napping.don)
resindscal<- indscal(napping.don, napping.words)
x11()
prefpls(cbind(resindscal$points, napping.words))
x11()
pmfa(napping.don, napping.words, mean.conf = resindscal$points)
Run the code above in your browser using DataLab