求助：无法从summary()/confint()提取指定配对对比对，代码返回NA

问题原因 #

你之前提取失败的核心原因是对比名称的格式不匹配——比如实际对比名称带有region前缀（和你的specs参数对应），或者你调用的对象结构不对：summary()返回的结果里，contrasts并不是可直接用$提取的命名列表元素，而是数据框中的一列。

正确提取步骤 #

1. 先确认目标对比的准确名称

首先得搞清楚系统生成的对比名称到底是什么，避免拼写/格式错误：

# 把对比结果转成数据框，查看所有对比名称
tuk_df <- as.data.frame(mN3tuk)
# 打印前几个对比名称确认格式
head(tuk_df$contrast)
# 或者直接搜索目标对比相关条目
grep("853|902", tuk_df$contrast, value = TRUE)

2. 方法一：转成数据框后筛选（最直观）

把mN3tuk转成普通数据框，就可以用常规的数据筛选方法提取目标对比：

# 转成数据框
tuk_df <- as.data.frame(mN3tuk)

# 提取单个对比（精确匹配名称）
target <- tuk_df[tuk_df$contrast == "region853 - region902", ]

# 如果要提取多个对比
targets <- tuk_df[tuk_df$contrast %in% c("region853 - region902", "region1573 - region1575"), ]

# 提取单个统计值（比如估计值）
estimate_val <- tuk_df$estimate[tuk_df$contrast == "region853 - region902"]
# 提取置信区间
ci <- tuk_df[, c("lower.CL", "upper.CL")][tuk_df$contrast == "region853 - region902", ]

3. 方法二：用emmeans自带的filter函数（更简洁）

emmeans的filter()函数可以直接对emmGrid对象（也就是mN3tuk）进行筛选：

# 提取单个对比
filter(mN3tuk, contrast == "region853 - region902")

# 提取多个对比
filter(mN3tuk, contrast %in% c("region853 - region902", "region1573 - region1575"))

4. 方法三：直接从emmGrid对象提取数值

如果你只需要某个统计量（比如p值、估计值），可以直接访问对象的内部元素：

# 提取目标对比的估计值
mN3tuk@estimates[names(mN3tuk@estimates) == "region853 - region902"]

# 提取目标对比的p值
mN3tuk@misc$p.values[names(mN3tuk@misc$p.values) == "region853 - region902"]

注意事项 #

• 对比名称必须完全匹配，包括前缀（比如region）、空格、大小写，建议直接从tuk_df$contrast里复制粘贴目标名称，避免手动输入错误。
• 如果你的region水平名称本身包含特殊字符，可能会导致对比名称格式异常，此时用grep()模糊匹配会更稳妥。

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