Question

努力将凌乱不等长的data.frame从宽表转换为长表，然后折叠（汇总）新变量。目前看起来像这样，Gene是一个变量，GO_terms是一个包含多个逗号分隔值的变量：

Gene            GO_terms                        
AA1006G00001    GO:0098655, GO:0008643, GO:0005351, GO:0005886, GO:0016021      
AA100G00001     GO:0098655, GO:0009944, GO:0009862, GO:0010075, GO:0010014, GO:0009855, GO:0010310
AA100G00002     GO:0098655, GO:0008643, GO:0005886

我想做的第一步是将其转换为“长”格式，因此看起来像这样：

Gene            GO_terms 
AA1006G00001    GO:0098655
AA1006G00001    GO:0008643
AA1006G00001    GO:0005351
AA1006G00001    GO:0005886
AA1006G00001    GO:0016021
AA100G00001     GO:0001666
AA100G00001     GO:0009944
AA100G00001     GO:0009862
AA100G00001     GO:0010075
AA100G00001     GO:0010014
AA100G00001     GO:0009855
AA100G00001     GO:0010310
AA100G00002     GO:0008270
AA100G00002     GO:0005634
AA100G00002     GO:0005886
AA100G00003     GO:0005488
AA100G00003     GO:0005634

然后，我希望通过切换两个变量来重组data.table，因为它的整理如下：

GO_terms    Genes
GO:0005351  AA1006G00001        
GO:0005886  AA1006G00001,   AA100G00002 
GO:0008643  AA1006G00001,   AA100G00002 
GO:0009855  AA100G00001     
GO:0009862  AA100G00001     
GO:0009944  AA100G00001     
GO:0010014  AA100G00001     
GO:0010075  AA100G00001     
GO:0010310  AA100G00001     
GO:0016021  AA1006G00001        
GO:0098655  AA1006G00001,   AA100G00001,      AA100G00002

包含变量的基因可以在一列内（以逗号分隔），也可以在多列中。

请问有人能够提供tidyr，reshape2或dplyr解决方案吗？

编辑：dput()表是：

structure(list(`Gene    ` = c("AA1006G00001\t", "AA100G00001\t", 
"AA100G00002\t"), `GO_terms                     ` = c("GO:0098655, GO:0008643, GO:0005351, GO:0005886, GO:0016021\t\t", 
"GO:0098655, GO:0009944, GO:0009862, GO:0010075, GO:0010014, GO:0009855, GO:0010310", 
"GO:0098655, GO:0008643, GO:0005886")), row.names = c(NA, -3L
), class = c("tbl_df", "tbl", "data.frame"), spec = structure(list(
    cols = list(`Gene   ` = structure(list(), class = c("collector_character", 
    "collector")), `GO_terms                        ` = structure(list(), class = c("collector_character", 
    "collector"))), default = structure(list(), class = c("collector_guess", 
    "collector"))), class = "col_spec"))

Answer 1

似乎您正在执行一些GO分析。您可以尝试使用inverseList中的topGO（Bioconductor中用于GO分析的最受欢迎的R软件包之一）：

library(topGO)

gene.to.go <- strsplit(gsub('\t', '', df$GO_terms), ', ', fixed = TRUE)
names(gene.to.go) <- gsub('\t', '', df$Gene)

go.to.gene <- inverseList(gene.to.go)

data.frame(GO_term = names(go.to.gene), Genes = sapply(go.to.gene, paste0, collapse = ', '),
           stringsAsFactors = FALSE, row.names = NULL)

#       GO_term                                  Genes
# 1  GO:0005351                           AA1006G00001
# 2  GO:0005886              AA1006G00001, AA100G00002
# 3  GO:0008643              AA1006G00001, AA100G00002
# 4  GO:0009855                            AA100G00001
# 5  GO:0009862                            AA100G00001
# 6  GO:0009944                            AA100G00001
# 7  GO:0010014                            AA100G00001
# 8  GO:0010075                            AA100G00001
# 9  GO:0010310                            AA100G00001
# 10 GO:0016021                           AA1006G00001
# 11 GO:0098655 AA1006G00001, AA100G00001, AA100G00002

实际上，如果使用readMappings中的topGO导入GO映射文件，则对数据进行操作会更容易。

Answer 2

这是一个解决问题的方法：

library(tidyr)
library(dplyr)

#allow up to seven Genes per GO_term if there is more increase the letters expression
long<-df %>% separate(GO_terms, into=paste0("a", 1:100), sep=", ", extra="merge") %>% 
  gather( key="key", value="GO_terms", -Gene)

#filter data frame, remove the NA and keep the desired columns
long<-long[!is.na(long$GO_terms), c("Gene", "GO_terms")]

final<-long %>% group_by(GO_terms) %>% summarize( Gene=toString(Gene) )

重塑然后折叠数据帧

2 个答案: