好久以前就觉得这种分布图非常好,好像origin可以实现;曾经试图用SAS or R实现它,也没有思路;
今天无意中看到jitter函数,成功实现,呵呵
效果图
code
## box-dotplot 散点和箱式图的合体
no.sample <- 100
y <- factor(c(rep("ctrl",no.sample),rep("treat",no.sample)))
x <- c(rnorm(no.sample,mean=0,sd=1),rnorm(no.sample,mean=2,sd=1))
boxplot(x ~ y, border=4, col="light grey", boxwex=0.5)
points(jitter(rep(1:2, each=no.sample), 1),unlist(split(x, y)),col=2,cex=0.5, pch=16)
## 两组样本不均衡的时候,需要更改jitter