問題タブ [tukey]

私はアンバランスなデザイン/サンプルを扱っており、もともと学んだaov(). lm()ANOVA テストでは、 を使用するのではなく、 を使用してフィッティングを使用するタイプ III 平方和を使用する必要があることがわかりましたaov()。

問題は、 を使用して事後テスト (特に Tukey の HSD) を取得することlm()です。私が行ったすべての調査simintでは、multcompパッケージでの使用は機能すると言われていましたが、更新されたので、そのコマンドは利用できないようです. aov()また、計算するために通過することに依存しているようです。

基本的に、R について私が見つけたすべての Tukey HSD テストaov()は、lm(). タイプ III の二乗和を取得するには、不均衡な設計に必要な以下を使用する必要があります。

を使用して mod で Tukey HSD テストを使用するにはどうすればよいlm()ですか? または逆に、Type III を使用して ANOVA を計算しても、Tukey HSD テストを実行できますか?

ありがとう！

タイプ III エラーと種と生息地内のペアワイズ比較を使用しながら、ホームレンジのサイズに対する種と生息地の影響を比較したいデータセットがあります。
データのサブセットを次に示します。

現在、私は 3 つの種のデータを分割し、それぞれについて別々の ANOVA を実行していますが、1 つの ANOVA で同時に種と生息地について尋ねる方が理にかなっていると思います。ある種に対して実行した ANOVA の例を次に示します。

aov() はタイプ I エラーを使用しているようです。. . 適切ではないと思います。さらに、テューキーの検定はペアワイズ比較のアプローチとして保守的すぎる可能性があると思います。ホームレンジに対する種と生息地の両方の影響を考慮した1つのANOVAを実行できるアプローチを誰かが助けてくれますか?

StackOverflow で関連する質問への回答を読んでいるときに、matplotlibでウィスカーの位置を計算し、外れ値を検出するために使用されるコードを見ました。

これで、このelse部分は完全に理にかなっています。Tukey boxplotsの仕様に従って、上位四分位数の 1.5 IQR 内で最高のデータを見つけます。実際、それは-以下max(wish_hi)にある最大のデータ エントリです。 Q3+1.5*IQR

しかし、そのor部分は...私には理解できません。if len(wisk_hi) == 0を翻訳すると...

この条件はどのように適用できますか? Q3 は、中央値でデータを分割し、上半分の中央値を取り、その上に 1.5*IQR を追加することによって検出されます。この値よりも低いデータが存在しないのはなぜですか?

これが空のデータセットに関するものである場合、の 2 番目の部分もor意味がありません (Q3 または IQR はデータなしでは意味がないため)。

おそらく明らかな何かが欠けている - 助けて？

次のような csv ファイルが複数あります。

実際には次のようになります: CSV テーブル

各 .csv ファイルは、この場合「logtPSA.csv」などの化学パラメーターにちなんで命名されます。ご覧のとおり、ヘッダーは「ライブラリ」と「値」です。FDA、SMMRNA、VD_SM、および VD_MV の 4 つの異なるライブラリがあります。

私はそれらすべてに対してANOVAテストとTukeyHSDテストを実行してきましたが、Rは結果がセルに分割される形式で提供しないため、Excelにコピーして貼り付けるだけでよいため、手動で行う場合はかなり退屈になります。

すべてのファイルを一度に 1 つずつ調べてテストを実行し、すべてを 2 つのきれいなテーブル (ANOVA 用に 1 つ、TukeyHSD 用に 1 つ) に出力する関数を作成できる方法はないかと考えていました。

Tukey HSD テストの後に同種のサブセットを表示することについての私の質問は、過去に既に質問されています ( https://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukeyを参照してください)。 -hsd-post-hoc-test-in-a-table-showing-groupe、議論の最後に）しかし、答えられたことはありません。私は現在同じ問題を抱えており、このようなグループを表示する R で使用できる関数はこの関数だけであるように思われるため、解決策がどうしても必要です。そうでない場合は、教えてください。要するに、これは私が得るものです：

Groups, Treatments and means a 140095-001 36.79 b 150004-001 32.1 b 136936-021 31.97 bc 137219-004 31.39 bc 136673-017 31.27 bc 136963-009 30.79 bcd 147328-016 30.63 bcd 0147592-01 30.55 bcde 140094-001 30.02 cde 136730-007 29.7 de 136963-066 29.49 ef 136936-004 28.4 efg 147414-004 28.2 efg 137109-036 28.2 efg 136765-001 28.06 efg 140089-001 27.82 efg 137186-020 27.8 fg 136936-006 27.48 fgh 147350-014 27.43 gh 136992-001 27.36 gh 136730-015 27.18 ghi 0147785-01 27.08 ghi 0147691-01 26.98 ghi 136891-010 26.7 ghij 0147792-01 26.49 ghijk 136947-014 26.3 ghijkl 140097-001 25.8

そして、これは私が欲しいものです:

私はすでに関数の背後にあるスクリプトをチェックしましたが、私の知識がかなり低いため、何も見つけることができませんでした。少なくとも、まっすぐで明確なものは何も見つかりませんでした。誰かが解決策を手伝ってくれますか? お答えいただければ幸いです。

前もって感謝します、ジェル

R には、4 つの数値ベクトルがあります。これらの 4 つのベクトルの平均が同じかどうかをテストしたいと思います。そうでない場合は、どの平均が大きく、どの平均が小さいかをテストします (基本的に、それらの平均で 4 つのベクトルをランク付けしたいと思います)。

この投稿は私に非常に適しています：

http://brownmath.com/stat/anova1

私はこのガイドに従っています：

https://www.r-bloggers.com/analysis-of-variance-anova-for-multiple-comparisons/

ただし、この投稿は 4 つの平均が等しいかどうかのテストで終わりますが、答えが等しくない場合、どの平均が R で大きく、どの平均が小さいかをどのように知ることができますか?

どうもありがとうございました、

投稿のコードは次のとおりです。

次のデータがあります（dat）

上記のデータ セットに対して TukeysHSD ペアワイズ テストを実行したいと考えています。テストの結果から、有意な比較をグラフに組み込みたいと思います (有意なグループ間に「*」または「**」記号を表示)。

これは試みられたコードです:

各結果を外部配列として保存してからggplot2に呼び出すことなく、テストの結果（「重要性」）をグラフに直接組み込む方法（できるかどうか）と方法がわかりません。

次のデータがあります（dat）

次のデータがあります（dat）

上記のデータ セットに対して TukeysHSD ペアワイズ テストを実行したいと考えています。

ただし、実行しようとするたびに、実行が継続され、出力が得られません。これを修正する方法について何か提案はありますか?

また、関数 TukeyHSD() を使用して同じテストを実行したいと思います。ただし、ワイド/ロング形式を使用しようとすると、次のようなエラーが発生します