R编程模拟

Question

我一直在研究我在 R 中的模拟问题。问题是

一名矿工被困在矿井中，他有 3 扇门可供选择。
门 1 = 他可以从矿井中出来，需要 2 小时，
Door2 = 他将在 3 小时内旅行，但只会返回矿井，
门 3 = 与门 2 相同的情况，但 5 小时。

我已经创建了一个随机数生成器，如果它得到 1 就会停止（使用 while 循环），但我的问题是如何重复此操作 100 次并将结果放入表格列中。做这类问题有没有更简单的解决方案？ （这里是 R 初学者）

这是我的代码：

library(purrr) 
storage <-c() 
d = NULL 
x1 <- rdunif(1, 3, a=1) 

while ( x1 != 1) { 
  x1<- rdunif(1, 3, a=1) 
  storage <-c(storage, x1) 
} 

storage

Answer 1

• 将 100 个实验存储在一个包含 100 个部分的列表中
• n表示1号门的编号，如果出现1号门，则进入下一个实验。
• trial 表示实验顺序

---

storage <- vector("list", 100)
n <- 0
trial <- 1
while(n < 100){
    x <- sample(1:3, 1)
    storage[[trial]] <- c(storage[[trial]], x)
    if(x == 1){
        n <- n + 1
        trial <- trial + 1
    }
}

storage

之后，我建议您不要将列表转换为数据框，因为每个试验的长度不同。它将使用数据框生成许多 NA 值。

library(dplyr)
hour <- lapply(storage, recode, `1` = 2, `2` = 3, `3` = 5)
total <- lapply(hour, sum) %>% unlist

lapply是处理列表数据的好工具。