事後測試雙向混合模型ANOV

我正在做一個混合模型重複措施ANOVA。我想運行事後測試來查看交互TREAT * TIME的p值，但是我只設法使用下面的ghlt Tukey測試，它不會給我所尋找的交互。事後測試雙向混合模型ANOV

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")))

什麼，我要尋找的東西，如：

summary(glht(oi, linfct=mcp(TIME="Tukey",TREAT="Tukey",TREAT*TIME="Tukey")))

來源

2016-11-27 kumbu

使用snk.test(model, term="TREAT*TIME", among="TREAT", within="TIME")從包裝GAD如果你有從LSMEANS均衡模型和summary(lsmeans(oi, pairwise ~ TIME*TREAT), infer=TRUE)如果你的模式是不平衡

來源

2016-12-07 11:11:03 kumbu

我也有這個問題。似乎不存在雙向ANOVA的直接事後檢驗。但是，您可能想嘗試引導，這是一種雙向ANOVA的穩健估計形式。我發現以下鏈接非常有幫助。

http://rcompanion.org/rcompanion/d_08a.html

它包含使用rcompanion，WRS2，psych，並multcompView封裝來執行你的自舉方差分析和跟進事後一步一步的教程。祝你好運。

來源

2016-11-30 04:32:39

嗨！感謝您的回答。不過，我更深入地瞭解了這一點，並問了一些人，實際上這些功能是存在的。如果你有一個平衡模型和'summary（lsmeans（oi，pairwise）'，從包GAD中使用'snk.test（model，term =「TREAT * TIME」，在=「TREAT」中，在=「TIME」〜TIME * TREAT），推斷= TRUE）'如果你的模型不平衡，請享用！ – kumbu

@ kumbu，這太好了！我沒有問這個問題，但現在我有一個答案。 :) –

因爲你可以從AFEX包，而不是lme()找到與aov_ez()功能的替代混合模式，然後使用lsmeans()進行事後分析。

您將在這裏找到一個詳細的教程：

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

來源

2017-02-19 14:17:41

事後測試雙向混合模型ANOV

回答

相關問題