如何基于多列对应关系使用S4Vectors::expand展开GRanges对象行

要展开GRanges对象并保持txid与txname的对应关系，可通过以下两种方式实现：

方法1：基础Bioconductor操作 #

# 计算每行需要重复的次数（等于对应txid列表的元素数量）
repeat_n <- lengths(mcols(gr)$txid)

# 重复GRanges的基因组区域行
gr_expanded <- gr[rep(seq_len(nrow(gr)), repeat_n)]

# 展开元数据列，将列表转为匹配长度的向量
mcols(gr_expanded)$txid <- unlist(mcols(gr)$txid)
mcols(gr_expanded)$txname <- unlist(mcols(gr)$txname)

方法2：tidyverse工具链 #

若习惯使用tidyverse语法，可借助tibble和unnest完成操作：

library(tidyverse)
library(GenomicRanges)

gr_expanded <- as_tibble(gr) %>%
  unnest(c(txid, txname)) %>%
  as_granges()

关于S4Vectors::expand()的问题 #

expand()的核心逻辑是生成元数据列的笛卡尔积，因此无法满足行内元素对应展开的需求。上述两种方法都是先根据元数据列表的长度重复对应的基因组行，再将列表列拆分为与行数匹配的向量，从而确保txid和txname的对应关系始终一致。

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