在双因素方差分析的上下文中，我想在另一个因子（SpeciesA，SpeciesB）的水平内输入一个因子（10,20,30度）的对比度水平。这假设相互作用效果显着，因此不能简单地对温度主效应进行对比。
我试图在EZanova，multcomp和phia中做到这一点，但一直无法弄清楚有效的语法。我错过了允许这个的套餐吗？
对比度和系数
对比1：对于物种A，10度是否与20和30不同？ 种类：Temp（2，-1，-1,0,0,0） Contrast2：物种B的10度与20度和30度不同吗？ 种类：Temp（0,0,0,2，-1，-1） 对比3：物种A的10度是否与物种B的20度和30度不同？ 种类：Temp（2,0,0,0，-1，-1）
species<-rep(c("speciesA","speciesB"),each=12)
temp<-rep(rep(c("10","20","30"),each=4),2)
y<-rnorm(24,5,3)
(result<-anova(lm(y~species*temp)))