r-计算一段时间内已知死亡率的预期幸存者人数

Question

总人数为N，日死亡率为p1。我想计算几天后的幸存者人数。

day1 : N - N*p1 
day2 : (N - N*p1) -  (N - N*p1)*p1
day3 : ((N - N*p1) -  (N - N*p1)*p1) - ((N - N*p1) -  (N - N*p1)*p1)*p1
...

这就是我到目前为止所做的，开始获得第1天的数字。如有建议，不胜感激。

df <- data.frame(day = c(1:30))  
N <- 1000
p1 <- 0.06
apply(df,1, function(x) N-N*p1) 

Answer 1

基本上，您可以使用复利方程。

foo = function(N, p1, d){
for (i in 1:d){
N = N - N * p1
}
N
}
foo(1000, 0.06, 30)

或

bar = function(N, p1, d){
N * (1 - p1)^d
}
bar(1000, 0.06, 1:30)

或者如果你只是想要一个绘图

curve(1000 * (1 - 0.06) ^ x, 0, 30, n = 31)