可用于R中Tufte boxplots的功能？

本文介绍了可用于R中Tufte boxplots的功能？的处理方法，对大家解决问题具有一定的参考价值，需要的朋友们下面随着小编来一起学习吧！ 问题描述 我有一些数据，我已经分成了足够多的分组，标准箱图显得非常拥挤。 Tufte有他自己的盒子，其中你基本放弃了全部或部分盒子，如下所示：一些示例数据： code $ c $变换（ChickWeight， Time = cut（ChickWeight $ Time，4）） cw $小鸡< - as.factor（sample（LETTERS [ sew（3）]，nrow（cw），replace = TRUE）） levels（cw $ Diet） 我想为每个Diet * Time * Chick分组设置一个箱子重量。 几年前，我遇到了这个问题，并使用网格图形将解决方案拼凑在一起，我将稍后发布。但是在解决这个新的（和类似的）问题时，我想知道是否有一种股票方式来做这些事情，而不是修复我一起出现的问题。 似乎是Tufte的创作中不太受欢迎的人，但我真的很喜欢他们在大量的分组中密集地展示分布模式，如果他们在ggplot2或格子中具有良好的功能，我会更多地使用它们。 / p> 解决方案您显然只想要一个垂直版本，所以我使用了panel.bwplot代码，剥离了所有非必需品，例如盒子和帽子，并在参数中设置horizo​​ntal = FALSE并创建了一个panel.tuftebxp函数。同时将点的cex设置为默认值的一半。还有相当多的选项可以根据你的口味进行调整。 时间的数字因素名称看起来马虎，但我认为概念验证很清楚，您可以清理对您来说重要的事情： pch = box.dot $ pch，col = box.dot $ col， alpha = box.dot $ alpha，cex = box.dot $ cex，font = box.dot $ font， fontfamily = box.dot $ fontfamily，fontface = box.dot $ fontface，$ b $ fill = box.rectangle $ fill，varwidth = FALSE，notch = FALSE， notch.frac = 0.5，.. 。，levels.fos = if（horizo​​ntal）sort（unique（y））else sort（unique（x））， stats = boxplot.stats，coef = 1.5，do.out = TRUE，identifier =bwplot ） { if（all（is.na（x）| is.na（y））） return（）x< - as.numeric（x） y< - as.numeric（y） box.dot< - trellis.par.get（box.dot） box.rectangle< - trellis.par.get （box.rectangle） box.umbrella< - trellis.par.get（box.umbrella） plot.symbol< - trellis.par.get（plot.symbol） fontsize.points< - trellis.par.get（ fontsize）$ points cur.limits< - current.panel.limits（） xscale< - cur.limits $ xlim yscale< - cur.limits $ ylim if（！notch） notch.frac ＃移除横码 blist blist。 < - lapply（blist，[[，out） blist.height< - box.width if（varwidth）{ maxn< - max （x）） blist.n< - sapply（blist，[[，n） blist.height< - sqrt（blist.n / maxn）* blist.height } blist.conf< - if（notch） sapply（blist，[[，conf） else t（blist.stats [，c（2 ，4），drop = F ALSE]） ybnd< - cbind（blist.stats [，3]，blist.conf [2，]，blist.stats [， 4]，blist.stats [，4]，blist .conf [2，]，blist.stats [， 3]，blist.conf [1，]，blist.stats [，2]，blist.stats [， 2]，blist.conf [1，]，blist.stats [，3]） xleft< - levels.fos - blist.height / 2 xright< - levels.fos + blist.height / 2 xbnd< - cbind（xleft + notch.frac * blist.height / 2，xleft， xleft，xright，xright，xright - notch.frac * blist.height / 2， xright，xright ，xleft，xleft，xleft + notch.frac * blist.height / 2） xs ys blist.stats [，4]），rep（levels.fos，2），c（ blist.stats [， 1]，blist.stats [，5]），col = box.umbrella $ col，alpha = box.umbrella $ alpha， lwd = bo x.umbrella $ lwd，lty = box.umbrella $ lty，identifier = paste（标识符，whisker，sep =。）） if（all（pch == |））{ mult 1 - notch.frac else 1 panel.segments（levels.fos - mult * blist.height / 2， blist.stats [，3]，levels.fos + mult * blist.height / 2， blist.stats [，3]，lwd = box.rectangle $ lwd，lty = box .rectangle $ lty， col = box.rectangle $ col，alpha = alpha，identifier = paste（标识符，dot，sep =。））} else { panel.points（x = levels.fos，y = blist.stats [，3]， pch = pch，col = col，alpha = alpha，cex = cex，标识符=粘贴（标识符，点，sep =。））} panel.points（x = rep（levels.fos，sapply（blist.out，length ））， y = unlist（blist.out），pch = plot.symbol $ pch，col = plot.symbol $ col， alpha = plot.symbol $ alpha，cex = plot.symbol $ cex * 0.5，标识符=粘贴（标识符，离群值，sep =。）） } bwplot（重量〜饮食+时间+小鸡，数据= cw，面板= 函数（x，y，...）panel.tuftebxp（x = x，y = y，...）） I have some data that I've divided into enough groupings that standard boxplots look very crowded. Tufte has his own boxplots in which you basically drop all or half of box, like this:Some sample data:cw <- transform(ChickWeight, Time = cut(ChickWeight$Time,4) )cw$Chick <- as.factor( sample(LETTERS[seq(3)], nrow(cw), replace=TRUE) )levels(cw$Diet) <- c("Low Fat","Hi Fat","Low Prot.","Hi Prot.")I want a boxplot of weight for every Diet * Time * Chick grouping.I had this problem come up years ago, and kludged together a solution using grid graphics, which I'll post in a bit. But in solving this new (and similar) problem I'm wondering if there's a stock way to do them rather than fixing my kludged together example.As an aside, these seem to be amongst the less-beloved of Tufte's creations, but I really like them for densely displaying patterns of distributions across a large number groupings, and I'd use them more if there was a good function for them in ggplot2 or lattice. 解决方案 You apparently wanted just a vertical version, so I took the panel.bwplot code, stripped out all the non-essentials such as the box and the cap, and set horizontal=FALSE in the arguments and created a panel.tuftebxp function. Also set the cex of the points at half of the default. There are still quite a few of options left that could be adjusted to your tastes. The "numeric" factor names for "Time" look sloppy but I figure the "proof of concept" is clear and you can clean up what is important to you:panel.tuftebxp <- function (x, y, box.ratio = 1, box.width = box.ratio/(1 + box.ratio), horizontal=FALSE, pch = box.dot$pch, col = box.dot$col, alpha = box.dot$alpha, cex = box.dot$cex, font = box.dot$font, fontfamily = box.dot$fontfamily, fontface = box.dot$fontface, fill = box.rectangle$fill, varwidth = FALSE, notch = FALSE, notch.frac = 0.5, ..., levels.fos = if (horizontal) sort(unique(y)) else sort(unique(x)), stats = boxplot.stats, coef = 1.5, do.out = TRUE, identifier = "bwplot") { if (all(is.na(x) | is.na(y))) return() x <- as.numeric(x) y <- as.numeric(y) box.dot <- trellis.par.get("box.dot") box.rectangle <- trellis.par.get("box.rectangle") box.umbrella <- trellis.par.get("box.umbrella") plot.symbol <- trellis.par.get("plot.symbol") fontsize.points <- trellis.par.get("fontsize")$points cur.limits <- current.panel.limits() xscale <- cur.limits$xlim yscale <- cur.limits$ylim if (!notch) notch.frac <- 0 #removed horizontal code blist <- tapply(y, factor(x, levels = levels.fos), stats, coef = coef, do.out = do.out) blist.stats <- t(sapply(blist, "[[", "stats")) blist.out <- lapply(blist, "[[", "out") blist.height <- box.width if (varwidth) { maxn <- max(table(x)) blist.n <- sapply(blist, "[[", "n") blist.height <- sqrt(blist.n/maxn) * blist.height } blist.conf <- if (notch) sapply(blist, "[[", "conf") else t(blist.stats[, c(2, 4), drop = FALSE]) ybnd <- cbind(blist.stats[, 3], blist.conf[2, ], blist.stats[, 4], blist.stats[, 4], blist.conf[2, ], blist.stats[, 3], blist.conf[1, ], blist.stats[, 2], blist.stats[, 2], blist.conf[1, ], blist.stats[, 3]) xleft <- levels.fos - blist.height/2 xright <- levels.fos + blist.height/2 xbnd <- cbind(xleft + notch.frac * blist.height/2, xleft, xleft, xright, xright, xright - notch.frac * blist.height/2, xright, xright, xleft, xleft, xleft + notch.frac * blist.height/2) xs <- cbind(xbnd, NA_real_) ys <- cbind(ybnd, NA_real_) panel.segments(rep(levels.fos, 2), c(blist.stats[, 2], blist.stats[, 4]), rep(levels.fos, 2), c(blist.stats[, 1], blist.stats[, 5]), col = box.umbrella$col, alpha = box.umbrella$alpha, lwd = box.umbrella$lwd, lty = box.umbrella$lty, identifier = paste(identifier, "whisker", sep = ".")) if (all(pch == "|")) { mult <- if (notch) 1 - notch.frac else 1 panel.segments(levels.fos - mult * blist.height/2, blist.stats[, 3], levels.fos + mult * blist.height/2, blist.stats[, 3], lwd = box.rectangle$lwd, lty = box.rectangle$lty, col = box.rectangle$col, alpha = alpha, identifier = paste(identifier, "dot", sep = ".")) } else { panel.points(x = levels.fos, y = blist.stats[, 3], pch = pch, col = col, alpha = alpha, cex = cex, identifier = paste(identifier, "dot", sep = ".")) } panel.points(x = rep(levels.fos, sapply(blist.out, length)), y = unlist(blist.out), pch = plot.symbol$pch, col = plot.symbol$col, alpha = plot.symbol$alpha, cex = plot.symbol$cex*0.5, identifier = paste(identifier, "outlier", sep = "."))}bwplot(weight ~ Diet + Time + Chick, data=cw, panel= function(x,y, ...) panel.tuftebxp(x=x,y=y,...)) 这篇关于可用于R中Tufte boxplots的功能？的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！