R 使用命名列表重命名列表列表列表

因此，我正在处理一个包含其他列表的列表，其结构如下：

library(graph)
library(RBGL)
library(Rgraphviz)

show(tree)

$`SO:0001968`
$`SO:0001968`$`SO:0001622`
$`SO:0001968`$`SO:0001622`$`SO:0001624`
$`SO:0001968`$`SO:0001622`$`SO:0001624`$`SO:0002090`
[1] 1

$`SO:0001968`$`SO:0001622`$`SO:0001623`
$`SO:0001968`$`SO:0001622`$`SO:0001623`$`SO:0002091`
[1] 1

$`SO:0001968`$`SO:0001969`
$`SO:0001968`$`SO:0001969`$`SO:0002090`
[1] 1

$`SO:0001968`$`SO:0001969`$`SO:0002091`
[1] 1


dput(tree)
list(`SO:0001968` = list(`SO:0001622` = list(`SO:0001624` = list(
    `SO:0002090` = 1), `SO:0001623` = list(`SO:0002091` = 1)), 
    `SO:0001969` = list(`SO:0002090` = 1, `SO:0002091` = 1)))

我用于构建列表的数据来自一个名为

g

的对象：

show(g)

A graphNEL graph with directed edges
Number of Nodes = 7 
Number of Edges = 8 


dput(g)
new("graphNEL",
nodes = c("SO:0001968", "SO:0001969", "SO:0001622", 
"SO:0001623", "SO:0001624", "SO:0002090", "SO:0002091"), edgeL = list(
    `SO:0001968` = list(edges = 3:2), `SO:0001969` = list(edges = 6:7), 
    `SO:0001622` = list(edges = 5:4), `SO:0001623` = list(edges = 7L), 
    `SO:0001624` = list(edges = 6L), `SO:0002090` = list(edges = integer(0)), 
    `SO:0002091` = list(edges = integer(0))), edgeData = new("attrData",

    data = list(`SO:0001968|SO:0001622` = list(weight = 1), `SO:0001968|SO:0001969` = list(
        weight = 1), `SO:0001969|SO:0002090` = list(weight = 1), 
        `SO:0001969|SO:0002091` = list(weight = 1), `SO:0001622|SO:0001624` = list(
            weight = 1), `SO:0001622|SO:0001623` = list(weight = 1), 
        `SO:0001623|SO:0002091` = list(weight = 1), `SO:0001624|SO:0002090` = list(
            weight = 1)), defaults = list(weight = 1)), nodeData = new("attrData",

    data = list(`SO:0001968` = list(label = "coding_transcript_variant"), 
        `SO:0001969` = list(label = "coding_transcript_intron_variant"), 
        `SO:0001622` = list(label = "UTR_variant"), `SO:0001623` = list(
            label = "5_prime_UTR_variant"), `SO:0001624` = list(
            label = "3_prime_UTR_variant"), `SO:0002090` = list(
            label = "3_prime_UTR_intron_variant"), `SO:0002091` = list(
            label = "5_prime_UTR_intron_variant")), defaults = list(
        label = NA_character_)), renderInfo = new("renderInfo",

    nodes = list(), edges = list(), graph = list(), pars = list()), 
    graphData = list(edgemode = "directed"))

每个

SO:000XXX

对应一个名称，我可以使用返回命名列表的函数

nodeData

查找名称：

nodeData(g, nodes(g), "label")

$`SO:0001968`
[1] "coding_transcript_variant"

$`SO:0001969`
[1] "coding_transcript_intron_variant"

$`SO:0001622`
[1] "UTR_variant"

$`SO:0001623`
[1] "5_prime_UTR_variant"

$`SO:0001624`
[1] "3_prime_UTR_variant"

$`SO:0002090`
[1] "3_prime_UTR_intron_variant"

$`SO:0002091`
[1] "5_prime_UTR_intron_variant"

我需要的是用

nodeData

