R-混合设计方差分析事后测试

时间:2018-08-18 11:37:14

标签: r anova mixed-models emmeans

我具有以下数据结构(带有示例值):

id     var1     var2   value
1      true     tr     1.34
2      true     ct     4.89
3      false    mm     2.38
4      true     tr     1.28

数据保存在'longData'中。因此,“ var1”介于可以为真或为假的主题变量之间,“ var2”是主题因子范围内的3个级别(tr,ct,mm),“值”是数字值。

我已经做出了如下混合设计方差分析:

anovaResult = ezANOVA(data=longData, 
                                dv=.("value"), 
                                wid=.("id"),
                                within=.("var2"),
                                between=.("var1"),
                                type=3)

结果显示var1和var2之间有显着的相互作用。现在,我想进一步研究这种交互,但是我不知道如何。我听说过emmeans软件包(估计的边际均值似乎是此处选择的统计数据,因为我是统计学的新手,请随时向我提出建议),但无法执行命令。这可能是因为我是R的新手,并且不完全理解语法。

有人能给我提供一个如何测试两个因素之间相互作用的有效例子吗?我也不会拒绝解释如何解释结果。

我知道这是我所要求的,但是我无法自己解决这个问题,并且不得不很快花时间介绍结果,而又没有太多时间来学习统计学和R。 谢谢。

1 个答案:

答案 0 :(得分:1)

这将有助于提供示例数据集。

但是,您可以运行Tukey测试:

Private Sub add_Click()
Sheets("Block B").Range("C8").Select
ActiveCell.EntireRow.Insert Shift:=xlDown Sheets("Block B").Range("C8:L8").Select  Selection.Borders.Weight = xlThin

或者通过交互在Anova上运行emmeans:

mod1<-aov(value~Factor1*Factor2, df)
TukeyHSD(mod1)

或者做一个混合模型,这似乎是您要做的:

mod1<-aov(value~Factor1*Factor2, df)
library(emmeans)
emmeans(mod1, pairwise~Treatment*Time)