如何在R Markdown中用循环生成含ggplotly的多标签页？

想要创建包含ggplotly交互式图表的R Markdown文档，用for/lapply/purrr::map循环生成标签页，但所有图表都集中显示在最后一个标签页里。如果去掉ggplotly只用ggplot，代码能正常运行；单独分块生成单个图表时，图表也能正确显示在对应标签页。试过Stack Overflow相关方案和plot.new()，但plot.new()会把最后一个标签页的图表也删掉。测试代码如下：

## **Count Distribution** {.tabset}
```{r, results='asis', warning=FALSE}
# List of dataframes for the first tabset
library(ggplot2)
library(reshape2)
library(plotly)
hist_func <- function(matrix) {
  hist_plot <- ggplot(melt(as.data.frame(matrix))) +
                        geom_line(stat = 'density', aes(x = value, ..scaled.., colour = variable)) +
                        labs(x = NULL) + theme(legend.position = 'right') + scale_x_log10()

  return(hist_plot)
}

# Set seed 
set.seed(42)

# Generate random gene count matrix
c1 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c2 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c3 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c4 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c5 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)

# List of dataframes for the  tabset
dataframes_tab1 <- list(counts = c1, rawlog2 = c2, TMM = c3, DESEQ2 = c4, VOOM = c5)

plt <- htmltools::tagList()

# Generate tabset content using loop
plotlist=list()
for (i in seq_along(dataframes_tab1)) {
  
  df_name <- names(dataframes_tab1)[i]
  plotlist[[i]] <- ggplotly(hist_func(dataframes_tab1[[df_name]]))
  cat(paste('### ', df_name, '\n \n'))
  cat(' \n')

  # trying multiple approaches
  plt[[i]] <- as.widget(ggplotly(plotlist[[i]]))
  
  # to check if either method works #as taglist or without
  plotlist[[i]]
  plt[[i]]
  cat(' \n\n')
}

# 解决方法

问题核心是**Markdown的{.tabset}语法与htmlwidgets（ggplotly属于此类）的输出机制不兼容**：用`results='asis'`结合`cat`生成Markdown标签时，ggplotly的HTML输出无法被自动关联到对应标签的内容区域，最终所有图表被追加到最后一个标签里。

推荐两种可行方案：

## 方案1：用htmltools直接构建标签页结构
放弃Markdown的{.tabset}语法，直接用`htmltools`包的`tabsetPanel`和`tabPanel`构建标签页，每个图表会被明确绑定到对应标签，彻底避免错位问题。

修正后的代码：
```r
library(ggplot2)
library(reshape2)
library(plotly)
library(htmltools)

hist_func <- function(matrix) {
  hist_plot <- ggplot(melt(as.data.frame(matrix))) +
    geom_line(stat = 'density', aes(x = value, ..scaled.., colour = variable)) +
    labs(x = NULL) + theme(legend.position = 'right') + scale_x_log10()
  return(hist_plot)
}

# 生成模拟数据
set.seed(42)
c1 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c2 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c3 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c4 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c5 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
dataframes_tab1 <- list(counts = c1, rawlog2 = c2, TMM = c3, DESEQ2 = c4, VOOM = c5)

# 构建每个标签页内容
tab_list <- lapply(names(dataframes_tab1), function(name) {
  p <- ggplotly(hist_func(dataframes_tab1[[name]]))
  tabPanel(title = name, p)
})

# 渲染标签页集合
tabsetPanel(
  tab_list,
  type = "tabs" # 指定标签样式，可选
)

方案2：用knitr::asis_output统一输出流 #

如果坚持使用Markdown的{.tabset}语法，需要用knitr::asis_output包裹所有输出内容，保证Markdown标题和ggplotly图表的输出顺序与关联性：

修正后的代码：

## **Count Distribution** {.tabset}
```{r, results='asis', warning=FALSE}
library(ggplot2)
library(reshape2)
library(plotly)
library(knitr)

hist_func <- function(matrix) {
  hist_plot <- ggplot(melt(as.data.frame(matrix))) +
    geom_line(stat = 'density', aes(x = value, ..scaled.., colour = variable)) +
    labs(x = NULL) + theme(legend.position = 'right') + scale_x_log10()
  return(hist_plot)
}

set.seed(42)
c1 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c2 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c3 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c4 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
c5 <- matrix(sample(0:10000, 10 * 1000, replace = TRUE), nrow = 10)
dataframes_tab1 <- list(counts = c1, rawlog2 = c2, TMM = c3, DESEQ2 = c4, VOOM = c5)

for (i in seq_along(dataframes_tab1)) {
  df_name <- names(dataframes_tab1)[i]
  # 输出标签标题
  asis_output(paste('### ', df_name, '\n'))
  # 生成并输出对应图表
  p <- ggplotly(hist_func(dataframes_tab1[[df_name]]))
  print(p)
  # 添加换行分隔
  asis_output('\n\n')
}

### 关键原因说明
- ggplotly生成的是htmlwidgets对象，这类对象的输出需要被knitr正确解析并插入到HTML的对应位置；
- 直接用`cat`生成Markdown标签时，ggplotly的输出会脱离Markdown解析流，导致所有图表被集中到最后一个标签的内容区域；
- 用htmltools构建标签页时，每个图表被明确绑定到对应的`tabPanel`；用`asis_output`则保证了Markdown标题和图表输出的顺序性与关联性，两种方式都能解决错位问题。

内容的提问来源于stack exchange，提问作者Aryan Neupane