r - car::Anova拥有不与主题内因素相互作用的协变量的方式

Question

我想使用car::Anova运行ANCOVA，但无法找出是否有一种方法可以仅将协变量添加为主要效果(即，不应与任何事物交互)。

据我了解的ANCOVA，协变量只是添加到模型中的另一个主要效果(即另一个效果)，从而控制了该协变量的总体累加影响。随后，协变量不与其他因素相互作用。但是，我无法在Anova中添加一个不会与主题内因素相互作用的变量(即，我的最终模型似乎没有ANCOVA)。

让我用?Anova中的示例来说明我的问题。 OBrienKaiser数据集在(treatment和gender)因数之间有2个，在(phase和hour)因数内有2个。现在假设我们还记录了参与者的age，并希望将其添加为任何分析的协变量。

require(car)
set.seed(1)

n.OBrienKaiser <- within(OBrienKaiser, age <- sample(18:35, size = 16, replace = TRUE))

# the next part is taken from ?Anova
# I only modified the mod.ok <- ... call by adding + age
phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)), levels=c("pretest", "posttest", "followup"))
hour <- ordered(rep(1:5, 3))
idata <- data.frame(phase, hour)

mod.ok <- lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5, post.1, post.2, post.3, post.4, post.5, 
          fup.1, fup.2, fup.3, fup.4, fup.5) ~  treatment*gender + age, data=n.OBrienKaiser)
(av.ok <- Anova(mod.ok, idata=idata, idesign=~phase*hour, type = 3)) 

结果表明，结果包含与协变量 age的交互作用，即对象内部(或重复测量)因素 phase， hour及其相互作用 phase:hour的交互作用:

Type III Repeated Measures MANOVA Tests: Pillai test statistic
                            Df test stat approx F num Df den Df Pr(>F)  
(Intercept)                  1     0.129     1.33      1      9  0.278  
treatment                    2     0.443     3.58      2      9  0.072 .
gender                       1     0.305     3.95      1      9  0.078 .
age                          1     0.054     0.52      1      9  0.490  
treatment:gender             2     0.222     1.28      2      9  0.323  
phase                        1     0.418     2.87      2      8  0.115  
treatment:phase              2     0.871     3.47      4     18  0.029 *
gender:phase                 1     0.084     0.37      2      8  0.703  
age:phase                    1     0.393     2.59      2      8  0.136  
treatment:gender:phase       2     0.545     1.69      4     18  0.197  
hour                         1     0.565     1.95      4      6  0.222  
treatment:hour               2     0.580     0.72      8     14  0.676  
gender:hour                  1     0.310     0.68      4      6  0.633  
age:hour                     1     0.508     1.55      4      6  0.301  
treatment:gender:hour        2     0.707     0.96      8     14  0.504  
phase:hour                   1     0.975     9.56      8      2  0.098 .
treatment:phase:hour         2     1.145     0.50     16      6  0.873  
gender:phase:hour            1     0.693     0.56      8      2  0.770  
age:phase:hour               1     0.974     9.40      8      2  0.100 .
treatment:gender:phase:hour  2     1.314     0.72     16      6  0.723  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

我的问题是: 可以使用car::Anova运行ANCOVA吗？如果可以，是否可以在不与age进行任何交互的情况下指定此ANCOVA？

更新(2012年7月22日):我在 R-help上问了这个问题，但到目前为止没有任何答复。如果有新闻，我会在这里发布。

Answer 1

我在R-help上问了这个问题，该问题开始了与John Fox(后来由Peter Dalgaard参加)的有益讨论。不幸的是，它分为两个线程:one，two。

重点是:
“对象内部的对比度是由Anova()构造的，在设计的行基础上是正交的，因此您应该能够安全地忽略对您(由于某种原因使我逃避)不感兴趣的效果。” (约翰·福克斯)

因此，该问题的答案是:没有人不能，但是没关系，因为这些交互作用不会改变其他影响，因为它们是正交的。