問題タブ [psych]

パッケージのdescribe関数を使用しようとしています。psychただし、次のエラーが発生します。

describe(ToothGrowth) のエラー: description は少なくとも長さ 1 の文字列でなければなりません

パッケージを再インストールしてロードしdata(ToothGrowth)、データセット ライブラリから例を添付するようにしましたが、それでもこのエラーが発生します。

psych パッケージの principal() 関数が $score 要素を計算する方法を理解したいと思います。

相関行列ではなく共分散行列を試してみたい。

基本的に、PCAのスコアは、重みとしてローディングマトリックスを使用して、元の中心データの線形結合である必要があるため、試しました：

principal()関数が負荷に何らかのスケーリングを使用することは理解していますが、比率は各列で同じである必要がありますが、ここではtest.

相関行列を使用すれば、これは問題になりません。例えば：

ヘルプ ドキュメントでは、因子分析の用語が使用されているため、さらに混乱しました。誰かがここで私を啓発できることを願っています。

前もって感謝します！

R 言語とパッケージは初めてです。行列形式で約 9000 の遺伝子のペアワイズ ピアソン相関分析を行うために、ここのリンクからの情報に従って R で psych パッケージを使用しました。

しかし、サイケマニュアルでは解決できない分析上の問題に直面しています。

1 つ目: 一般的なエラー "Error in cor(x, use = use, method = method) : 'x' must be numeric" . 要素名を削除して値のみを保持すると、機能します。ヘッダーも含めるにはどうすればよいですか?次のコードは上記のエラーを示しています

私の 2 番目の疑問: ピアソン相関 >=0.5 を持つペアをフィルタリングする最良の方法は何ですか? 私はそれを個別に行うべきですか、それともR自体にメソッドがありますか?

編集：

有意性検定を使用cor.ciしてポリコリック相関を取得しようとしていますが、エラー メッセージが表示され続けます。コードは次のとおりです。

これが与えるエラーメッセージは次のとおりです。

mydata を認識させたり、分析のためにこのデータセットから特定の変数を選択したりするにはどうすればよいですか? ここから上記のコードを取得しました：
Rで有意なポリコリック相関行列

ありがとう！

RStudioでknitrを使ってレポートを書いています。また、因子分析に psych パッケージを使用します。

このコード チャンクを実行すると、エラーが発生します。

TeX コードを保存すると、LaTeX では正常に動作しますが、knitr では動作しません。

これを解決する方法について何か提案はありますか?

スタック コミュニティの皆様、こんにちは。

私の質問はコマンドに関するものです

パッケージからerror.bars.by{psych}。

コマンドの仕組みに精通しています：

私の質問は、このコマンドが 2 つのグループ化変数を処理できるかどうかです。

もしそうなら、どのようにコーディングしますか？

データテーブル

PCA ローディング と の 2 つの行列がloa.origあり、それが (手動またはその他の方法で) の回転loa.rotであることを知っているとします。loa.rotloa.orig

(また、バリマックスなどによってすでにloa.orig直交回転されている可能性もありますが、それは問題ではないと思います)。

が回転して に達する角度を知りたいのはわかっています。loa.origloa.rot

この別の質問に対するコメントから、「回転は通勤しない」ことを理解しているため、ペアワイズ (プレーンワイズ) 回転の順序も重要です。

したがって、から再現loa.rotするloa.origには、一連の必要な回転を、理想的にrotsは以下の順序で知る必要があります。

MWE は次のとおりです。

明らかに、この場合、私は角度と順序を知っていますが、知らないと仮定しましょう。また、実際の使用例では、3 つだけでなく、多くのコンポーネントが保持およびローテーションされる可能性があると仮定しましょう。

コンポーネント ペアごとの (平面ごとの) 回転角度の順序を報告する従来の方法がどうなるかは少しわかりませんが、考えられる組み合わせのリスト(~~順列ではなく~~) を使用する必要があると思います。

rotsこれらの順列を与える.

loa.rotfromloa.origに到達する方法を知っておくとよいでしょうrots。

更新:以下の回答に基づく試み

以下の回答に基づいて、関数をまとめて、varimax回転と実際のデータセットでテストしようとしました。(特に理由はありません。実際に角度がわからないvarimax回転が欲しかっただけです。)

次に、抽出した角度を使用して、バニラの負荷からバリマックス回転を実際に再作成できるかどうかをテストします。

残念ながら、それらは同じではなく、似ていません:(

明らかに、私は間違いを犯しています。