We have a project on which we need to compare various methods for computing factor scores. Are there any R routines available that do parametric bootstrap or compute factor scores?