3

を使用して ANCOVA を実行したいのcar::Anovaですが、主効果としてのみ共変量を追加する方法があるかどうかを確認できません (つまり、何とも相互作用しないでください)。

私が ANCOVA を理解している限り、共変量はモデルに追加されたもう 1 つの主効果 (つまり、もう 1 つの効果) であり、それによってこの共変量の全体的な相加的影響を制御します。したがって、共変量は他の要因と相互作用しません。ただし、Anova被験者内要因と相互作用しない変数を追加することはできません (つまり、私の最終モデルは ANCOVA に適合していないようです)。

の例で私の問題を説明しましょう?Anova。OBrienKaiserデータセットには、( と ) の間に 2 個、(およびtreatment)gender内に 2個phaseありますhour。ここ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

私の質問は:で ANCOVA を実行できますか?car::Anovaもしそうなら、この ANCOVA を相互作用なしで指定する方法はありますageか?

更新(2012 年 7 月 22 日): R-helpでこの質問をしましたが、今のところ回答がありません。ニュースがあれば、ここに投稿します。

4

1 に答える 1

1

R-help でこの質問をしたところ、John Fox (後に Peter Dalgaard が参加) との有益な議論が始まりました。残念ながら、2 つのスレッド ( one、two )に分割されました。

オチは次のとおりです。無関心。」(ジョン・フォックス)

したがって、質問に対する答えは次のとおりです。誰もできませんが、これらの相互作用は直交しているため、他の効果を変更しないため、問題ではありません。

于 2012-08-10T13:45:50.617 に答える