質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.51%

  • Python 3.x

    6361questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

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

受付中

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 221

ak_miyamoto

score 23

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

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

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

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

左端に観測変数がありますが、省略してあります

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

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

メッシュコードごとに、因子得点の計算を回帰法によって推定する。
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公式ライブラリ

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

0

Factor analysis in R and Python – Kohei Watanabe

によると、入力のスケーリングが必要なようです。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.51%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

  • 解決済

    PyQtにおけるQWidgetとQGraphicsViewの組み合わせ

    http://d.hatena.ne.jp/mFumi/20141112/1415806010 このサイトをお手本にQWidgetとQGraphicsViewの組み合わせを実

  • 受付中

    windows上でmecabとpython3をバインディングさせる

    前提・実現したいこと windows上でmecabとpython3をバインディングさせたい 発生している問題・エラーメッセージ こちらの方の記事を参考に( http://

  • 解決済

    vim-flake8を導入した際エラー行の色を変えたい

    実現したいこと macのデフォルトターミナル(novel)上でvimを開き、 nvie/vim-flake8のプラグインを導入しました。 F7を押すとpep8違反となる行を表

  • 解決済

    pythonの図示 matpotlib

    pythonで四角形が格子状に区切られそれぞれの点に数字が振られています。それをコンターマップのように(ただし、等高線はいらない)図示したいのですがどのように書けばよいでしょうか

  • 解決済

    Pythonのmatplotlibを使って三次元散布図をシンプルに描写したい

    前提・実現したいこと Python初心者です。 1人3属性からなる1000件ほどのデータを、3次元の散布図として描写したいです。 属性一つ一つがそのままX,Y,Zに対応するイメージ

  • 解決済

    python エラー

    weights, params = [], [] for c in np.arange(-5, 5): lr = LogisticRegression(C=10**c,

  • 解決済

    python-shogiの使い方

    python3(anaconda)にpython-shogiモジュールを入れて shogi.CSAを使いたいのですがエラーがでて使えません。どのようにすればよいのでしょうか sho

  • 解決済

    経路問題を図示[python]

    前提・実現したいこと 以下の画像のように経路問題を図示したいです (番号はついてなくてもいいです) 辺は以下のように格納されています #全部の辺 lines = set(

同じタグがついた質問を見る

  • Python 3.x

    6361questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。