Question

我有一个数据集，我的响应变量是离合器大小（TCL），两个预测因子是因素：位置（TYPE2）和年份（AN）。我通过交互进行了双向ANOVA分析。但是，在我的结果中，我无法看到交互的结果。数据看起来像这样（一个例子）：

是因为AN独立于LOC吗？（即，如果AN = 2000属于LOC = N，则意味着它永远不在LOC = S中）。每个AN只能在一个类别中，N或S.

这是我的模特和结果：

INIcrp <- aov(INI ~ TYPE2 * AN,data=todos)
summary(INIcrp)

#              Df Sum Sq Mean Sq F value Pr(>F)  
#TYPE2          1   3437    3437   380.9 <2e-16 *** 
#AN            14  29713    2122   235.2 <2e-16 ***

#Residuals   5168  46628       9

Answer 1

我不确定为什么缺少交互项，但有一些方法可以明确地添加它，例如：

aov(INI ~ TYPE2 * AN + I(TYPE2*AN),data=todos)

或

todos$TYPE2_AN <- with(todos,TYPE2*AN)
aov(INI ~ TYPE2 + AN + TYPE2_AN,data=todos)

在R中缺少相互作用的双向ANOVA

1 个答案: