问题描述
我有一个具有两个固定影响的glm模型，即“处理”和“日期”，可以根据时间序列中收集的数据估算温度。在“处理”中，有三种不同的类别：悬钩子，毛圈布或对照，并且在这些遮篷下测量温度。这样创建模型mod1 <- glm(Temp ~ Treatment * Date,data = aveTerry.df )
我正试图说明Terrycloth是否具有与Fucus冠层相似的效果（即复制它）。
我找到了emmeans软件包，并相信它可以通过使用我的模型来帮助我比较治疗中的这些水平，并以此来找到估计的边际均值terry.emmeans <- emmeans(modAllTerry,poly ~ Treatment | Date)并绘制比较结果通过plot(terry.emmeans.average,comparison = TRUE) +theme_bw()
我正在寻找有助于理解此图形输出是什么的信息，尤其是比较究竟是什么（由红色箭头显示）。我有些理解，蓝框是一天中每种处理的温度平均值的置信区间（基于模型），但是想知道如何进行比较吗？ 为什么有些日子里只有一个箭头？？
解决方法
如plot.emmGrid的文档中所述，比较箭头的创建方式是：当且仅当两个箭头各自的均值在所述水平上显着不同时，两个箭头才脱节。
该集合中的最低均值仅具有一个向右箭头，因为该均值将不会与任何较小的均值进行比较，从而无需使用向左箭头。由于类似的原因，最高均值只有一个向左的箭头。这些箭头没有定义间隔。他们的唯一目的是描绘比较。
在成对比较的SE差异很大的情况下，可能无法构造比较箭头。如果发生这种情况，则会显示错误消息。
也可以使用置信区间，但不应将这些配置项用于比较均值。
更多信息和示例可通过vignette("comparisons","emmeans")找到。另外，vignette("xplanations","emmeans")