如何在散点图上叠加箱线图？

Question

假数据

set.seed( 123)
x<-rnorm(1000, mean=60,sd=20)
y <-  exp(-10 +  .95*log(x^3)) + rnorm(1000,mean=1,sd=1)
df <- data.frame(x,y)
cls.x <- quantile(df$x, seq(.1, .9, by=.1))
df$x.class <- findInterval(df$x, cls.x)
df$x.class <- as.factor(df$x.class)
head(df)

以下都不行

plot(df$x,df$y,col=3)
par(new=T)
boxplot(y~x.class, data=df,xlab="",ylab="",xaxt="n")

也不是这个

boxplot(y~x.class, data=df,xlab="",ylab="",xaxt="n")
points(df$x,df$y,col=3)

使用 ggplot，我得到的最接近的是使用类似的东西

library(ggplot2)
ggplot(df,aes(x.class,y))+geom_boxplot() + geom_point()

不幸的是，它没有显示轴的真实可变性。 我尝试使用 jitter 选项，但无法强制绘图使用 X 变量的真实可变性

非常感谢任何建议。

Ps：我知道 Rlab 中的 bplot.xy() 函数，但是，该函数不允许我更改箱线图的颜色，或者先绘制点。

library(Rlab)
bplot.xy( x,y,  N=10)
points( x,y, pch=".", col=3, cex=3)

Answer 1

以下是你想要的吗？

library(ggplot2)

ggplot(df, aes(x, y, fill = x.class)) + 
  geom_point(alpha = 0.10) +
  geom_boxplot(alpha = 0.50)  

Answer 2

您的df$x 从 3 到 124 不等，而您的 x 轴从 1 到 10。在基本图形中，您可以这样做：

plot(jitter(as.integer(df$x.class)), df$y, col=3, type="p", xlab = "", ylab = "", xaxt = "n")
boxplot(y~x.class, data=df,xlab="",ylab="",xaxt="n", add = TRUE)

我添加了jitter 以帮助打破分发。您也可以尝试使用pch=16 使点变成实心，或者使用透明度（例如，col="#aaffaa22" 用于点）。