Python scikit learn: Factor_Analysisで因子負荷量を計算したい

環境　python 3.6
初学者です。scikitlearnのFactor_analysisを用いて因子分析を行いたいと考えておりますが、以下の2点でつまってしまっています。

1
因子負荷量を計算したいと思うのですが、こちらの例にのっとって以下のように計算したところ

fa = FactorAnalysis()
factors = fa.fit_transform(matrix)

結果としては絶対値１をこえるものが多くありました。

因子負荷量の定義上、おそらくは上の計算で因子負荷量でないものが計算されてしまっていると思うのですが、因子負荷量を計算するにはどのようにしたらいいのかわかりません。教えて下さると幸いです。

2
１で算出した因子負荷量を利用して因子スコアを計算したいと思うのですが、どのように計算したらよいでしょうか。以下は自分で書いてみたコードですが、何かよいライブラリ・書き方がありましたら教えて下さると幸いです。

メッシュコードごとに、因子得点の計算を回帰法によって推定する。
F = A.T ・R(-1)・ X
F：因子得点
A：因子負荷(転置行列使用)
R：相関係数行列(逆行列使用)
X：元データ行列

AT = np.array(factors).T
Rinv = np.linalg.inv(np.array(corr_mat))
X = np.array(data).T
tmp = np.dot(AT,Rinv)
score = np.dot(tmp,X)

参考：scikit-learn FactorAnalys公式ライブラリ