xci <- do.call(rbind, tapply( d.pizza$delivery_min, d.pizza$driver,
MeanCI, conf.level=0.99, na.rm=TRUE))
PlotDotCI(xci, main="delivery_min ~ driver",
args.legend=list(y=-1.5, legend=c("estimate", "99%-CI")))
Run the code above in your browser using DataLab