aov对象的类和错误术语是R中的事后测试的问题

Question

我在R中运行以下代码作为重复测量anova分析的一部分：

model1 = aov(GazeCueTarget.RT ~ ProbeState * RecallType * CueType 
             + Error(Subject/(ProbeState * RecallType * CueType)), edata.cue.only)

这三个因素anova有两个ProbeState级别，两个用于RecallType，两个用于CueType。我遇到的问题是，我的学术主管告诉我包括一个错误术语（见上文），以说明存在的任何主题内错误。然而，当我这样做时，我遇到了进一步测试的问题。我可以生成结果摘要，但无法运行leveneTest(model1)或

posthoc = TukeyHSD(x=model1, 'ProbeState', conf.level=0.95)

这是一个问题，因为我发现了ProbeState的主要影响，但我无法区分显着差异所在的位置。我认为这个问题源于这样一个事实，即在我的模型中包含错误项会改变R anova对象的类。如果我运行class(model1)，我可以看到输出显示为“aovlist”和“listof”。然而，如果我运行模型1不包含错误项的class(model1)，那么我得到“aov”和“lm”。在TukeyHSD()的帮助部分中，它表示该函数将用于x，它通常是具有aov拟合的模型对象。

我的问题是我是否在anova中错误地写了我的错误术语？我应该使用另一个事后测试吗？我想继续使用anova。我不想将分析更改为回归/线性混合模型。

非常感谢任何明确的帮助。