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

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()

不幸的是，它没有显示出轴上的实际变化。 我尝试使用抖动选项，但无法强制绘图使用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"作为点）。