私は繰り返しモデルanovaを混合モデルで使っています。ポストホックテストを実行して、TREAT * TIMEインタラクションのp値を確認したいのですが、私が探しているインタラクションを与えてくれない次のghlt Tukeyテストを使用することしかできませんでした。二方向混合モデルのためのポストホックテストanova
summary(glht(oi, linfct=mcp(TIME="Tukey",TREAT="Tukey",TREAT*TIME="Tukey")))
使用snk.test(model, term="TREAT*TIME", among="TREAT", within="TIME")モデルがあなたの代わりにlme()のAFEXパッケージからaov_ez()機能を持つ代替品を見つけることができる混合モデルについて
私もこの問題を抱えている：私が探していられるもの
library(multcomp)
library(nlme)
oi<-lme(total ~ TREAT * TIME, data=TURN, random = ~1|NO_UNIT)
anova(oi)
summary(glht(oi, linfct=mcp(TIME="Tukey", TREAT="Tukey")))
のようなものです。 2方向ANOVAのためのストレートフォローホックテストは存在しないようです。 しかし、二方向ANOVAのための強力な推定の一形態であるブートストラッピングを試してみることができます。私は非常に有用な次のリンクを発見した。
http://rcompanion.org/rcompanion/d_08a.html
は、それはあなたのブートストラップANOVAを行い、事後のフォローアップするrcompanion、WRS2、psych、およびmultcompViewパッケージを使用してステップバイステップのチュートリアルが含まれています。がんばろう。パッケージGADから
アンバランスである場合は、LS平均からバランスの取れたモデルとsummary(lsmeans(oi, pairwise ~ TIME*TREAT), infer=TRUE)を持っている場合は、 lsmeans()を使用してポストホック分析を実行します。
あなたはここに詳細なチュートリアルを見つける：
https://www.psychologie.uni-heidelberg.de/ae/meth/team/mertens/blog/anova_in_r_made_easy.nb.html
こんにちは！ご回答有難うございます。しかし、私は少し深くそれを見て、いくつかの人々に尋ね、実際にそれらの機能が存在します。 バランスの取れたモデルと 'summary（lsmeans（oi、pairwise）がある場合は、GADパッケージの' snk.test（model、term = "TREAT * TIME"、= "TREAT" 〜TIME * TREAT）、infer = TRUE） 'を返します。楽しい！ – kumbu
@kumbu、これは素晴らしいです！私はこの質問をしなかったが、今私は答えがある。 :) –