雙向混合模型的事後tukey測試ANOVA

Question

我和我的一些學生在許多地方搜索了一個解決方案，在沒有運氣的字面上，幾個月。我一直參考我不想使用的lme命令。所提供的輸出不是我的同事或我自己使用超過15年的。此外，由於我使用R作爲教學工具，因此它不會像t檢驗一樣流暢，並且對於介紹統計學生的學生來說也不是單向的anovas。我正在進行雙因素RM ANOVA，其中一個因素重複。我成功地讓R複製了Sigmaplot爲主要效果提供的內容。然而，由R給出的事後分析與Sigmaplot中的相同事後分析顯着不同。這裏是我使用的代碼 - 與筆記（因爲我用這也教學生）。

#IV between: IVB1 - Independent variable - between subject factor 
#IV within: IVW1 - Independent variable - within subject factor 
#DV: DV   - Dependent variable. 

aov1= aov(DV ~ IVB1*IVW1 + Error(Subject/IVW1)+(IVB1), data=objectL) 
summary(aov1) 

# post hoc analysis 
ph1=TukeyHSD(aov(DV ~ IVB1*IVW1, data=objectL)) 
ph1 

我希望有人能幫忙。 謝謝！

Answer 1

我也有這個問題，我覺得從AFEX包，而不是AOV（在aov_ez（）函數方便的選擇），然後進行使用LSMEANS事後分析（）而不是TukeyHSD（）：

model <- aov_ez(data,       
      id="SubjID",      
      dv="DV",       
      within=c("IVW1", "IVW2"), 
      between = "IVB1") 

# Post hoc 
comp = lsmeans(model,specs = ~ IVB1: IVW1: IVW2, adjust = "tukey") 
contrast(comp,method="pairwise") 

你會發現一個詳細的教程在這裏：

https://www.psychologie.uni-heidelberg.de/ae/meth/team/mertens/blog/anova_in_r_made_easy.nb.html