# defining Population Expansion Index PEI <- function(X, CLS=1, MODE=1) { # Xはゼロ歳からの年齢階級別人口ベクトル。 # CLSは年齢階級幅。デフォルトは年齢各歳。1か5か10のみ取れる。 # MODE=1でふくらみ指数改,他のMODEでオリジナルのふくらみ指数 N1 <- 20/CLS + 1 - ifelse(MODE==1, 0, 5/CLS) N2 <- N1 + 20/CLS - 1 D1 <- 10/CLS + 1 - ifelse(MODE==1, 0, 5/CLS) D2 <- D1 + 10/CLS - 1 D3 <- D1 + 30/CLS D4 <- D2 + 30/CLS return(sum(X[N1:N2]) / (sum(X[D1:D2]) + sum(X[D3:D4])) * 100) }