主页 / coding / 问题 / 78942585
Accepted
DaniCee
Asked: 2024-09-03 11:33:57 +0800 CST

为 cowplot 中的多个图创建一个公共 x 轴

  • 772

我使用 制作了多个图，并将它们ggplot2排列在一起cowplot。它们的长度都相同，并且 x 轴在相同的位置有相同的断点，但标签不同。

这是一个 MWE：

mylist <- list()
mylist[["ID1"]][["sequence"]] <- "AGAATATTATACATTCATCT"
mylist[["ID2"]][["sequence"]] <- "GCTAGCGTTTAGTTTAGCTG"
mylist[["ID3"]][["sequence"]] <- "AACCCTTTAAACTCGAAGTA"
set.seed(123)
mylist[["ID1"]][["data"]] <- data.frame(time=1:100, value=rnorm(100, mean=10, sd=2))
mylist[["ID2"]][["data"]] <- data.frame(time=1:100, value=rnorm(100, mean=2, sd=1))
mylist[["ID3"]][["data"]] <- data.frame(time=1:100, value=rnorm(100, mean=5, sd=3))
indexes <- seq(5, 100, length.out=20)
all_df <- NULL
for (id in names(mylist)){
  seqsplit <- unlist(strsplit(mylist[[id]][["sequence"]], ""))
  ind_df <- data.frame(call=seqsplit, time=indexes)
  final_df <- dplyr::left_join(mylist[[id]][["data"]], ind_df, by="time")
  xcolors <- ifelse(seqsplit=="A", "green", ifelse(seqsplit=="C", "blue", ifelse(seqsplit=="G", "black", "red")))
  P <- ggplot2::ggplot(final_df, ggplot2::aes(x=time, y=value)) +
    ggplot2::geom_line(linewidth=0.5) +
    ggplot2::scale_x_continuous(breaks=indexes, labels=seqsplit, expand=c(0,1)) +
    ggplot2::scale_y_continuous(breaks=seq(-5, 15, 5), limits=c(-5,15)) +
    ggplot2::theme_light() +
    ggplot2::theme(axis.title=ggplot2::element_blank(),
                   axis.text.x=ggtext::element_markdown(face="bold", color=xcolors))
  mylist[[id]][["plot"]] <- P
}
plot_list <- sapply(mylist, "[", "plot")
grDevices::pdf(file="test.pdf", height=4, width=10)
print(
  cowplot::plot_grid(plotlist=plot_list, ncol=1)
)
grDevices::dev.off()

生成结果：

图1

我现在想要完成的是为所有人制作一个共同的 x 轴，具有共同点indexes（或者更确切地说从 1 到序列的长度），我会将其放在cowplot；类似于：

图2

这可能涉及将 x 轴索引标尺制作成ggplot2相同宽度的单独标尺（考虑 y 轴水平空间和所有），然后将其放置在cowplot具有适当的相对高度上......欢迎任何帮助！

1 个回答

  1. Best Answer

    一种方法是在第一个图的顶部添加第二个 x 轴：

    for (id in names(mylist)){
...
  P <- ggplot(final_df, ggplot2::aes(x=time, y=value)) +
    geom_line(linewidth=0.5) +
    scale_y_continuous(breaks=seq(-5, 15, 5), limits=c(-5,15)) +
    theme_light() 
    

  if(id=="ID1") {
    P <- P + scale_x_continuous(breaks=indexes, labels=seqsplit, expand=c(0,1),
                                     sec.axis = sec_axis(~.,
                                                         breaks=indexes, labels=seq_along(indexes))) +
      theme(axis.title=element_blank(),
            axis.text.x.top = ggtext::element_markdown(color='white', fill="steelblue", 
                                                       padding=unit(c(1, 8, 1, 8), "pt")))
  }
  else {
    P <- P + scale_x_continuous(breaks=indexes, labels=seqsplit, expand=c(0,1)) +
      theme(axis.title=element_blank(),
            axis.text.x=ggtext::element_markdown(face="bold", color=xcolors))
  }

  mylist[[id]][["plot"]] <- P
}

    在此处输入图片描述

    • 2

相关问题