如何在ggplot分组图表中添加p值？报错找不到对象mp_1怎么解决

你触发报错的核心原因是stat_pvalue_manual()函数无法将统计检验结果和ggplot的分组映射匹配：你在ggplot中指定了fill=mp_1作为亚组分组变量，但手动生成的统计检验结果stat.test3缺少和绘图映射匹配的位置、分组对应信息，导致函数寻找mp_1变量时匹配失败。

无需手动指定y坐标和x范围，直接用rstatix自带的add_xy_position()函数生成匹配绘图映射的位置参数即可，具体修改两处代码：

1. 调整统计检验代码，自动生成匹配位置 #

将原来的统计检验代码替换为以下内容，add_xy_position()会自动匹配你绘图时的x轴变量TBI、亚组dodge宽度0.9的设置，补全所有stat_pvalue_manual()需要的映射信息：

stat.test3 <- MP %>%
  group_by(TBI) %>%
  wilcox_test(age ~ mp_1) %>%
  adjust_pvalue(method = "bonferroni") %>%
  add_significance("p.adj") %>%
  add_xy_position(x = "TBI", dodge = 0.9)

2. 简化stat_pvalue_manual()调用 #

原来的调用参数配置错误，直接替换为以下代码即可，函数会自动识别匹配的分组和位置：

C2 + stat_pvalue_manual(stat.test3, label = "p.adj.signif")

如果你想显示具体校正后p值而非显著性星号，把label参数改为"p.adj"即可。

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