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

Question

У меня есть набор данных, где моя переменная ответа - это размер муфты (TCL) и два предиктора, которые являются факторами: местоположение (TYPE2) и year (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). Каждый может быть только в одной категории, либо 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)