ANOVA двусторонний с отсутствующим взаимодействием в R

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)