問題タブ [hypothesis-test]

このようなデータフレームがあります。2 つの個別のグループと 5 つの異なる変数があります。

T 検定を使用して 2 つのグループを区別するために、各変数の有意性を計算しています (以下を参照)。ただし、これは非常に小さなデータセットであるため、順列テストを実装して p 値を計算したいと思います。Rでこれを行うための最良の方法は何ですか?

私は現在、1 つの大きな行列に結合した複数のデータ セットに Wilcoxon ランクサム テストを実装しようとしています。Aつまり、ランクサム テスト時間705x17635を実行したいのです。 17,635for ループの使用はlapplyであり、次のように実行しました。

ここで、bはネガティブ コントロール データで、 は20000x1ベクトルです。ただし、これを実行するには非常に時間がかかります (私は 30 分後にあきらめました)。これを実行するためのより迅速な方法があるかどうか疑問に思っています。分ですが、さまざまな理由で R を使用する必要があります。

Rを使用して、40000行と4列のコンポストである区切られたテーブルからの入力を使用して多数を実行しています。iam t test ,p value を実装しようとしていますが、(データは本質的に一定である) というエラーがあります。for ループを使用し、両方のケースに適用します。コードは次のとおりです。

何が問題なのかわからない??

各係数の p 値 (有意性) を見つけるにはどうすればよいですか?

いくつかの変数を含むデータセットがあります:

X は数値変数、Y と Z は 2 つの因子のみを含む因子変数 (Y=1,2 Z=3,4)

因子変数が Y (t.test(X~Y)) の場合、t 検定は簡単に実行できます。しかし、たとえば、Z (3 と 4) の間で Y==2 の X 値のみを比較する t 検定を行う方法がわかりませんか?

私の言い方が正しいかどうか自信がないので、表で見た方が分かりやすいかもしれません。そこで、因子変数が Z で Y==2 である X のテストを行いたいと思います。どうすればこれを行うことができますか? STATA では簡単です: ttest var1 if var3==3, by(var2)

しかし、私はRでそれを取得しません:(

次のようなデータフレーム名を持つRのリストオブジェクトがあります。

等々 ...

次の平均値の差の t 検定を適用したいと思います。

についても同じ

等々...

リストには 84 個のデータフレームがあります。計算は、84 個のデータフレームすべてに対して上記のように実行されます。

ご提案いただきありがとうございます。歌。

「pav_DJF_histo.csv」および「pav_DJF...」という名前のすべてのデータフレームの再現可能な例。

各 df には 1 行 * 120 列があります。

私はやりたいです：

t.test(pav_DJF_histo.csv,all_dataframes)

R では多重回帰に対して返されることは知っていますが、βi=0 の仮説検定を返しますが、βi=1 のような検定をテストしたい場合はどうすればよいでしょうか。これには簡単なコマンドがありますか、そうでない場合は、係数の標準誤差、係数の値、回帰の自由度をどのように呼び出すので、t 分布 cdf を使用して p 値を計算できますか。一般的なプログラムで複数のデータを実行するためにこれを行いたい