r - 复制、重复还是其他？

Question

我有如下详述的代码。我需要的基本描述是重复此代码 30 次，并尽可能将最终值分配为向量。每次运行模拟时，我还需要将种子增加一。所以我在考虑使用复制或重复。但是作为复制需要rep(vector, times)，我不确定，因为我没有要分配的特定向量。

为了更好地解释代码，A 是值的第一个最佳猜测，然后通知 B。B 然后用作 C 的输入，通知 D。C 和 D 然后迭代直到收敛三个值。

rep(vector,40){
set.seed(100)

...............
A
..............
B
..............

for (n in 1:50){
...............
C
..............
D
..............
}
}

Answer 1

您可以将所有内容包装在一个函数中，然后使用 sapply 将索引提供给您的函数。

yourfun <- function(seed){
    set.seed(seed)
    # your code here

    ans <- rnorm(1)
    # just so we can return something
    return(ans)
}

sapply(seq(500, 530), yourfun)