误差棒点阵图

Question

我想用误差线和箱线图来表示一个树因子点图。我已经完成了几乎所有的图形，除了我的错误栏是重叠的。具体来说，我的想法是抖动点和误差线

 set.seed(1)

data<-sample(rnorm(50,2,2),40)
sd<-sample(rnorm(50,0.2,2),40)
factor<-sample( LETTERS[1:3], 40, replace=TRUE, prob=c(0.3,.3, 0.3) )
frame<-data.frame(data,sd,factor)
frame<-as.data.frame(frame)
up=frame$data+frame$sd
lo=frame$data-frame$sd
bwplot(data ~ factor, frame,
       ylim=c((min(lo)+0.5*min(lo)),(max(up)-0.5*min(lo))),
       up=frame$data+frame$sd,
       lo=frame$data-frame$sd,
       panel=function(x,y,...){
         panel.bwplot(x,y,...)
         panel.stripplot(x,y,pch=16 ,alpha=0.7, jitter=TRUE, factor=0.2,
                         col="grey",cex=1.2,...)
         panel.arrows(x0=x, y0=lo,
                      x1=x, y1=up,code=3,
                      angle=90, length=0.05)

       }
)

谢谢！

Answer 1

单独抖动

set.seed(1)

data<-sample(rnorm(50,2,2),40)
sd<-sample(rnorm(50,0.2,2),40)
factor<-sample( LETTERS[1:3], 40, replace=TRUE, prob=c(0.3,.3, 0.3) )
frame<-data.frame(data,sd,factor)
frame<-as.data.frame(frame)
up=frame$data+frame$sd
lo=frame$data-frame$sd
bwplot(data ~ factor, frame,
       ylim=c((min(lo)+0.5*min(lo)),(max(up)-0.5*min(lo))),
       up=frame$data+frame$sd,
       lo=frame$data-frame$sd,
       panel=function(x,y,...){
         panel.bwplot(x,y,...)
xj=jitter(as.numeric(x))
         panel.stripplot(xj,y,pch=16 ,alpha=0.7, factor=0.2,
                         col="grey",cex=1.2,...)
         panel.arrows(x0=xj, y0=lo,
                      x1=xj, y1=up,code=3,
                      angle=90, length=0.05)

       }
)

然后条状图和误差条（箭头）的抖动量相同。

误差棒点阵图

Error bar lattice dotplot

r

panel

lattice