ggplot2技术求助：替换geom_point为geom_boxplot及绘图设置问题

1. 替换geom_point()为geom_boxplot()的报错解决 #

原代码报错是因为x=ID为数值型变量，geom_boxplot()默认无法直接对连续型x轴生成箱线图（需明确分组依据）。推荐将ID转换为因子类型，让ggplot识别为离散分组，彻底解决分组问题。

2. 移除x轴20和40刻度 #

数据中ID仅包含10、30、50三个值，只需通过刻度设置指定显示的刻度即可：若ID转为因子型，使用scale_x_discrete()会自动只显示存在的因子水平，自然移除20、40这类不存在的刻度。

3. 子图标题加入数学表达式 #

借助facet_wrap()的labeller参数，搭配label_bquote()可以在子图标题中插入数学表达式（如$N=48$）。label_bquote()支持用.引用变量值，同时兼容R的plotmath数学语法。

完整优化代码 #

n=48
p=7
df=data.frame(BIAS = c(rnorm(n*p,2,.01),rnorm((n+10-2),4),rnorm(n*p,1.5,.01),rnorm((n+30-2),1),rnorm(n*p,1,.01),rnorm((n+50-2),0)),
              ID   = rep(c(10,30,50),c((n*p+(n+10-2)),(n*p+(n+30-2)),(n*p+(n+50-2)))),
              TYPE = rep(c("DELTA","PSI","DELTA","PSI","DELTA","PSI"),c(n*p,(n+10-2),n*p,(n+30-2),n*p,(n+50-2))))

# 将ID转换为因子，解决箱线图分组问题
df$ID <- factor(df$ID)

ggplot(df,aes(x=ID ,y=BIAS, color=TYPE)) + 
  geom_boxplot() +  # 替换为箱线图
  facet_wrap(~TYPE, scales = "free", 
             labeller = label_bquote(.(TYPE) ~ N == .(n))) +  # 子图标题加入数学表达式
  theme_bw() + 
  xlab("T") + 
  ggtitle("N = 48") +
  scale_x_discrete()  # 自动移除不存在的20、40刻度

额外说明 #

• 若想使用填充色区分组别，可将color=TYPE改为fill=TYPE，同时添加color="black"设置箱线边框，让图形更清晰
• 数学表达式支持希腊字母、上标下标等复杂格式，只需遵循R的plotmath语法即可

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