All functions extract the indicated matrix from the fitted model, and
will work with any object having a suitably named component.
The default scores and loadings methods also handle
prcomp objects (their scores and loadings components are called
x and rotation, resp.), and add an attribute
"explvar" with the variance explained by each component, if
this is available. (See explvar for details.)