R语言如何镜像两个组图合成

mkshixfv 于 2023-01-22 发布在其他

关注(0)|答案(1)|浏览(145)

我已经有了下面的情节，我已经很满意了

library(ggplot2)
library(wesanderson)
library(ggsignif)
library(ggdist)

 
dat1 <- as.data.frame(c(175.070, 279.760, 310.120, 454.370, 356.190, 398.720, 394.790, 427.940, 450.090, 452.240, 344.630, 374.450, 401.810,
          366.250, 455.150, 455.150, 284.820, 325.640, 437.840, 455.150, 455.150, 454.900, 433.970, 346.830, 455.050, 455.150,
          448.740, 455.150, 455.000, 455.150, 243.920, 295.660, 282.430, 305.190, 368.580, 348.440, 437.100, 453.180, 454.100,
          44.464, 420.280, 455.150, 427.250, 101.490, 454.100, 454.810, 323.690, 454.410, 235.160, 294.870, 455.150, 330.460,
          289.610, 236.310, 454.810, 299.020,  29.289,  39.988, 348.340, 455.150, 439.540, 409.910, 440.310, 433.090, 455.150,
          198.890, 454.100, 346.370, 261.780, 280.760, 427.180, 455.150, 233.740, 334.040, 452.240, 411.890, 378.310, 354.410,
          347.670, 439.970))
colnames(dat1) = "value"

dat1$group <- rep(NA,nrow(dat1))
for (i in 1:nrow(dat1)){
  if (i <= 40){
    dat1[i,]$group <- print('C')
  }else{
    dat1[i,]$group <- print('I')
  }
}

dat1$group <- factor(dat1$group, levels = c('C','I'))

pltt <- rev(wes_palette(n=2, name="Darjeeling1"))

#Model 1

#M1 Loglikelihood
p1 <- ggplot(dat1, aes(x=group, y=value, color=group, fill=group)) +
  scale_fill_manual(values=pltt, guide = "none") +
  scale_color_manual(values=pltt, guide = "none") +
  stat_halfeye(adjust = .33, width = .6, alpha=0.3, position = position_nudge(x=0.3)) +
  geom_jitter(shape=16, position=position_jitter(0.1), alpha=0.5) +
  geom_signif(comparisons = list(c("C", "I")), tip_length = -0.01, y_position = -0.5, 
              map_signif_level=TRUE, test="wilcox.test", color = "black") +
  geom_boxplot(position = position_nudge(x=0.2), width=0.05,outlier.shape = NA, fill=alpha(0.2)) +
  ylim(0,460) +
  theme(legend.position = "none",
        axis.text=element_text(size=10),
        axis.title=element_text(size=12),
        plot.title = element_text(hjust = 0.5,size=14,face="bold")) +
  theme_light() +
  ggtitle("Data distribution")

现在我想知道是否有一种方法，可以颠倒顺序，或者把这个组合镜像成这样。

我读过这个post，我看到有一种方法可以用密度图来做，但我想知道是否有人知道一种方法可以用像我这样的复合图来做。
先谢了！

r

来源：https://stackoverflow.com/questions/75171150/how-to-mirror-two-group-plot-compositions

1条答案

按热度按时间

cotxawn71#

一个选项是将每个geom与组中的一个subset一起使用。要获得右侧的stat_halfeye，可以将参数side用于两者。下面是一个可重现的示例：

library(ggplot2)
library(wesanderson)
library(ggsignif)
library(ggdist)

ggplot() +
  stat_halfeye(data = subset(dat1, group == "C"),
               aes(x=group, y=value, color=group, fill=group), 
               adjust = .33, width = .6, alpha=0.3, position = position_nudge(x=-0.3), side = "left") +
  geom_jitter(data = subset(dat1, group == "C"),
              aes(x=group, y=value, color=group, fill=group),
              shape=16, position=position_jitter(0.1), alpha=0.5) + 
  geom_boxplot(data = subset(dat1, group == "C"),
               aes(x=group, y=value, color=group, fill=group), 
               position = position_nudge(x=-0.2), width=0.05,outlier.shape = NA, fill=alpha(0.2)) +
  stat_halfeye(data = subset(dat1, group == "I"),
               aes(x=group, y=value, color=group, fill=group), 
               adjust = .33, width = .6, alpha=0.3, position = position_nudge(x=0.3), side = "right") +
  geom_jitter(data = subset(dat1, group == "I"),
              aes(x=group, y=value, color=group, fill=group),
              shape=16, position=position_jitter(0.1), alpha=0.5) +
  geom_boxplot(data = subset(dat1, group == "I"),
               aes(x=group, y=value, color=group, fill=group), 
               position = position_nudge(x=0.2), width=0.05,outlier.shape = NA, fill=alpha(0.2)) +
  geom_signif(data = dat1,
              aes(x=group, y=value, color=group, fill=group),
              comparisons = list(c("C", "I")), tip_length = -0.01, y_position = -0.5, 
              map_signif_level=TRUE, test="wilcox.test", color = "black") +
  ylim(0,460) +
  theme(legend.position = "none",
        axis.text=element_text(size=10),
        axis.title=element_text(size=12),
        plot.title = element_text(hjust = 0.5,size=14,face="bold")) +
  theme_light() +
  ggtitle("Data distribution")

创建于2023年1月19日，使用reprex v2.0.2

赞(0）回复(0）举报 2023-01-22

我来回答

R语言如何镜像两个组图合成

1条答案

相关问题

热门标签

最新问答

R语言 如何镜像两个组图合成

1条答案

相关问题

热门标签

最新问答

R语言如何镜像两个组图合成