在R中按分组基于基准年和增长率计算变量连续值

针对你的需求，这里提供三种常用的R语言实现方案，均能按Scode分组，基于初始GDP值逐年计算连续GDP：

方案1：使用dplyr（tidyverse风格） #

适合习惯 tidy 语法的用户，代码可读性强：

library(dplyr)

# 确保数据按国家和年份排序，避免计算顺序出错
dat <- dat %>%
  arrange(Scode, Year) %>%
  group_by(Scode) %>%
  mutate(
    # 提取每个国家的初始GDP基准值
    base_gdp = first(OfficialGDP),
    # 计算累积乘数：基准年设为1（无需相乘），后续年份按顺序累积乘积
    cum_multiplier = cumprod(ifelse(row_number() == 1, 1, GDPMultiplier)),
    # 生成连续GDP值
    OfficialGDP = base_gdp * cum_multiplier
  ) %>%
  # 移除临时辅助列
  select(-base_gdp, -cum_multiplier) %>%
  ungroup()

方案2：使用data.table（高效大数据处理） #

适合处理大规模数据集，运行效率远高于基础循环：

library(data.table)

# 将数据框转换为data.table格式
setDT(dat)

# 按国家分组，按年份排序后计算连续GDP
dat <- dat[order(Scode, Year), 
           OfficialGDP := first(OfficialGDP) * cumprod(ifelse(seq_len(.N) == 1, 1, GDPMultiplier)),
           keyby = Scode]

方案3：基础R分组循环（贴合你的原始思路） #

完全匹配你设想的循环逻辑，按组逐行计算：

# 先确保数据按国家和年份升序排列
dat <- dat[order(dat$Scode, dat$Year), ]

# 获取所有唯一的国家代码
unique_scodes <- unique(dat$Scode)

# 遍历每个国家分组计算
for (sc in unique_scodes) {
  # 定位当前国家的所有行索引
  group_rows <- which(dat$Scode == sc)
  # 从第二行开始，用上一行GDP乘当年乘数
  for (i in 2:length(group_rows)) {
    prev_row <- group_rows[i-1]
    curr_row <- group_rows[i]
    dat$OfficialGDP[curr_row] <- dat$OfficialGDP[prev_row] * dat$GDPMultiplier[curr_row]
  }
}

结果验证 #

以AGO为例，计算结果会和你预期一致：

• 1991年GDP：57.96 * 1.009 ≈ 58.48
• 1992年GDP：58.48 * 0.94 ≈ 54.97
• 后续年份以此类推

内容的提问来源于stack exchange，提问作者BaharrÖ