Объедините соответствующие элементы списков в вектор или список

Я собираюсь использовать %dopar% и foreach, и мне нужно объединить выходные данные.

Функция, которая будет вызываться параллельно, имеет на выходе список постоянной длины для каждого вызова. Однако длины элементов этого списка не всегда постоянны.

После объединения я хочу, чтобы результат был максимально упрощен, но при этом позволял мне идентифицировать список (итерацию) происхождения каждого элемента.

Длина B/c различается на этом более глубоком уровне, подобные ответы меня не совсем понимают.

Вот некоторые примеры данных:

list1 <- list(rnorm(1), rnorm(1), rnorm(1), rnorm(8))
list2 <- list(rnorm(1), rnorm(1), rnorm(1), rnorm(8))
list3 <- list(rnorm(1), rnorm(1), rnorm(1), rnorm(14))

do.call(Map, c(c, list(list1, list2, list3))) дает:

    [[1]]
    [1] -0.2923462  0.4891224 -0.5080176

    [[2]]
    [1]  0.3229466  0.9511572 -0.9815504

    [[3]]
    [1] -1.160413  0.707568 -1.564874

    [[4]]
     [1] -1.13093146  0.06791923  0.65380844  1.01829862  0.47360903  0.68616334 -1.07166155 -1.54018814 -0.60860430  1.64524185  0.40222817 -0.54747627
    [13] -1.73420011  0.67861611  0.55527953  1.36454409  0.40215155 -0.65706184 -0.71008434 -1.11484886 -0.69811408 -0.45451101 -0.85574891 -0.79241329
    [25]  0.31018144 -0.03212242 -1.55192430 -2.19142725 -1.85528112  0.85204097

и do.call(Map, c(list, list(list1, list2, list3))) дает:

[[1]]
[[1]][[1]]
[1] -0.2923462

[[1]][[2]]
[1] 0.4891224

[[1]][[3]]
[1] -0.5080176


[[2]]
[[2]][[1]]
[1] 0.3229466

[[2]][[2]]
[1] 0.9511572

[[2]][[3]]
[1] -0.9815504


[[3]]
[[3]][[1]]
[1] -1.160413

[[3]][[2]]
[1] 0.707568

[[3]][[3]]
[1] -1.564874


[[4]]
[[4]][[1]]
[1] -1.13093146  0.06791923  0.65380844  1.01829862  0.47360903  0.68616334 -1.07166155 -1.54018814

[[4]][[2]]
[1] -0.6086043  1.6452418  0.4022282 -0.5474763 -1.7342001  0.6786161  0.5552795  1.3645441

[[4]][[3]]
 [1]  0.40215155 -0.65706184 -0.71008434 -1.11484886 -0.69811408 -0.45451101 -0.85574891 -0.79241329  0.31018144 -0.03212242 -1.55192430 -2.19142725
[13] -1.85528112  0.85204097

EDIT, правильный ответ должен выглядеть примерно так (простите за ГСЧ):

part1 <- do.call(Map, c(c, list(list1, list2, list3)))
part2 <- do.call(Map, c(list, list(list1, list2, list3)))
correct <- list(part1[[1]], part1[[2]], part1[[3]], part2[[4]])
correct
[[1]]
[1]  1.80341685 -0.06408827  0.07004951

[[2]]
[1]  0.4389224 -0.1821140  0.2538133

[[3]]
[1]  0.008303713 -1.004631075  1.936738072

[[4]]
[[4]][[1]]
[1] -0.86790931  1.20414809  0.04373068 -0.49097606  1.12826503 -0.76263091 -0.93364770  0.13392904

[[4]][[2]]
[1] -1.0823008 -0.4382813  1.4328709 -0.8961412  0.8350054  1.4855032 -1.3800748  1.4300227

[[4]][[3]]
 [1]  0.02126034  0.30640618  0.49420442  0.72107997  0.97666620 -0.48049810  1.22227279 -1.00918452 -0.23290645 -1.27834163  2.55142878  1.07120297
[13]  1.37473759  0.72308135

Я также должен отметить, что элементы списка не обязательно должны быть числовыми — это могут быть выходные данные модели, например, из jags().

Использование c дает правильную первую часть (относительно того, что я хочу), а использование list дает правильную последнюю часть. Как получить лучшее из обоих миров?