問題タブ [posthoc]

溶けた のサブグループでテューキーのテストを実行していdata.frameます。結果には、CLASS と comparison の 2 つの列でグループ化されたペアごとの比較が表示されます。私は現在multcompLetters、コンパクトな文字表示を取得するために使用しようとしていますが、この機能は 2 つのグループ化列の存在に苦労しているようで、機能させることができません。助言がありますか？

非正規応答変数に Anova を適用し、残差が正規かどうかを検証し、事後分析を使用したいのですが、どのように作成できますか?

データフレームの砂丘 (ページの下) に3 つの列があり、3 つの異なる砂丘生態系について記録されたマラムグラスの被覆率を説明しています。

(1) 復元された; (2)劣化。(3) ナチュラル。

エコシステム間の有意差を確立するために、テスト 1 とテスト 2 という 2 つの異なる一元配置分散分析テスト (以下) を実行しました。テスト 1 は、生態系間の大きな違いを明確に示しています。ただし、テスト 2 では有意差は見られません。ボックス プロット (以下) は、生態系間の差異の著しい違いを示しています。

その後、データフレームを溶かして、応答変数でもある階乗列 (つまり、Ecosystem.Type という見出し) を作成しました。アイデアは、glm モデル (テスト 3 - 以下) を適用して、一方向 Anova でテストすることです。ただし、この方法は失敗しました (以下のエラー メッセージを参照してください)。

問題

各一元配置分散分析テストを実行するコードが正しいかどうか、および有意に異なる生態系のペアを区別するための事後テスト (トルコ HSD、シェッフェなど) を実行する正しい手順であるかどうか、私は混乱しています。どなたかお役に立てることがありましたら、アドバイスをいただければ幸いです。どうもありがとう....

テスト 1

事後テスト

テスト 2

テスト 3

溶けたデータフレーム

データ

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

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

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

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

混合モデルで反復測定分散分析を行っています。事後検定を実行して相互作用 TREAT*TIME の p 値を確認したいのですが、探している相互作用が得られない次の ghlt Tukey テストしか使用できませんでした。

私が探しているのは次のようなものです：

線形混合モデルの事後テストを実行したいと考えています。私は最初に を使用しましたがlsmeans∷lsmeans、これは正常に動作します。問題は次lsmeansのとおりです。少なくとも私が知っていることではありません。そのため、グループを手動で並べ替えました。私の実験計画には多くの処理があり、多くの事後テストを実行するため、これはなんとなく苦痛です (そしてエラーが発生しやすくなります)。

私は を見てきましたagricolae∷HSD.test。これは、事後テスト用のグループを提供します。欠点は、HSD.test混合モデルでは実行されないことです ( multcomp::cld、multcomp::hsd、またはmultcomp::multcompLettersそのこともありません)。

編集：答え