Почему среднее значение стабильного распределения по возрасту совпадает с ожидаемой продолжительностью жизни?

Если предположить, что уровень опасности (или смертности) в популяции постоянен, это приводит к экспоненциальному распределению времени выживания. Обратная величина коэффициента смертности дает вам ожидаемую продолжительность жизни. Ниже я моделирую популяцию индивидуумов (отслеживая возраст каждого индивида), которые умирают с постоянной скоростью (т.е. независимо от их возраста). Почему среднее значение стабильного распределения по возрасту совпадает с ожидаемой продолжительностью жизни?

pop.size <- 5000

mort.rate <- 0.02

mean.age.dist <- c()

ages.in.pop <- rep(0, pop.size)

for(i in 1 : 5000)
 {
   ages.in.pop <- ages.in.pop + 1

   death.vec <- rbinom(pop.size, 1, mort.rate)

   ages.in.pop[which(death.vec == 1)] <- 0

   mean.age.dist[i] <- mean(ages.in.pop)
 }

plot(mean.age.dist, ylab='mean age in population', xlab='iteration', type      ='l')
print('life expectancy')
print(1 / mort.rate)

0 ответов

Другие вопросы по тегам