如何将列表值存储在矩阵中

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进行多个估计，在第二列中对j = 1进行一个估计，在第三列中对j = 2进行mul估计，在第二栏中对j = 2进行估计第四列，依此类推？

还有其他方法吗？谢谢您的帮助。