对emmeans比较结果的解释

Question

我有一个glm模型，有两个固定的效应，治疗和日期，从一个时间序列收集的数据来估计温度。在治疗过程中，有三种不同的类型: Fucus，terry布或Control，在这些檐篷下测量温度。该模型的创建方式类似于so mod1 <- glm(Temp ~ Treatment * Date, data = aveTerry.df )

我正试着判断一下，红衣是否有类似于Fucus冠层的效果(也就是复制它)。

我找到了emmeans软件包，相信它可以帮助我通过使用我的模型来比较治疗中的这些水平，并且使用它来找出估计的边际均值terry.emmeans <- emmeans(modAllTerry, poly ~ Treatment | Date)，并通过plot(terry.emmeans.average, comparison = TRUE) +theme_bw()绘制比较图。

给我这个输出链接在这里。

​

​

我正在寻找一些帮助来理解这个图形输出是什么，特别是，比较到底是什么(用红色箭头显示)。我有点理解蓝色盒子是一天内每一次处理的温度平均值的置信区间(基于模型)，但我想知道比较是如何进行的？和为什么有些日子只有单向箭头？

Answer 1

正如plot.emmGrid文档中所描述的那样，创建比较箭头的方式是两个箭头不相交的当且仅当它们各自的手段在指定的级别上有很大的不同。

集合中最低的平均值只有一个右指向箭头，因为这意味着不能与任何较小的箭头相比较，从而避免了对左箭头的需要。出于类似的原因，最高的平均值只有一个左箭头.这些箭头没有定义间隔；它们的唯一目的是描述比较。

在两两比较的be差异很大的情况下，可能无法构造比较箭头。如果发生这种情况，将显示一条错误消息。

置信区间也是可用的，但不应将这些置信区间用于比较手段。

更多的信息和例子可以通过vignette("comparisons", "emmeans")找到。此外，在vignette("xplanations", "emmeans")中给出了箭头实际构造的详细信息。