OPAG_simple: redistribute an open age group count over higher ages...

In timriffe/DemoTools: Standardize, Evaluate, and Adjust Demographic Data

redistribute an open age group count over higher ages proportional to an arbitrary standard

Description

This method could be useful whenever a reasonable standard is available. At present the standard must be supplied by the user.

Usage

OPAG_simple(Pop, Age, OAnow = max(Age), StPop, StAge, OAnew = max(StAge))

Arguments

Pop	numeric vector of population counts
Age	integer vector of single age lower bounds
OAnow	integer. The lower age bound above which counts will be redistributed
StPop	numeric vector of standard population counts
StAge	integer vector of single age lower bounds for the standard population
OAnew	integer. The desired new open age, must be no higher than max(StAge).

Details

In this implementation both the original population counts and the standard must be in single ages.

References

\insertRef

PASDemoTools

Examples

 Pop        <- c(38129,38382,38824,39275,39500,37304,35152,
 34061,33911,32875,31599,30376,29822,29691,28765,
 28695,28917,28203,29209,30316,29062,26977,26577,
 27727,28599,30513,31774,32347,34093,33736,32085,
 30807,28279,26873,25612,23503,22207,21388,20122,
 18014,15626,15006,14158,11195,7931,7640,9053,
 13276,17226,18918,17697,18424,17723,16706,14410,
 13342,14787,15183,15727,16045,14777,14267,13102,
 10866,9311,6933,5030,3785,3551,2848,3080,
 2874,2368,2681,3165,3010,3009,2721,2705,
 2492,2244,1971,1644,1565,1307,5027)
 Age        <- 0:85
 # standard pop taken from ages 55+
 StPop      <- c(6258,6177,6089,5995,5894,5787,5672,5552,
 5423,5286,5140,4985,4824,4652,4477,4293,
 4107,3912,3712,3502,3282,3055,2823,2591,
 2360,2138,1921,1710,1502,1297,1098,910,
 741,592,463,353,265,192,137,95,
 63,42,26,17,9,13)

 StAge      <- 55:100

PopExtended <- OPAG_simple(
 		Pop = Pop,
 		Age = Age,
 		StPop = StPop,
 		StAge = StAge)

## Not run: 
 plot(Age, Pop, type = 'l',xlim=c(80,100),ylim=c(0,1e4))
lines(0:100, PopExtended, col = "red", lty = 2)

## End(Not run)
stopifnot((sum(PopExtended[86:101]) - Pop[86]) == 0)

timriffe/DemoTools documentation built on Oct. 14, 2024, 12:53 p.m.