R中缺少交互的ANOVA双向答案

【问题标题】：ANOVA two-way with interaction missing in RR中缺少交互的ANOVA双向
【发布时间】：2015-01-16 02:47:30
【问题描述】：

我有一个数据集，其中我的响应变量是离合器尺寸 (TCL)，以及两个作为因子的预测变量：位置 (TYPE2) 和年份 (AN)。我做了一个带有交互作用的双向方差分析。但是，在我的结果中，我无法看到交互的结果。数据看起来像这样（一个例子）：

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

【问题讨论】：

你试过anova(lm(INI ~ TYPE2 * AN,data=todos))吗？
我无法用示例数据重现问题。

标签： r interaction anova

【解决方案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)

【讨论】：

对不起，我很难理解这里出了什么问题“这是因为年份不同”。你的意思是 AN 专栏？它们有何不同，导致您的结果中缺少交互？
我来到这里是因为我遇到了一个类似的问题，即交互项没有出现在嵌套模型中，并且想指出其他任何通过仔细检查使用 I(term1:术语 2)。我认为模型最初运行是因为 I(term1:term2) 出现在输出中，但结果与我只测试了 term1 一样。我最终遇到了一个不同的问题，导致交互不被包括在内，但只是想提醒任何来到这里的人仔细检查你使用 I() 的任何输出。