函数的相应字符串替换（或重命名）树

列表中的数据

例如，将
nodeData
函数中的
编码抄本变量的
树
列表中的
'SO:0001968'
替换为
nodeData（）
函数中的
。
如果将
nodeData（）
的输出保存到向量，则可以使用
名称（）
函数将名称分配给
列表（）

为列表元素指定名称的示例：

x <- 1:5
y <- 11:20
z <- 21:25

theList <- list(x,y,z)

listNames <- c("element1","element2","element3")
names(theList) <- listNames
# access first element by name, using $ form of extract operator
theList$element1

您可能需要
unlist（）
节点数据（）的输出，如下所示：

theNames <- unlist(nodeData(g, nodes(g), "label"))
names(g) <- theNames 

theNames此递归函数应实现以下功能：

# you will do this but I couldn't install your packages
# nodeD <- nodeData(g, nodes(g), "label")

nodeD <- list(`SO:0001968` = "coding_transcript_variant",
              `SO:0001969` = "coding_transcript_intron_variant",
              `SO:0001622` = "UTR_variant",
              `SO:0001623` = "5_prime_UTR_variant",
              `SO:0001624` = "3_prime_UTR_variant",
              `SO:0002090` = "3_prime_UTR_intron_variant",
              `SO:0002091` = "5_prime_UTR_intron_variant")

rename_items <- function(item){
  if (is.list(item)){
    item <- lapply(item,rename_items)
    names(item) <- unname(nodeD[names(item)])
  }
  item
}

tree2 <- rename_items(tree)

您使用的是什么软件包，
data.tree
？另外，将您的
树
和
g
对象，
dput（树）
。请提供可复制的数据。@zx8754刚刚添加了它。无法复制
g
对象，错误：error in.local（.object，…）：未使用的参数（数据=列表（
SO:0001968 | SO:0001622
=列表（权重=1）…@zx8754是的，我刚检查过，你是对的，似乎是dput（）的内部错误，因为我尝试了另一个graphnel，但也不起作用。我需要名为
树
的特定列表中的名称。而且，如果我这样做，向量甚至没有相同的长度，因此它不会起作用。@JFernandez-如果无法访问您的数据来复制您的问题，我所能做的最好是提供一个我自己生成的数据示例.
# you will do this but I couldn't install your packages
# nodeD <- nodeData(g, nodes(g), "label")

nodeD <- list(`SO:0001968` = "coding_transcript_variant",
              `SO:0001969` = "coding_transcript_intron_variant",
              `SO:0001622` = "UTR_variant",
              `SO:0001623` = "5_prime_UTR_variant",
              `SO:0001624` = "3_prime_UTR_variant",
              `SO:0002090` = "3_prime_UTR_intron_variant",
              `SO:0002091` = "5_prime_UTR_intron_variant")

rename_items <- function(item){
  if (is.list(item)){
    item <- lapply(item,rename_items)
    names(item) <- unname(nodeD[names(item)])
  }
  item
}

tree2 <- rename_items(tree)

# $coding_transcript_variant
# $coding_transcript_variant$UTR_variant
# $coding_transcript_variant$UTR_variant$`3_prime_UTR_variant`
# $coding_transcript_variant$UTR_variant$`3_prime_UTR_variant`$`3_prime_UTR_intron_variant`
# [1] 1
# 
# 
# $coding_transcript_variant$UTR_variant$`5_prime_UTR_variant`
# $coding_transcript_variant$UTR_variant$`5_prime_UTR_variant`$`5_prime_UTR_intron_variant`
# [1] 1
# 
# 
# 
# $coding_transcript_variant$coding_transcript_intron_variant
# $coding_transcript_variant$coding_transcript_intron_variant$`3_prime_UTR_intron_variant`
# [1] 1
# 
# $coding_transcript_variant$coding_transcript_intron_variant$`5_prime_UTR_intron_variant`
# [1] 1