如何做R循環

Question

我有一個關於R循環的問題。

例如，目前在t = 0時，有100人活着。基本上，每個人都會以μ= 0.1的指數（-mu）的概率存活。

我想要生成10個樣本來獲得t = 1時存活的人數。所以我已經完成並獲得以下內容。

命令：

set.seed(123) 
alive <- 100  
mu <- 0.1 
sample <- 10 
alive1 <- rbinom(sample,alive,exp(-mu)) 
alive1 

# [1] 92 88 91 87 86 95 90 87 90 91 

，現在，我想保持繼續這樣做直到時間t = 20。

命令：

alive2 <- rbinom(10,alive1,exp(-mu)) 
alive2 
alive3 <- rbinom(10,alive2,exp(-mu)) 
alive3 
.... 

alive20 <-rbinom (10,alive19,exp(-mu)) 
alive20 

輸出：

alive2 <- rbinom(10,alive1,exp(-mu)) 
alive2 

# [1] 78 80 81 78 81 82 83 83 83 77 

alive3 <- rbinom(10,alive2,exp(-mu)) 
alive3 

# [1] 67 71 72 63 72 73 75 75 77 72 

...

然而

，我不想繼續重複命令，特別是如果我想我的時間延長至更長的時間。我如何在r中爲我的問題做循環？

謝謝！

Answer 1

set.seed(123) 
alive <- vector("list", 20) 
mu <- 0.1 
n <- 10 

alive[[1]] <- rbinom(n, 100, exp(-mu)) 
for(i in 2:20) 
    alive[[i]] <- rbinom(n, alive[[i-1]], exp(-mu)) 

我改名變量sample到n以避免與常用的功能sample()混亂。

Answer 2

set.seed(123) 
alive <- 100  
mu <- 0.1 
sample <- 10 
alive1 <- rbinom(sample,alive,exp(-mu)) 


for (i in 2:20) 
{ 
    assign(
     paste0("alive",i), 
     rbinom(10,get(paste0("alive",(i-1))),exp(-mu)) 
    ) 
} 

或@把它在一個列表Backlin的建議 -

set.seed(123) 
alive <- 100  
mu <- 0.1 
sample <- 10 
Aliveset <- vector(mode = "list", length = 20) 
Aliveset[[1]] <- rbinom(sample,alive,exp(-mu)) 
for (i in 2:20) 
{ 
Aliveset[[i]] <- rbinom(10,Aliveset[[i-1]],exp(-mu)) 
}