One way is to use the loadingplot() function in the package `pls': molprop.pc <- princomp(whatever) library(pls) loadingplot(molprop.pc, scatter = TRUE, labels = "names") (If you want comp 2 vs. comp 1: loadingplot(molprop.pc, comps = 2:1, scatter = TRUE, labels = "names") -- Bjørn-Helge Mevik