我想在下面的混合效果模型中可视化重要交互项的影响。
c1 <- lmer(log.weight ~ time + I(time^2) + temp + precip + time:precip + time:temp + (1|indiv), data = noctrl)
此模型包含“时间”的固定效果。 （简单和二次项），＆＃39;温度＆＃39;，＆＃39;降水＆＃39;和对数响应的两个相互作用＆＃39;权重＆＃39;所有术语都很重要，并且满足了模型对正态性和同质性的假设。
我一直在使用＆＃39;效果＆＃39;包以生成交互图以显示交互的效果。当我尝试用以下代码显示时间和温度（时间：温度）的相互作用时，我不确定结果图是否正确显示了这种相互作用。
ef2 <- effect(term="time:temp", c1, multiline=TRUE)
y <- as.data.frame(ef2)
ggplot(y , aes(time, fit, color=temp)) + geom_point() + geom_errorbar(aes(ymin=fit-se, ymax=fit+se), width=0.4)
我需要帮助理解resulting plot。为什么尽管这种相互作用的术语非常重要，但是在每个x的值上SE都会重叠？ 我正确使用效果包吗？或者是因为我需要在交互条件中包含时间I（时间^ 2）的二次项？
非常感谢您的澄清。