0

私はRを初めて使用し、データに対してwilcox.testを使用しようとしています。probeIDとして行名を持つデータフレーム36021X246があり、最後の行は、サンプルが属するグループを示すラベルです。最初の140の「コントロール」と「最後の106の「扱われた」。

テストを実行するときに2つのグループを定義する方法を知っていただければ幸いです。...「式」の引数に関する多くの情報をオンラインで見つけることができません。ただし、-

"方式

lhs〜rhsの形式の式。ここで、lhsはデータ値を与える数値変数であり、rhsは対応するグループを与える2つのレベルの係数です。」

誰かがlhs〜rhsの意味と、この式を定義する方法を説明できれば、本当にありがたいです。

ありがとう!

4

1 に答える 1

1

R は通常、各行がケースであり、列が関連付けられた変数であると想定します。両方のサンプルのケースが同じデータ フレームで発生する場合、1 つの列がサンプル メンバーシップの指標変数になります。IndSample を呼び出しましょう。Wilcoxon は単変量検定であるため、検定対象の応答値を含む別の列が存在します。それを Y としましょう。

wilcox.test(y ~ IndSample, data=MyData, .....)

テスト用の残りのパラメーター: 両面ですか? 正確な統計が必要ですか?(おそらく、あなたの場合はそうではありません。)

私には、あなたのデータが横になっているように見えます。マトリックスの場合のように、データ フレームから行を引き出すことはできないため、これはデータ フレームでは問題になります。

最後の行を取得して、それを要因に変換する必要があります-次のようなものです

factor(c(MyData[lastrow,]))

次に、応答を含む行を引き出します。

Y <- as.numeric(c(MyData[ResponseRow,]))

次にウィルコクソンを行います。

しかし、私はあなたの状況を正しく理解しているかどうか確信が持てません。これは、控えめなウィルコクソン テストの非常に大きなデータ マトリックスのようです。

于 2013-02-08T20:38:48.627 に答える