MASS里面就可以，配套书籍《Modern applied statistics with s 》中的R代码：
ability.FA <- factanal(covmat = ability.cov, factors = 1)
> ability.FA
更多信息参考该书，或者其他多元统计书籍（with R）。

SPSS

去找spss书

可见附件。