ccp_closed0: Cohort Component Population Projection Model based on a...

In gjabel/demoproj: Demographic Projection Models

Description

A cohort component projection model based on a closed population,

N(t+1) = L[t, t+1] N(t)

where the Leslie matrix, eqnL, is created given user defined age specific fertility and survivorship rates.

Usage

ccp_closed0(n = NULL, x = NULL, p = NULL, Nx_f = NULL, Nx_m = NULL,
  sx_f = NULL, sx_m = NULL, fx = NULL, sn_f = NULL, sn_m = NULL,
  sex_ratio = 1/(1 + 1.05), tidy_output = TRUE, age_lab = x,
  gender_lab = c("Female", "Male"), ...)

ccp_closed(n = NULL, x = NULL, p = NULL, Nx_f = NULL, Nx_m = NULL,
  sx_f = NULL, sx_m = NULL, fx = NULL, sn_f = NULL, sn_m = NULL,
  sex_ratio = 1/(1 + 1.05), tidy_output = TRUE, age_lab = x,
  gender_lab = c("Female", "Male"), ...)

Arguments

n	Numeric value for the number of projection steps.
x	Vector containing a character string of age group labels.
p	Numeric value for step size of the population projection.
Nx_f, Nx_m	Vectors containing numeric values of the initial female and male population size in each age group (x).
sx_f, sx_m, fx	Vectors containing numeric values of the age specific female and male survival and fertility rates. If ccp_closed0 is used a single vector is required. If ccp_closed is used a matrix of period specific rates is required, where rows represent age groups (females in top rows, males in bottom rows) and columns future periods. If a single vector is passed to ccp_closed a matrix based on constant assumptions in all future rates will be constructed.
sn_f, sn_m, sex_ratio	Numeric value of the female and male survivorship of new-born babies from birth to the end of the interval and the sex ratio at birth of new-born babies. If ccp_closed0 is used a single value is required. If ccp_closed is used a vector of period specific values is required. If a single value is passed to ccp_closed a vector based on constant assumptions in all future rates will be constructed.
tidy_output	Logical value to indicate if projection output should be in a tidy data format (TRUE, the default) or as a matrix where rows represent age and gender groups and columns the projection year.
age_lab, gender_lab	Vector containing a character string of age and gender group labels. Only used if projection output is in a tidy data format. See tidy_pp for more details.
...	Additional arguments passed to tidy_pp to build a tidy data frame (if chosen from tidy_output).

Value

Projected populations by age and gender for n future steps, given the age specific fertility and survivorship rates. Depending on the tidy_output value the projections will be returned as either a matrix or a tibble. Both versions contain the initial population sizes given in Nx.

ccp_closed0 produces population projections based strictly on constant future rates.

ccp_closed produces population projections based non-constant future rates.

Examples

df0 <- sweden1993

# matrix output
ccp_closed0(n = 5, x = df0$x, p = 5, Nx_f = df0$Nx_f, Nx_m = df0$Nx_f,
            sx_f = df0$sx_f, sx_m = df0$sx_f,
            fx = df0$fx, sn_f = df0$Lx_f[1]/(5*100000), sn_m = df0$Lx_m[1]/(5*100000),
            tidy_output = FALSE)

# tidy data frame output
ccp_closed0(n = 5, x = df0$x, p = 5, Nx_f = df0$Nx_f, Nx_m = df0$Nx_f,
            sx_f = df0$sx_f, sx_m = df0$sx_f,
            fx = df0$fx, sn_f = df0$Lx_f[1]/(5*100000), sn_m = df0$Lx_m[1]/(5*100000),
            year0 = 1993, age_lab = df0$age)

# setting up non-constant future age specific fertility rates
ff <- matrix(df0$fx, nrow = length(df0$fx), ncol = 5)
ff <- sweep(ff, 2, seq(from = 1, to = 1.5, length = 5), "*")
# tfr increase
5 * colSums(ff)
# run projection with increasing fx, sx remains constant
ccp_closed(n = 5, x = df0$x, p = 5, Nx_f = df0$Nx_f, Nx_m = df0$Nx_f,
            sx_f = df0$sx_f, sx_m = df0$sx_f,
            fx = ff, sn_f = df0$Lx_f[1]/(5*100000), sn_m = df0$Lx_m[1]/(5*100000),
            tidy_output = FALSE)

gjabel/demoproj documentation built on May 17, 2019, 6:01 a.m.