r - 如何使用 R 中的网格图形绘制直方图？

Question

我想制作一个类似于beanplot或violin绘图的绘图，但不是对称形状，而是我想在中心线的不同侧绘制两个不同的分布。本文 [pdf] http://www.jstatsoft.org/v28/c01/paper中的图 4 类似。

我想grid在R. 任何指针都会很有用。我查看了lattice包和包中的histogram功能lattice，但这不是我想要做的。

任何帮助/指针表示赞赏。

Answer 1

panel.violin在 Sarkar 的 lattice 包中，通过对函数进行相当简单的修改，您可以相当容易地获得半小提琴图。该函数中有四行可以在grid.polygon调用中更改，以从“双面”密度绘图更改为单面密度绘图。首先require(lattice)和require(grid)。然后键入panel.violin。我将向您展示“水平”更改：

require(grid)
panel.violin2  <-  
    # snipped all the arguments and processing

       grid.polygon(x = c(dx.list[[i]] ), 
   # Notice I removed: ... ,rev(dx.list[[i]]) 
               y = c(2*dy.list[[i]]    ), default.units = "native", 
   # Removed: ...  , -rev(dy.list[[i]])
               name = trellis.grobname(identifier, type = "panel", 
                 group = group), gp = gpar(fill = col, col = border, 
                 lty = lty, lwd = lwd, alpha = alpha))

也将其从 else{ ...} 子句的相应部分中删除。现在您可以使用 help(panel.violin) 中的示例运行它

bwplot(voice.part ~ height, singer,
       panel = function(..., box.ratio) {
           panel.violin2(..., col = "transparent",
                        varwidth = FALSE, box.ratio = box.ratio)
           panel.bwplot(..., fill = NULL, box.ratio = .1)
       } )

如果您想要另一边的密度，您需要做的就是删除 dx.list[[i]] 和 dy.list[[i]] 并留在rev(dx.list[[i]])and中-rev(dy.list[[i]])。