这是做你想做的吗?它需要一些数据处理和两把小提琴的绘制。
set.seed(1)
dat <- data.frame(x=1, y=rnorm(10^5))
#calculate for each point if it's central or not
dat_q <- quantile(dat$y, probs=c(0.025,0.975))
dat$central <- dat$y>dat_q[1] & dat$y < dat_q[2]
#plot; one'95' violin and one 'all'-violin with transparent fill.
p1 <- ggplot(data=dat, aes(x=x,y=y)) +
geom_violin(data=dat[dat$central,], color="transparent",fill="gray90")+
geom_violin(color="black",fill="transparent")+
theme_classic()
编辑:圆边困扰着我,所以这里是第二种方法。如果我这样做,我会想要直线。所以我做了一些演奏与密度(这是小提琴地块是基于)
d_y <- density(dat$y)
right_side <- data.frame(x=d_y$y, y=d_y$x) #note flip of x and y, prevents coord_flip later
right_side$central <- right_side$y > dat_q[1]&right_side$y < dat_q[2]
#add the 'left side', this entails reversing the order of the data for
#path and polygon
#and making x negative
left_side <- right_side[nrow(right_side):1,]
left_side$x <- 0 - left_side$x
density_dat <- rbind(right_side,left_side)
p2 <- ggplot(density_dat, aes(x=x,y=y)) +
geom_polygon(data=density_dat[density_dat$central,],fill="red")+
geom_path()
p2
啊,刚才打我! – Axeman
@Axeman伟大的思想相似吗?我已经添加了第二种方法。 – Heroka
@贺冈,这太棒了!按照你的猜测,我自己偶然发现了你的第一个方法,但并不满意。你的第二种方法正是我想要的。非常感谢! – dewarrn1