3 方向の交互作用 (factorA、factorB、factorC) の 3 つのカテゴリ因子と、単純にモデルに追加される 4 番目の連続因子 (factorD) を持つ一般化線形モデルで分析した一連のデータがあります。モデルから一連の Tukey 文字グループ (つまり、コンパクトな文字表示) を取得しようとしていますが、相互作用をうまく含める方法が見つかりません。factorD を含めることには興味がありません。相互作用の 3 つだけです。
これでTukey調整されたペアごとの比較が得られました：
lsmeans(my.glm, factorA*factorB*factorC)
しかし、そこからコンパクトな文字表示を作成する方法がわかりませんでした。パッケージで実行multcompできますが、相互作用ではなく、そのパッケージの主効果でのみ実行する方法を見つけることができました。
それで、agricolaeこの投稿のようにパッケージを試してみました ( https://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukey-hsd-post-hoc-test-in -a-table-showing-groupe )は、それが機能するはずであることについて説明しています。ただし、その回答の指示に従うと、HSD.test から機能しない応答が返されました。具体的には、主効果テストを正常に機能させることができました。たとえばHSD.test(my.glm,"factorA")、交互作用を機能させることができませんでした。私はこれを試しました：
intxns<-with(my.data, interaction(factorA,factorB,factorC))
HSD.test(my.glm,"intxns",group=TRUE)
しかし、HSD.test 関数が "intxns" を有効なオブジェクトとして認識しなかったことを示すエラーが表示されます。私のGLMの）：
Name: inxtns
factorA factorB factorC factorD
HSD.test 関数呼び出しの factor フィールドにナンセンスを入力すると、同じエラーが発生します。inxtns オブジェクトをチェックしたところ、良さそうに見え、行数が residu の数と一致していました。agricolaeノートでは、実際には HSD.test での相互作用の使用については説明していませんが、機能すると思います。
HSD.test をインタラクションで動作させる方法を知っている人はいますか? または、相互作用のある glm 用のコンパクトな文字表示を作成するために機能するようになった他の機能はありますか?
私はこれに何日も取り組んできましたが、解決策を見つけることができませんでした。
ありがとう！