R — заполнение пустого списка операциями в другом списке

Я пишу программу R для анализа древовидной структуры. В моем примере ниже в дереве 10 узлов, и предки каждого узла (родители этого узла и родители родителей этого узла и т. д.) хранятся в списке под названием Ancestors. Пользователь будет запрашивать вектор имен узлов, и я пытаюсь создать список, который будет заполнен предками этого запроса. Каждый элемент в списке будет содержать список потомков запроса для каждого из вызываемых предков. Пожалуйста, смотрите ниже для примера

Допустим, у меня есть следующая структура.

Таким образом, список предков будет выглядеть так.

Ancestors <- list()
Ancestors$'p1' <- c('p2', 'p3', 'p4', 'p5', 'p8', 'p9', 'p10')
Ancestors$'p2' <- c('p4', 'p5', 'p8', 'p9', 'p10')
Ancestors$'p3' <- c('p4', 'p5', 'p9', 'p10')
Ancestors$'p4' <- c('p5', 'p9', 'p10')
Ancestors$'p5' <- c('p9', 'p10')
Ancestors$'p6' <- c('p4', 'p5', 'p9', 'p10')
Ancestors$'p7' <- c('p5', 'p9', 'p10')
Ancestors$'p8' <- c('p5', 'p9', 'p10')
Ancestors$'p9' <- NA
Ancestors$'p10' <- NA

Допустим, запрос

query <- c('p5', 'p4', 'p1')

Тогда список, который я хотел бы создать,

# lst <- list()
# 
# lst$'p2'
#   'p1'
# lst$'p3'
#   'p1'
# lst$'p4'
#   'p1'
# lst$'p5'
#   'p1', 'p4'
# lst$'p8'
#   'p1'
# lst$'p9'
#   'p1', 'p4', 'p5'
# lst$'p10'
#   'p1', 'p4', 'p5'

(2,3,4,5,8,9,10) — это все предки, существующие для терминов запроса. Вот список, который я хотел бы составить. Затем для каждого из этих именованных элементов я хотел бы составить список терминов запроса, которые являются потомками элемента списка. Прошу прощения за запутанный пример. Я надеюсь, что это имеет смысл.

Вот что я пробовал до сих пор

lst <- list()

lapply(query, function(x) {
  theAncestors <- Ancestors[[x]]
  sapply(theAncestors, function(y) {
    lst[[y]][[1]] <- c(lst[[y]][[1]], x)
  })
})

Но это не заполняет список lst. Все, что происходит, это то, что он распечатывает

[[1]]
  p9  p10 
"p5" "p5" 

[[2]]
  p5   p9  p10 
"p4" "p4" "p4" 

[[3]]
  p2   p3   p4   p5   p8   p9  p10 
"p1" "p1" "p1" "p1" "p1" "p1" "p1" 

что немного отличается от того, что я хочу. Кроме того, когда я пытаюсь вывести lst, он все еще пуст. Так что этот код даже не влияет на lst. Итак, как я могу получить желаемый результат? Я думал об использовании цикла for, но я думаю, что он очень медленный в R. Моя реальная проблема, вероятно, будет состоять из сотен или тысяч терминов запроса и многих других терминов-предков. Так что это будет очень долго. Поэтому я думаю, что цикл for, вероятно, не сработает.

Редактировать: я понял это. Мой код сейчас:

lst <- list()

aLst <- unlist(lapply(query, function(x) {
  theAncestors <- Ancestors[[x]]
  sapply(theAncestors, function(y) {
    lst[[y]][1] <- c(lst[[y]][[1]], x)
  })
}))
aLst <- split(unname(aLst), names(aLst))

Это распечатывает

$p10
[1] "p5" "p4" "p1"

$p2
[1] "p1"

$p3
[1] "p1"

$p4
[1] "p1"

$p5
[1] "p4" "p1"

$p8
[1] "p1"

$p9
[1] "p5" "p4" "p1"

что я и хотел