r - 我需要多次运行公平骰子模拟

Question

我有一个模拟掷骰子 100 次的程序。我需要知道如何运行这个程序 10^5 次，我认为这与数字有关。

 set.seed(123)

 x <- sample(1:6, size=100, replace = TRUE)

 hist(x,
 main="10^6 fair rolls",
 xlab = "Dice Result",
 ylab = "Probability",
 xlim=c(0.5,6.5),
 breaks=-1:100+.5,
 prob=TRUE )

Answer 1

正如@markus 所建议的，您可以使用replicate：

set.seed(123)

nTime <- 10^5

x <- replicate(nTime, sample(1:6, size=100, replace = TRUE))

hist(x,
     main="10^6 fair rolls",
     xlab = "Dice Result",
     ylab = "Probability",
     xlim=c(0.5,6.5),
     breaks=-1:100+.5,
     prob=TRUE )