我需要把这个代码中的值放在一个20x20矩阵中，以便能够在以后的计算中引用每一行

Question

我需要能够在以后的计算中使用这20个向量，所以我想创建一个20x20矩阵，我可以将代码中的每一行存储到一行中，或者如果你可以建议一种更好的方法来引用代码。

lambda=12
mu=2
answer = numeric(20)
k = 0
for (i in 1:20) {
Pi_i <- numeric(i)
m <- 1:i

k = k+((1/factorial(i)*(lambda/mu)^i))
for (j in m) {

Pi_i[j] = (k^-1)*((lambda/mu)^j)/factorial(j)
}
cat(Pi_i,"nn")
}

Answer 1

您可以创建一个空的20 x 20矩阵，并将每个向量写入i行。矩阵中不能有空单元格，因此必须将正确数量的NA值附加到每个向量的末尾：

lambda=12
mu=2
answer = numeric(20)
k = 0
result <- matrix(0, ncol = 20, nrow = 20)
for (i in 1:20) {
Pi_i <- numeric(i)
m <- 1:i

k = k+((1/factorial(i)*(lambda/mu)^i))
for (j in m) {

Pi_i[j] = (k^-1)*((lambda/mu)^j)/factorial(j)
}
result[i,] <- c(Pi_i, rep(NA, 20 - i))

}
result