R语言准备金流量三角形下三角0值批量替换为NA或空白的方法

核心逻辑 #

标准流量三角形以事故年为行、进展年为列，n行n列结构中，满足行索引 + 列索引 > n + 1的位置即为未观测的下三角区域，通常默认输出为0值，可通过定位该区域批量替换目标值。

实现方法 #

方法1：基础R实现（无需额外安装包） #

# 构造示例5*5流量三角形，下三角默认填充0
triangle <- matrix(c(
  100, 150, 180, 200, 210,
  120, 170, 205, 220, 0,
  110, 165, 195, 0, 0,
  130, 185, 0, 0, 0,
  140, 0, 0, 0, 0
), nrow = 5, byrow = TRUE, dimnames = list(paste0("AY",2019:2023), paste0("DY",1:5)))

n <- nrow(triangle)
# 批量替换下三角0为NA
for(i in 1:n){
  triangle[i, (n - i + 2):n] <- NA
}
# 如果需要替换为空白，替换为以下代码：
# triangle[i, (n - i + 2):n] <- ""

方法2：ChainLadder包实现（准备金评估常用工具包） #

针对准备金评估场景，直接使用行业通用包的内置逻辑处理，无需手动定位下三角：

# 安装加载包（已安装可跳过安装步骤）
install.packages("ChainLadder")
library(ChainLadder)

# 构造长格式准备金原始数据示例
raw_data <- data.frame(
  ay = rep(2019:2023, each = 5),
  dy = rep(1:5, 5),
  loss = c(100,150,180,200,210,120,170,205,220,0,110,165,195,0,0,130,185,0,0,0,140,0,0,0,0)
)

# 转换为标准三角格式，自动将下三角无效0值替换为NA
triangle_cl <- as.triangle(raw_data, origin = "ay", dev = "dy", value = "loss")

注意：如果已有现成的矩阵格式流量三角，可直接转换：as.triangle(your_triangle_matrix)，输出结果自动将下三角0值替换为NA。

内容的提问来源于stack exchange，提问作者Tiberius Nyabigo