R中排列的并排图

Question

我试图在R中并排放置两个图，并有以下示例。

library(vioplot)
x <- rnorm(100)
y <- rpois(100,1)
plot(x, y, xlim=c(-5,5), ylim=c(-5,5),type='n')
vioplot(x, col="tomato", horizontal=TRUE, at=-4, add=TRUE,lty=2, rectCol="gray")
vioplot(y, col="cyan", horizontal=TRUE, at=-3, add=TRUE,lty=2)
vioplot(y, col="cyan", horizontal=TRUE, at=-2, add=TRUE,lty=2)

有了这些数据，我就可以对我的x和y变量进行vioplot了。 现在，例如，我想开发与左侧每个vioplot相关的单独计数数据的条形图。

counts <- c(10, 20, 30)
barplot(counts, main="Car Distribution", horiz=TRUE)

我使用了mtcars示例，但它可以是任何计数数据。 我想知道是否可以并排生成这些图，以便计数图与vioplot正确对齐。 计数图不需要任何y-axis标签。

Answer 1

根据您的规格ggplot是我的建议

library(tidyverse)
p1 <- lst(x, y, y1=y) %>% 
  bind_cols() %>% 
  pivot_longer(1:3) %>% 
  ggplot(aes(name, value)) + 
   geom_violin(trim = FALSE)+
   geom_boxplot(width=0.15) + 
   coord_flip()
p2 <- mtcars %>% 
  count(gear) %>% 
  ggplot(aes(gear, n)) + 
   geom_col()+
   coord_flip()
cowplot::plot_grid(p1, p2)

在基础R中你可以做（​​请注意，我使用了箱线图，但也应该与 viopülot 一起使用）

par(mfrow=c(1,2))
counts <- table(mtcars$gear)
boxplot(cbind(x,y,y), col="tomato", horizontal=TRUE,lty=2, rectCol="gray")
barplot(counts, main="Car Distribution", horiz=TRUE,
        names.arg=c("3 Gears", "4 Gears", "5 Gears"))

Answer 2

如果你想使用 ggplot，另一个选择是来自 ggpubr 的函数 ggarrange()。

library(dplyr)
library(ggplot2)
library(ggpubr)

# Create a sample dataset
dt <- tibble(group = rep(c("x", "y"), each = 100)) %>% 
        mutate(value = if_else(group == "x", rnorm(200),
                               as.double(rpois(200, 1))))

# Combined violin/Box plot
violins <- dt %>% 
        ggplot(aes(value, group)) +
        geom_violin(width = 0.5) +
        geom_boxplot(width = 0.1)

# Bar chart
bars <- dt %>% 
        ggplot(aes(group)) +
        geom_bar(width = 0.1) +
        coord_flip()

# Combine
ggpubr::ggarrange(violins, bars + rremove("ylab") + rremove("y.text"), ncol = 2)

输出：

Answer 3

您可以使用以下代码：

library(vioplot)
x <- rnorm(100)
y <- rpois(100,1)
par(mfrow=c(1,2))
plot(x, y, xlim=c(-5,5), ylim=c(-5,-1),type='n')
vioplot(x, col="tomato", horizontal=TRUE, at=-4, add=TRUE,lty=2, rectCol="gray")
vioplot(y, col="cyan", horizontal=TRUE, at=-3, add=TRUE,lty=2)
vioplot(y, col="cyan", horizontal=TRUE, at=-2, add=TRUE,lty=2)

counts <- table(mtcars$gear)
barplot(counts, main="Car Distribution", horiz=TRUE,
        names.arg=c("3 Gears", "4 Gears", "5 Gears"))

输出：

Answer 4

感谢您提出有趣的问题，这促使我探索基本 R 图形功能。 我试图找到一个案例，小提琴图和条形图之间的并排配置提供了有意义的关系。 情况是我有一个iris数据的子集，其中包含不同的物种数量。 我想显示三个统计数据：

1. 通过显示条形图来采样每个物种的计数；
2. 通过显示小提琴图，每个采样物种中萼片长度的分布； 和
3. 每个采样物种的中值花瓣宽度，通过定位小提琴图。

我在这里按照@GW5 的想法创建可以控制轴上位置的条形图。 我在这里遵循@IRTFM 的想法来调整轴的原点。

这是完整的代码：

library(vioplot)

some_iris <- iris[c(1:90, 110:139), ]
ir_counts <- some_iris |> with(Species) |> table()
ir_counts
#    setosa versicolor  virginica 
#        50         40         30 
ir_names <- names(ir_counts)
ir_colors <- c("cyan", "green", "pink")

x_vio1 <- some_iris |> subset(Species == ir_names[1]) |> with(Sepal.Length)
x_vio2 <- some_iris |> subset(Species == ir_names[2]) |> with(Sepal.Length)
x_vio3 <- some_iris |> subset(Species == ir_names[3]) |> with(Sepal.Length)

y_vio1 <- some_iris  |> subset(Species == ir_names[1]) |> with(Petal.Length) |> median() 
y_vio2 <- some_iris  |> subset(Species == ir_names[2]) |> with(Petal.Length) |> median() 
y_vio3 <- some_iris  |> subset(Species == ir_names[3]) |> with(Petal.Length) |> median() 

# `xpd = FALSE` to keep the grid inside the plotting boxes.

par(mfrow = c(1, 2), xpd = FALSE)

# The violin plots, put on the left side.

plot(NULL, 
     xlim = c(0, 10), ylim = c(0, 10), type = "n", las = 1, xaxs =  "i", yaxs = "i", 
     xlab = "Sepal Length (cm)", ylab = " Median Petal Width (cm)")

vioplot(x_vio1, col = ir_colors[1], horizontal = TRUE, at = y_vio1, add = TRUE, lty = 2)
vioplot(x_vio2, col = ir_colors[2], horizontal = TRUE, at = y_vio2, add = TRUE, lty = 2)
vioplot(x_vio3, col = ir_colors[3], horizontal = TRUE, at = y_vio3, add = TRUE, lty = 2)
grid()

# The texts that informs the names of the species

text(labels = ir_names, y = c(y_vio1, y_vio2, y_vio3), 
    x = c (min(x_vio1), min(x_vio2), min(x_vio3)) - 1)

# The barplots, put on the right side.

plot(NULL,
     xlim = c(0, 60), ylim = c(0, 10), yaxt = "n", type = "n", 
     las = 1, xlab = "Counts", ylab = "", xaxs = "i", yaxs = "i"
)
rect(xleft = 0, xright = ir_counts[1],
     ybottom = y_vio1 - 0.3, ytop = y_vio1 + 0.3, col = ir_colors[1])
rect(xleft = 0, xright = ir_counts[2],
     ybottom = y_vio2 - 0.3, ytop = y_vio2 + 0.3, col = ir_colors[2])
rect(xleft = 0, xright = ir_counts[3],
     ybottom = y_vio3 - 0.3, ytop = y_vio3 + 0.3, col = ir_colors[3])
grid()

结果如下：

如果您想在条形图（右侧）上放置标签，您可以使用mtext ，如下所示：

# ... (The same code above)
mtext(text = ir_names, side = 2, at = c(y_vio1, y_vio2, y_vio3), 
      line = 0.2, las  = 1 )

结果标签：