Rで年齢調整

# データ作成

# A市（高齢者が多い）
city_A_data <- data.frame(
  age_group = c("0-19", "20-39", "40-59", "60+"),
  population = c(15000, 20000, 30000, 35000), # 合計10万人、60+が多い
  deaths = c(15, 40, 150, 700) # 高齢者の死亡者数が非常に多い
)

# B市（若者が多い）
city_B_data <- data.frame(
  age_group = c("0-19", "20-39", "40-59", "60+"),
  population = c(30000, 40000, 20000, 10000), # 合計10万人、若年層が多い
  deaths = c(36, 96, 120, 250) # 全体的に死亡率がA市より高めに設定
)

# 基準人口（比較の基準とする標準的な人口構成）
standard_pop_data <- data.frame(
  age_group = c("0-19", "20-39", "40-59", "60+"),
  standard_population = c(25000, 35000, 25000, 15000) # 合計10万人
)

cat("--- A市のデータ ---\n\n")
print(city_A_data)
cat("\n--- B市のデータ ---\n\n")
print(city_B_data)
cat("\n--- 基準人口のデータ ---\n\n")
print(standard_pop_data)

# A市の粗死亡率を計算
crude_rate_A <- sum(city_A_data$deaths) / sum(city_A_data$population)

# B市の粗死亡率を計算
crude_rate_B <- sum(city_B_data$deaths) / sum(city_B_data$population)

# 結果を表示（人口10万人あたりに換算）
cat("--- 年齢調整なし（粗死亡率）の比較 ---\n\n")
cat(sprintf("A市の粗死亡率: %.1f (人口10万人あたり)\n", crude_rate_A * 100000))
cat(sprintf("B市の粗死亡率: %.1f (人口10万人あたり)\n", crude_rate_B * 100000))

# 年齢調整死亡率を計算する関数を定義
calculate_adjusted_rate <- function(city_data, standard_pop) {
  # 1. 年齢階級別の死亡率を計算
  city_data <- city_data %>%
    mutate(age_specific_rate = deaths / population)

  # 基準人口データと結合
  merged_data <- left_join(city_data, standard_pop, by = "age_group")

  # 2. 期待死亡者数を計算
  merged_data <- merged_data %>%
    mutate(expected_deaths = age_specific_rate * standard_population)

  # 3. 期待死亡者数を合計
  total_expected_deaths <- sum(merged_data$expected_deaths)

  # 4. 基準人口の総人口で割る
  total_standard_pop <- sum(standard_pop$standard_population)
  adjusted_rate <- total_expected_deaths / total_standard_pop

  return(adjusted_rate)
}

# A市とB市の年齢調整死亡率を計算
adjusted_rate_A <- calculate_adjusted_rate(city_A_data, standard_pop_data)
adjusted_rate_B <- calculate_adjusted_rate(city_B_data, standard_pop_data)

# 結果を表示（人口10万人あたりに換算）
cat("--- 年齢調整あり（年齢調整死亡率）の比較 ---\n\n")
cat(sprintf("A市の年齢調整死亡率: %.1f (人口10万人あたり)\n", adjusted_rate_A * 100000))
cat(sprintf("B市の年齢調整死亡率: %.1f (人口10万人あたり)\n", adjusted_rate_B * 100000))