我想使用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及其互动{{ 1}}：
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
的ANCOVA吗？如果有的话，有没有办法在car::Anova没有任何互动的情况下指定此ANCOVA？
更新（2012年7月22日）：我在R-help上提出了这个问题，但到目前为止还没有回复。如果有新闻，我会在这里发布。