0

我正在尝试从目录中读取文件，然后将它们转换为列表，然后尝试将它们转换为二进制矩阵以进行打乱图。到目前为止，我已经完成了。

步骤1

    setwd("/run/media/punit/data2/PRIMARY_CELL_ATAC_SEQ/RNA_seq/literature/subtype_list/")
    
    
    files <- list.files(pattern = ".*tsv")
    files
    dflist <- list()
    for(i in 1:length(files)){
      dflist[[i]] <- read_tsv(files[i])
    }

my_names = c("Gene")
result_abd = lapply(dflist, "[",  my_names)

result_abd的输出

尝试了第一个 list_to_matrix

它没有用

list_to_matrix(result_abd)
Error in `[<-`(`*tmp*`, unique(lt[[i]]), i, value = 1) : 
  invalid subscript type 'list'

然后尝试了这个答案 它有效

out <- t(splitstackshape:::charMat(listOfValues = result_abd, fill = 0L))

使用上面的对象我运行了 uptR

upset(out,nsets = 8,nintersects = NA,matrix.color = "blue", 
      number.angles = 30, point.size = 3.5, line.size = 2,
      sets.bar.color = "#56B4E9",mb.ratio = c(0.55, 0.45),group.by = "degree",sets.x.label = "Gene Size",
      mainbar.y.label =  "DisGeNET",text.scale = 2,
      #boxplot.summary = c("minAgeRange"),
      order.by = "freq", empty.intersections = "off",keep.order = TRUE,
      queries = list(list(query = intersects, 
                          params = list("Dis"), 
                          color = "orange")))

错误

r in names(elem_data) <- c(names, "freq") : 
  'names' attribute [9] must be the same length as the vector [3]

任何建议或帮助将不胜感激

4

0 回答 0