問題タブ [statistical-test]

R で複数の Wilcoxon テストを実行することについて質問があります。7 つのデータセットがあり、それぞれについて、10 の分類子を使用して 9 つの異なる特徴抽出方法を比較する必要があります。それらは確率論的な方法であるため、各分類器に 10 個の値/繰り返しがあります。単一のデータセットの簡単な例として、3 つのメソッド、3 つの分類子、3 つの繰り返しがあるとします。

ここで、異なる方法、つまり、M1.C1 対 M2.C1、C2、C3、M1.C1 対 M3.C1、C2、C3 など、これらの分類器の各ペア間で Wilcoxon 検定を実行する方法が必要です。 . 私はRに比較的慣れていないので、誰かがデータをロードして複数のテストを簡単な方法で実行する適切な方法を提案できますか?

私は主成分分析を実行していますが、スケーリング 1 のベクトルは長さ 1 であるはずだと言われました。ここでは、それらは 1 よりも非常に大きくなっています。スケーリング 2 では、1 未満であると想定されています。

私は何か間違ったことをしていますか？

また、データを標準化して、スケーリング 2 で記述子の平衡円を描くことは可能ですか?

データは、さまざまな種の鳥のくちばしの測定値です。

データセットは次のようになります (名前は pca.bird.clean):

PCA:

これは、画像を生成するためのコードです (こちらを参照)

完全に再現可能な例を次に示します。

対応する平日 (1 = 月曜日、2 = 火曜日、7 = 日曜日など) の数値 (0 の場合もあります) があるデータがあります。さまざまな統計テストによって、自分の値 (たとえば、気温など) が平日に依存するかどうか (月曜日の値が大きく、日曜日の値が小さいなど) を調べたいと考えています。私はシミュレートしたデータを持っていますが、測定もしています。この可能性のある依存関係をクラスカル・ウォリス検定でテストするにはどうすればよいですか?

そして第二に; この表では、平日または「気温」が私の従属変数であり、それらはどのような性質（順序、カテゴリなど）ですか？

データ:

プロポーションの違いを比較するためにマラスキーリオ法を行っています。私は次のコードを使用しています（このチュートリアルからコピーして適応させました：

出力にカテゴリのラベルも出力する必要があります (たとえば、どのカテゴリが正確に比較されているか)。

したがって、カテゴリが分離されたベクトルでリストされている場合、たとえばcategories <- c("cat1", "cat2", cat"3)、比較はcat1-cat2、cat1-cat3、およびcat2-cat3です。

これらのラベルを出力に追加するにはどうすればよいですか?

(相対的な) 少量の繰り返し回数 = 10,100 で完全に機能する for ループがあります。しかし、「時間」の値が大きい場合、行列を埋めることでエラーが発生します: 下付き文字が範囲外です... (以下のコードとコードの説明を参照)

エラー： Error in M_zp_var[j, (1:n)] : subscript out of bounds

行列を埋めるために、最初の for ループで 2 番目の for ループを使用します

そして私もそれを試しました

代わりに、それも機能しません。

私が言ったように、私のouter-for-loopが小さい場合、コードは機能します。変数の「倍」を大きく（1000、5000）選択すると、なぜ異なる動作をする必要があるのか​​ わかりません。

「範囲外の添字」の意味を知っています

あなたが私を助けてくれることを願っています。:(

コード：

説明：

ループは、順列テスト (モンテカルロ) を「回」回シミュレートします。1 つのループは 2 つのランダム サンプルを生成し、それらを B=16000 回ランダムに並べ替えます (並べ替えは関数 で行われますZ_pi())。次に、異なる並べ替え ( unique()) のみを取得し、各並べ替えからテスト統計を計算します... コードは停止します関数を含むマークされた行 (太字) で作業しますZ_pi()。