問題タブ [posthoc]

Rにパッケージ/関数が存在するかどうか疑問に思っています：「比率の均一性のカイ2乗検定のための事後ペアワイズ比較」（またはそれに相当するもの）ここで説明されています： http:// epm.sagepub.com/cgi/content/abstract/53/4/951

私の状況は、2×X行列でカイ検定を作成するだけです。違いを見つけましたが、どの列がその違いに「責任がある」かを知りたいです。

ありがとう、タル

いくつかのデータがあり、ANOVA と事後検定を行う場合、R の外部で図を編集するのではなく、事後分類を自動的に追加する箱ひげ図を作成するにはどうすればよいですか?

たとえば、開始するためのデータを次に示します。

簡単な一元配置 ANOVA とすべての予定外の比較事後テストを実行するためのコードを次に示します。

この時点で、x をグループ "a" に、y をグループ "b" に、z をグループ "a,b" に分類したいと思います。箱ひげ図は作成できますが、文字で注釈を付けるにはどうすればよいですか?

R で ANOVA と事後検定を使用する方法を理解しようとしています。これまでのところ、aov() と TukeyHSD() を使用してデータを分析しました。例：

私の問題は、ペアごとの比較の膨大なリストを持っているが、それを使って何もできないことです:

このデータセットには 40 行あります...理想的には、次のようなデータセットを取得したいと考えています。

• Act_Glu2 : a
• Act_Ala2 : a
• NegI_Ala2: バ...

要点を理解していただければ幸いです。これまでのところ、オンラインで匹敵するものは何も見つかりませんでした... TukeyHSD から生成されたファイルで重要なペアのみを選択しようとしましたが、ファイルは行と列で構成されていることを「認識」せず、選択が不可能です.. .

たぶん、私のアプローチに根本的な問題がありますか？

次の行を実行します。

しかし、MEANSステートメントはクラス変数で機能します（「エラー：この効果ではクラス変数のみが許可されています。」 SASは言います）。次に、応答 (試行) のみを定義しました: MODEL AA BB CC DD EE = / NOUNI;. 誰でも私を助けることができますか？

PS
(基本的に私のデザインはCrossValidatedで説明されています。)

3つのサイト（L、M、H）metで、4つの異なる植生タイプ（a、b、c、d）のパラメーター（）を測定する不均衡な実験があります。すべての植生タイプが3つのサイトすべてに存在します。植生タイプは、LとMで4回、Hで8回複製されます。

したがって、単純な分散分析とTukeyHSDは機能しません。パッケージAgricolae（HSD.test）とDTK（DTK.test）は一方向の設計でのみ機能し、次にmultcompがあります...mcp関数のTukey検定は、Tukey-Kramerの対比を計算しますか、それとも通常のTukeyの対比を与えますか？パッケージは不均衡な設計の多重比較をテストすることを目的としているため、最初のケースであると思いますが、両方のアプローチで生成されたp値は実質的に同じであるため、確信が持てません。その場合、どのテストが適切でしょうか？

また、不均衡なデータセットに対してこのような双方向分散分析を行うためのより適切なアプローチはありますか？

R のパッケージの関数bwplotから箱ひげ図に、ポストホック分析に続く有意性ラベルを追加したいと思います。フレームワーク内のソリューション。これを以前に試したことがあり、私にいくつかの指針を与えることができる人はいますか? bwplot オブジェクトから位置を導出する方法がある場合 ( R メーリング リストのこの記事で通常の boxplot オブジェクトに対して提案されているように)、その方向に目を向けるべきだと思いますか? 前もって感謝します！latticeggplot2lattice

敬具、

FM

私は3つのグループで prop.test() をやっています:

p.value は非常に低い (1.041e-13)。どのグループのペアが比率の有意差を示しているかを知りたいです。

次のコードを実行できると思いました....

...しかし、HSD.test は aov または lm オブジェクトしか処理できないため、機能していません。

どのような解決策が残っていますか? ソリューションが Tukey 補正ではなく Bonferroni (またはその他) の補正を意味するかどうかは関係ありません。

データセット (サンプルサイズが等しくない、分散が等しくない) を使用して、r で双方向分散分析を実行しました: 3 つの種で測定された 1 つの変数 (各種のオスとメス)。これは種間で有意な結果をもたらすため、どのペアワイズ比較が有意性を生み出すかを知りたい. Rで事後テストを実行するためのパッケージに関数があることを私は知っています：例えば

http://www.uwlax.edu/faculty/toribio/math305_fall09/multiple.txtからの Dunnett の事後検定。必要なパッケージ: 「multcomp」、「mvtnorm」、「survival」、「splines」

*注: glht は「multcomp」に記述されています

ただし、Dunett の検定は、すべてのグループを対照と比較するように設計されています。代わりに、すべてのグループを相互に比較したいと思います。Dunnett C です。Dunnett C を実行するパッケージを知っているか、それをコーディングする方法を知っている人はいますか? (式: http://pic.dhe.ibm.com/infocenter/spssstat/v21r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_posthoc_unequalvar.htm )