R中缺少交互的ANOVA双向

Question

我有一个数据集，其中我的响应变量是离合器尺寸（TCL），并且有两个预测因子：位置（TYPE2）和年份（AN）。 我使用交互进行了双向ANOVA分析。 但是，在我的结果中，我看不到交互的结果。 数据如下所示（示例）：

TCL
8
6
5
3
2
3
3
3
NA
4

LOC
N
S
S
N
N
N
N
S
N
N

AN
2000
2005
2005
2010
2010
2000
2010
2008
2000
2000

是因为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)