問題タブ [chi-squared]

次のようなデータフレームがあります。

したがって、2 つのグループがあります: alleleA と alleleB について、確率を観察し、予測しました。

このデータに対してカイ 2 乗検定を実行し、2 つのグループの期待確率と観測確率を比較するように言われました。多くのチュートリアルを検索しましたが、カイ二乗検定が 2 つのグループでどのように機能するかを理解できませんでした。明確にするために、私は人々が私のためにコードを投稿したくありません. 私はただ確信が持てないので

1. Web 上の既存のカイ 2 乗検定の例を見ても、自分のデータに対してカイ 2 乗検定を実行する方法を理解できませんでした。

2. 2 つのグループ間でまとめてテストを実行する例は見つかりませんでした。

3. データフレームの各行に対してカイ二乗検定を実行する例が見つかりませんでした。

各行のカイ二乗値と p 値が必要です。どこを見ても、彼らはデータフレーム全体に対してカイ二乗テストを実行しています。

この問題を解決するにはどうすればよいか、お好きな方法でリンクを共有するか、洞察を共有してください。しかし、私があなたに私のコードを書いてほしいと思い込まないでください。

ありがとう！

与えられたさまざまなデータに対して、いわゆる Ljung Box テストを実行しようとしています。データがカイ 2 乗分布の 0.95 分位以内にあるかどうかを確認したいと考えています。この値は、表 ( http://www.unc.edu/~farkouh/usefull/chi.htmlなど) で見つけることができます。データのサンプル サイズがさまざまであるため、カイ 2 乗分布の自由度が異なり、すべての値を手動で調べるには多くの時間がかかります。

これらの値を自動的に見つけて、for ループでこれらの異なる値を使用する方法が matlab にあるかどうか疑問に思っていますか?

こんにちは、何かが欠けていると本当に思いますが、データフレームでchisqを計算する方法がわかりません。テストのようなデータフレームがあり、chisq を計算して、カテゴリの頻度がセット全体と比較して高いかどうかを判断します。また、すべてのカテゴリで実行したいと考えています。

ご協力いただきありがとうございます！

28 df が 7.5 までの x であるカイ 2 乗分布の密度関数のグラフを R でプロットしようとしています。

今まで、私が集めてきたものからこれを得ました：

しかし、プロットは機能していないようです._.

全て、

各行が 7 つの関心領域 (AoI) の 1 つの領域から他の各領域への視線の集計を表す遷移行列で、目のデータに対して chisq.test() を実行しようとしています。この分析では、ある AoI からそれ自体への移行があることは意味がありません。したがって、これらのフィールドには NA が含まれます。

8 列と行の基本的な表形式の入力 (一番上の行がヘッダーで、左の列が "from") から単純な 3 列のデータ from (from、to、values) まで、さまざまな形式を試しました。 ）。

私のdata.frameは次のようになります：

(これは実際のデータではなく、Holmqvist 氏らのアイ トラッキングに関する教科書から取得したサンプル データであることに注意してください。)

これをマトリックスとしても試したことに注意してください

しかし、問題はデータの形式ではなくNAにあると思いますが、そうであれば、それをどのように処理するかわかりません。

親愛なるプログラミングの神様へ

R で一連のカイ 2 乗検定 (私の種の存在/不在 data.frame の列ごとに 1 つの検定) を実行したいと思います。種 (列名)、カイ 2 乗検定統計量、df、および p.value を出力します。

私の種のデータ スニペット (実際のサイズ = 50x131):

次の適用関数コードは、登った岩と登っていない岩 (Env$Climbed) での特定の種の出現回数を含む分割表を最初に作成することにより、種 (列) ごとにカイ 2 乗検定を実行します。

これにより、種 (列) ごとに個別の data.frame が生成されます。各種の列名も含む 1 つの data.frame を作成したいと思います。このようなもの：

これは ddply または adply で行うべきですか? それともただのループ？（私は試しましたが、失敗しました）。同様のトピックに関する投稿を確認しました ([ R で for ループを使用したカイ 2 乗分析) が、私の目的のために機能させることができませんでした。

あなたの時間と専門知識をありがとう！TC

私は R の初心者です。プロジェクトの必要性により、10 万エントリの Chisq テストを行う必要があります。

私は数日間独学で学び、ループで chisq.test を実行するためのコードをいくつか書きました。コード:

このコードにはいくつかの問題がある可能性があります。しかし、それは機能します。

ただし、実行速度は非常に遅いです。

「適用」を使用して改善しようとします

「for」を使用する代わりに、apply を 2 回使用する予定です

ただし、行列が関数ではないというエラーがあります。zsd chisq.test の出力はリストです。write.table を使用してデータを出力することはできません。

the.data はこんな感じです。

私は昼夜を問わず努力してきました。誰かが私を助けてくれることを願っています。どうもありがとうございました。

カイ二乗検定のコードを MATLAB で作成しました。0.897 や 0.287 などの P 値を取得したいのですが、結果が小さすぎます。以下は私のコードです：

同様の結果でADテストを使用してみました：

以下は、ワイブル密度関数を当てはめたデータのヒストグラムです (ワイブル パラメーターはA=4.0420とですB=2.0853) 。

1 つのコピー領域に 4 つの遺伝子型の観測カウントのサンプルがあります。私がやりたいことは、これらの遺伝子型の対立遺伝子頻度を計算し、これらの頻度のテストが、R のカイ 2 乗を使用して 25%:25%:25%:25% の期待値から大きく外れることです。

これまでのところ、私は得ました：

次に、合計数を取得します。

現在の周波数:

ここで私は今迷っています。af1、af2、af3、af4 が 0.25、0.25、0.25、0.25 から大きく外れているかどうかを知りたい

Rでこれを行うにはどうすればよいですか?

ありがとう、エイドリアン

編集：

申し分なく、私は示唆されているように chisq.test() を試しています:

警告メッセージは何を伝えようとしていますか? 近似が正しくないのはなぜですか?

この方法論をテストするために、予想される 0.25 からかけ離れた値を選びました。

この場合、値が予想される 0.25 値からかなり離れていても、H0 はまだ拒否されません。

カイ二乗統計を取得するために scipy.stats.chi2_contingency メソッドを使用しています。パラメータとして度数表、つまり分割表を渡す必要があります。しかし、特徴ベクトルがあり、頻度表を自動的に生成したいと考えています。そのような機能はありますか？私は現在このようにやっています：

データ シリーズとターゲット シリーズは列の値で、他の 2 つはインジケーターの名前です。誰でも助けることができますか？ありがとう