分组小提琴图/箱线图的趋势线添加、刻度设置及p值标注问题

嘿，我来帮你搞定这三个ggplot2的问题，每个问题都给你讲清楚原因和解决办法，最后附上完整的修正代码：

1. 移除图例中的透明度（alpha）部分 #

你之前把alpha=0.5放在了全局的aes()里，这会让ggplot把透明度当成一个数据映射美学，从而自动在图例里生成alpha调节滑块。解决方法很简单：把alpha从美学映射里移出来，直接作为geom_violin()和geom_boxplot()的图层参数（因为这只是设置图层的固定透明度，不是映射数据）。如果之后图例还残留alpha项，再加一句guides(alpha = "none")强制隐藏。

修改示例：

# 把alpha从全局aes移到图层参数里
ggplot(df,aes(x=group_age,y=value,fill=age,color=age)) + 
  geom_violin(alpha=0.5) + 
  geom_boxplot(width=0.1,aes(middle=mean(value)), alpha=0.5) + 
  # 其他图层...
  guides(alpha = "none")  # 保险起见，强制隐藏alpha图例

2. 将x轴刻度设置为group并居中 #

要实现每个group对应一个居中的刻度，我们需要调整x轴的映射逻辑：不再用group_age（group+age的组合）作为x，而是直接用group作为x，然后通过position_dodge()让同一个group下的不同age的小提琴图/箱线图并排显示，这样x轴自然就是group的刻度，自动居中（哪怕某个group缺失age类别，这个逻辑依然成立）。

另外注意：你原来的geom_smooth()是基于离散的group_age拟合，这没有意义，我们需要用连续的time变量来拟合时间趋势，所以要在geom_smooth()的aes里手动指定x=time。

修改后的核心代码：

ggplot(df,aes(x=group,y=value,fill=age,color=age)) + 
  # 用position_dodge让同一group下的age分组错开
  geom_violin(position = position_dodge(width = 0.8), alpha=0.5) + 
  geom_boxplot(width=0.1, aes(middle=mean(value)), 
               position = position_dodge(width = 0.8), alpha=0.5) + 
  # 基于连续的time变量拟合，按group分组
  geom_smooth(aes(x=time, group=group), color="black", method='lm', size=1, se=T) + 
  theme_minimal()

3. 修复stat_poly_eq的报错并标注斜率p值 #

你遇到的报错是因为：之前用的x是group_age（离散因子），而stat_poly_eq()依赖线性模型lm()，要求x是连续变量。另外，你原来的mutate(time = as.integer(age))是错误的——age是字符串（比如"d3"），转成integer会返回NA，这也导致拟合失败！

解决步骤： #

1. 先修正time变量的生成逻辑，把"d3"/"d8"/"d24"正确映射为数值3/8/24；
2. 在stat_poly_eq()的aes里指定x=time（连续变量）、group=group，并自定义标注内容（比如只显示斜率的p值）；
3. 用label.x和label.y控制标注位置，确保每个group上方都能显示。

首先安装并加载ggpmisc：

install.packages("ggpmisc")
library(ggpmisc)

然后添加stat_poly_eq()的代码：

stat_poly_eq(aes(x=time, group=group, label = paste0("p = ", stat(p.value.label))),
             formula = y ~ x, parse = TRUE,
             label.x = "left", label.y = 10,  # 根据你的数据范围调整位置
             size = 4)

完整修正后的代码 #

library(dplyr)
library(plotly)
library(ggplot2)
library(ggpmisc)

set.seed(1)
df <- data.frame(value = c(rnorm(500,8,1),rnorm(600,6,1.5),rnorm(400,4,0.5),rnorm(500,2,2),rnorm(400,4,1),rnorm(600,7,0.5),rnorm(500,3,1),rnorm(500,3,1),rnorm(500,3,1)), 
                 age = c(rep("d3",500),rep("d8",600),rep("d24",400),rep("d3",500),rep("d8",400),rep("d24",600),rep("d3",500),rep("d8",500),rep("d24",500)), 
                 group = c(rep("A",1500),rep("B",1500),rep("C",1500))) %>% 
  # 修正time变量的生成逻辑，正确转为连续数值
  dplyr::mutate(time = case_when(
    age == "d3" ~ 3,
    age == "d8" ~ 8,
    age == "d24" ~ 24
  )) %>% 
  dplyr::arrange(group,time) %>% 
  dplyr::mutate(group_age=paste0(group,"_",age))
df$group_age <- factor(df$group_age,levels=unique(df$group_age))

# 最终绘图代码
ggplot(df,aes(x=group,y=value,fill=age,color=age)) + 
  geom_violin(position = position_dodge(width = 0.8), alpha=0.5) + 
  geom_boxplot(width=0.1, aes(middle=mean(value)), 
               position = position_dodge(width = 0.8), alpha=0.5) + 
  geom_smooth(aes(x=time, group=group), color="black", method='lm', size=1, se=T) + 
  stat_poly_eq(aes(x=time, group=group, label = paste0("p = ", stat(p.value.label))),
               formula = y ~ x, parse = TRUE,
               label.x = "left", label.y = 10,
               size = 4) +
  guides(alpha = "none") +
  theme_minimal()

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