首页 文章

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

提问于
浏览
0

我有一个数据集,我的响应变量是离合器大小(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

1 回答

  • 0

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

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

    要么

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

相关问题