数据是关于10只股票的超额收益率，给数据拟合三个因子的统计因子模型，如何基于所拟合的模型得到这10只股票的相关矩阵
S-plus中的命令是这样的
stat.fac=factanal(rtn,factor=3,method=’mle’)   
  corr.fit=fitted(stat.fac)
请教各位懂S-plus的，麻烦解释一下这个fitted(stat.fac)得到的值？在R 里面应该怎样做？