一个小例子，希望可以帮到你——  

data(USArrests)
fa <- factanal(~., factors = 1, data = USArrests, score = "Bartlett", rotation = "none") #不做旋转  
D <- diag (fa $ uniquenesses)   #特殊方差
  
A <- as.matrix(fa$loadings[,1])  #载荷矩阵

D1 <- solve(D)  #D的逆矩阵
x <- t(as.matrix(USArrests))

#Bartlett方法（最小二乘法）
  
因子得分为f <- solve(t(A) %*% D1 %*% A) %*% t(A) %*% D1 %*% x  
  
然后标准化就是fa$scores了
fa <- factanal(~., factors = 1, data = USArrests, score = "regression", rotation = "none")
r <- fa$correlation #x的相关矩阵
Thompson方法（回归方法）
f <- t(A) %*% solve(r) %*% x

可以使用网页版SPSSAU进行分析 结果简单方便操作