問題タブ [hypothesis-test]

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 値を計算できますか。一般的なプログラムで複数のデータを実行するためにこれを行いたい

pandas/statsmodels で t 検定を実行して 2 つのグループ間のパフォーマンスの違いを比較しようとしていますが、statsmodels が (合理的な方法で) 使用できる方法でデータをフォーマットするのが困難です。

私のパンダのデータフレームは現在次のようになっています。

そして、t検定を実行するには、次のように治療ごとに整理されたデータが必要であることを理解しています。

このコードはほとんどトリックを行います:

しかし、印刷すると、パフォーマンス値ではなく、その処理が発生したインデックスを取得しているように見えます。

[おそらく?] 次のようなものにする必要があります。

このデータフレームを t 検定を実行できる形式に変換する最良の方法は何ですか?

次の文を読みました。

機能的 MRI データは、サンプル数に比べて高次元です (通常、1000 サンプルに対して 50000 ボクセル)。この設定では、機械学習アルゴリズムのパフォーマンスが低下する可能性があります。ただし、単純な統計テストは、ボクセルの数を減らすのに役立ちます。

スチューデントの t 検定 (scipy.stats.ttest_ind) は、2 つの分布が統計的に異なるかどうかを判断する単純な統計検定を実行します。2 つの異なる条件でボクセルの時系列を比較するために使用できます (この場合、家または顔が表示されている場合)。時系列分布が 2 つの条件で類似している場合、ボクセルは条件を識別するのにあまり関心がありません。

この検定は、2 つの時系列が同じ分布から引き出される確率を表す p 値を返します。p 値が低いほど、ボクセルの識別性が高くなります。

から: http://nilearn.github.io/building_blocks/manipulating_mr_images.html

この t 検定は 4 つのクラス (条件) にも適用できますか。

利用可能なこれのMatlab実装はありますか?

次のように、マトリックスに 94 個の変数 (サンプル + タンパク質 + グループ) と 172 個の観測値があります。

サンプルの一部はグループ 0 にあり、一部はグループ 1 にあります。t 検定を使用してグループ 0 と 1 の間に違いがあるかどうかをテストしたいのですが、すべてのタンパク質に対してそれを実行したいと考えています。申請書を使おうと思ったのですが、使い方がわかりません。また、名前は Protein1, protein2... ではありません。かなり長いので、すべてを記述する必要はありません。

また、マトリックス内の各タンパク質の p 値のみが必要です。次のようなものです。

または、p 値が 0.05/92 未満のものだけを見つけることができるように、似たようなものを作成します。

編集：

長い形式で作業を開始しましたが、これはもう問題ではありません:

したがって、ここで概説されているように、KS テストを実行して、dtwo データセットの分布に違いがあるかどうかを評価できます。

では、以下のデータを取ってみましょう

次に、var3 の結果に基づいて分離できます

これで、コルモゴロフ–スミルノフ検定を実行して、個々の変数の分布の違いを検定できます。

当然のことながら、違いは得られず、異なるデータセットが同じ分布から引き出されたと想定できます。これは、次の方法で視覚化できます。

しかし、ここで、 &の両方を一緒 に考えたときにvar3==0、 との間の分布がvar3==1異なるかどうかを検討したいと思います。複数の予測子がある場合にそのようなテストを実行する R パッケージはありますかvar1var2

同様の質問がここに提起されましたが、回答はありません

いくつかの文献があるようです: 例 1 例 2

しかし、Rにリンクされているものは何もないようです

SPSS は、カスタム テーブルの z 検定に関しては、複数回答セットをカテゴリ変数とは異なる方法で扱います。この動作は、応答の重複に関連していると思いますが、その方法がわかりません。

では、多重応答セット (MRset) に関して、SPSS はどのように z 検定を行うのでしょうか?

私の目標は、R で MRsets の SPSS z-test を再現することですが、SPSS が実際に何をするのかわかりません。通常、SPSS カスタム テーブルの z-testing は、

しかし、明らかに MRsets の場合は異なります。

これを明確にするために、このカテゴリ対 MRset の比較を見てください。

カテゴリ変量 z 検定 (C 列と D 列は z 検定によると違いはありません)

• カテゴリー データセット (重複なし、3623 件):データセットのダウンロード
• カテゴリ重複行列 (重複なし):

• カテゴリカル z 検定 SPSS 構文

  /li>

• カテゴリー z 検定の出力:

• C<->D z-test の R 再現 (最初の行): http://www.r-fiddle.org/#/fiddle?id=p4gw9ftk

  /li>

MRset z-test (表の数値は同じですが、z-test の結果は異なります: C 列と D 列に有意差があります)

• MRset データセット (重複を含む、2404 件):データセットのダウンロード
• MRset オーバーラップ マトリックス:
• MRset z テスト出力:

• MRset z-test SPSS 構文:

  /li>

• C<->D z-test の R 再現 (最初の行): http://www.r-fiddle.org/#/fiddle?id=GAhnnrv0

  /li>

MRset R コードからわかるように、オーバーラップ ケースを差し引いても役に立ちません。多分それは重み付けか何かにリンクされていますか？アイデアをどうもありがとう。

おそらく役立つリンク:重みと多重応答セットに関する注意事項