R 中某一列的总百分比变化？

假设 LA 的 0.93 是打字错误：

library(dplyr)

city_pop <- tibble(
  city = c("NYC", "NYC", "LA", "LA"), 
  percentpopchange = c("0.1", "0.2", "0.5", "0.3"), 
  year = c(1980, 1981, 1980, 1981)
)

city_pop %>% 
  group_by(city) %>% 
  summarize(
    total_percentpopchange = as.character(
      prod(1 + as.numeric(percentpopchange)) - 1
  ))
#> # A tibble: 2 × 2
#>   city  total_percentpopchange
#>   <chr> <chr>                 
#> 1 LA    0.95                  
#> 2 NYC   0.32

^{创建于 2024-09-03，使用reprex v2.1.1}

如果我们假设 1980 年纽约州的人口为 1,000 人，那么 1980 年和 1981 年之后其人口变化如下：

1980 年：1000 美元 + 1000 美元*0.1 = 1100 美元

1981 年：1100 美元 + 1100 美元*0.2 = 1320 美元。

因此，在 1980 年和 1981 年期间，纽约州的人口增加了：（1320 美元 - 1000 美元）/1000 美元 = 0.32。

下面粘贴了一些计算该问题的 R 代码：

city_pop<-tibble(city=c("NYC", "NYC", "LA", "LA"), 
                 percentpopchange=c("0.1", "0.2", "0.5", "0.3"), 
                 year=c(1980, 1981, 1980, 1981))
city_pop$percentpopchange <- as.numeric(city_pop$percentpopchange)

cum_city_pop_final <- tibble(city=unique(city_pop$city))

unique_cities <- unique(city_pop$city)
total_percent_pop_changes <- c()
for (i in 1:length(unique(city_pop$city))){
  pop <- 100
  for (j in 1:length(unique(city_pop$year))){
    pop <- pop + pop * city_pop$percentpopchange[which(city_pop$city==unique(city_pop$city)[i] & city_pop$year==unique(city_pop$year)[j])]
    if (j==length(unique(city_pop$year))){
      total_percent_pop_changes <- c(total_percent_pop_changes, (pop-100)/100)
    }
  }
}

cum_city_pop_final$total_percentpopchange <- total_percent_pop_changes

顺便说一句，我认为 1980 年至 1981 年洛杉矶的总人口变化百分比应该是 0.95，而不是 0.93。