如何用marginaleffects获取三分类变量的全部两两比较结果？

默认情况下，marginaleffects::avg_comparisons()对分类变量采用treatment编码对比逻辑：仅返回每个水平与因子第一个水平（参考水平）的对比结果，因此三分类变量只会输出2组对比，而非全部3组两两组合。要获取所有两两比较结果，有两种简便方法：

方法1：使用pairwise快捷参数 #

直接通过comparisons参数指定变量采用两两对比逻辑：

avg_comparisons(
  lm1, 
  comparisons = list(species = "pairwise")
)

方法2：手动指定所有对比组合 #

如果需要自定义对比顺序或特定组合，可以手动列出所有两两对比项：

avg_comparisons(
  lm1,
  comparisons = list(species = c("Chinstrap - Adelie", "Gentoo - Adelie", "Gentoo - Chinstrap"))
)

两种方法都能输出和emmeans::pairwise一致的3组两两对比结果，包含Chinstrap与Gentoo之间的对比。

内容的提问来源于stack exchange，提问作者tomw