r - 如何在箱线图中绘制与水平线平行的条形图到箱线图的面积分数？

Question

我正在处理数据表示，我需要并行绘制带有堆叠条形图的箱线图。

用不同深浅的颜色按点覆盖面积的比例着色。我该如何绘制它？

Answer 1

这是一个通用解决方案，可用于将平行百分比框添加到格子弓形图中。它计算箱线图描述部分（晶须和主要部分）的确切百分比。

该解决方案基于网格包。lattice 的视口操作比 ggplot2 更容易，这就是我选择 bwplot 的原因。

我生成了一些示例数据：

df <- data.frame(cond = factor( rep(c("A"), each=200) ), 
                 rating = c(rnorm(200),rnorm(200, mean=.8)))

library(lattice)
bwplot(cond ~ rating, data=df,gp=gpar(fill='blue'),
       par.settings = list( box.umbrella=list(col= c( "red")), 
                            box.dot=list(col= c("green")), 
                            box.rectangle = list(fill= c( "blue"),alpha=0.6)))

我抓住了主视口

downViewport("plot_01.panel.1.1.vp")

然后我存储胡须的尺寸和四肢箱线图位置的位置。

segs <- grid.get('segments',grep=T)
cap.segs <- grid.get('cap.segments',grep=T)

函数 drawBox ，用百分比的文本绘制一个矩形。我为每个盒子调用了 3 次。

drawbox <- function(x0 = segs$x0[1] , col ='red',
                    width.vp = segs$x0[2] - segs$x0[1]){
  vpd <- viewport(
            x = x0 ,
            y = cap.segs$y1[2] + unit(0.5,'native'),
            width  = width.vp,
            height = unit(2,'cm'),
            just   = c('left','bottom'),
            name   = 'vpd')
  pushViewport(vpd)
    grid.rect(gp=(gpar(fill=col)))
    # The compute of percent is a little bit tricky due we can't apply '/'
    value <- as.numeric(convertUnit(width.vp,'native'))
    value <- value/as.numeric(convertUnit(cap.segs$x0[2]- cap.segs$x0[1],'native'))
    grid.text(label = paste(round(value*100),'%',sep='') , gp= gpar(cex=3))
  upViewport(1)
}

函数drawbox被调用了3次：

drawbox()
drawbox(col='yellow',width=segs$x0[1] - cap.segs$x0[1], x0= cap.segs$x0[1])
drawbox(col='green',width.vp= cap.segs$x0[2]- segs$x0[2],x0 = segs$x0[2])

Answer 2

一种可能性是结合boxplot()和rect()功能。坐标rect()被计算为您感兴趣的分位数。在这个例子中，我使用了分位数 25% 和 60%（25% 和 35% 的总和）。text()用于设置名称的函数。

  set.seed(123)
  x<-rnorm(100)
  boxplot(x,horizontal=TRUE,axes=FALSE)
  rect(min(x),1.5,quantile(x,0.25),1.4,col="red")
  rect(quantile(x,0.25),1.5,quantile(x,0.60),1.4,col="green")
  rect(quantile(x,0.60),1.4,max(x),1.5,col="yellow")
  text(-1.5,1.45,"25%")
  text(0,1.45,"35%")
  text(1.1,1.45,"40%")