Как сохранить значения списка в матрице

Question

set.seed(650)
library(maxLik)
y = c(rnorm(15,1,1), rnorm(15,3,1))
dat = data.frame(y)
B = 3 # number bootstrap sample
n = length(dat$y)
n1 = 15
boot.samples = matrix(sample(dat$y, size = B * n, replace = TRUE), n, B)
ml = list()
boot.l = 0
va.l = NULL
for (j in 1:B) {
  boot.l = boot.samples[, j]
  for (i in 1:n) {
    LLl <- function(param) {
      mul <- param[1]
      sigmal <- param[2]
      sum(log(dnorm(dat[1:i, ], mul, sigmal)))
    }
    ml[[i]] = coef(maxLik(logLik = LLl, start = c(mul = 1, sigmal = 1)))
  }
  va.l = matrix(unlist(ml), n-1, B*2, byrow = TRUE)
}
va.l

Ниже приведены мои выходные данные

Однако, когда я печатаю список, у меня появляется следующий вывод.

У меня вопрос: как я могу получить муль-оценки для j = 1 в 1-м столбце, сигма-оценки для j = 1 во втором столбце и муль-оценки для j = 2?в 3-м столбце сигма-оценки для j = 2 в 4-м столбце и т. д.

Есть ли другой способ сделать это?Спасибо за вашу помощь.