問題タブ [factor-analysis]

R でいくつかの Stata の結果を複製しようとしていますが、多くの問題が発生しています。具体的には、Stata が探索的因子分析で行うのと同じ固有値を回復したいと考えています。具体的な例を挙げると、factorStata のヘルプはbg2データ (医師の費用に関するもの) を使用し、次の結果を示します。

表の最初の列の固有値に興味があります。R で同じデータを使用すると、次の結果が得られます。

ご覧のとおり、これらの値は Stata の結果とはかなり異なります。2 つのプログラムは固有値の計算に異なる手段を使用している可能性がありますが、 R コマンドfaのオプションのほとんど (すべてではないにしても) を含め、さまざまな固有値の抽出方法を試しました。他のいくつかの R コマンド。Stata と同じ固有値を抽出することはできません。また、Stata のマニュアルを読んで、Stata が使用する方法を正確に把握しようとしましたが、十分に具体的に把握できませんでした。factanalprincipal

どんな助けでもいいです！質問に答えるために追加情報が必要な場合はお知らせください。

次のコードを使用してloadingsマトリックスを抽出しています。

どちらが機能しますか！

しかし、私が同じことをしようとするとcorrelation

私は得る：AttributeError: 'R' object has no attribute 'correlation'

もし - 代わりに - 私が試してみます:

私は得る：ValueError: All parameters must be of type Sexp_Type,or Python int/long, float, bool, or None

特に、相関と負荷の両方が有効である必要があるため、これは奇妙だと思います。

私に与えます：

次のコードでbiplot生成された赤い線を取り除こうとしています。誰かが助けてくれれば幸いです。

R の psych パッケージから fa.parallel() 関数を呼び出すときに問題が発生します。

スペースの都合上、データセットの先頭のみを示していますが、私の問題を理解していただければ幸いです。

次のサンプル データセットがあるとします。

欠損値には 1 ～ 5 のカテゴリと NA があります。

私が実行すると

コンソール出力には次のように表示されます。

上の表でわかるように、値 1 と 2 のみを持つ「SQ112」のような変数、または応答として 1 と NA を持つ「SQ109」のような変数があります。

ただし、これは観察されたデータに過ぎず、すべての変数で 1 ～ 5 または NA の応答選択肢が可能ですが、その場合は単に観察されません。

したがって、グローバルをTRUEに設定したいのですが、関数パラメータでTrueを指定しても、関数を呼び出すたびにグローバルがfalseに設定されます。

因子分析を実行するための統計factanalパッケージの関数を使用しています。

次のことを試しました。

Divide_item_Factor_Loading(fit) を実行しようとすると、エラーと呼ばれる

浮き出る。

Divide_item_Factor_Loading の私の完全なコード:

このエラーを回避するには、どのような手順を実行する必要がありますか?

mirtパッケージを使用して CFA を実行しようとしています。FA、FB、FC の 3 つの実質因子と、F1、F2 などのいくつかの一般因子があります。

実質的な要因を相互に相関させたいと思いますが、残りの要因を実質的な要因とは相関していないが、相互に相関しているものとしてモデル化します。

3+1 因子については、次のモデルを使用しました

これは完全に正常に機能し、次の summary() になります。

さらに因子 (F2、F3、...) を追加し、それらを「自由に相関させる」ように指定すると、同じ結果が得られません。モデル

次の summary() が生成されます。

COV定式化が機能しない理由について何か考えはありますか?

どうもありがとう！KH