epitools (version 0.09)

ageadjust.indirect: Age standardization by indirect method

Description

Calculates age standardized (adjusted) rates and "exact" confidence intervals using the indirect method

Usage

ageadjust.indirect(count, pop, stdcount, stdpop, stdrate = NULL, conf.level =
0.95)

Arguments

count

vector of age-specific count of events

pop

vector of age-specific person-years or population estimates

stdcount

vector of age-specific standard counts

stdpop

vector of age-specific standarad population

stdrate

vector of age-specific standard rates

conf.level

confidence level (default = 0.95)

Value

• $sir
  observed, expected, standardized incidence ratio, and confidence interval
• $rate
  crude.rate, adjusted rate, and confidence interval

Details

Calculates age standardized (adjusted) rates and "exact" confidence intervals using the indirect method

References

Anderson RN, Rosenberg HM. Age Standardization of Death Rates: Implementation of the Year 200 Standard. National Vital Statistics Reports; Vol 47 No. 3. Hyattsville, Maryland: National Center for Health Statistics. 1998, pp. 13-19. Available at http://www.cdc.gov/nchs/data/nvsr/nvsr49/nvsr49_09.pdf.

Selvin S. Statistical Analysis of Epidemiologic Data (Monographs in Epidemiology and Biostatistics, Vol 17), Oxford University Press, (1991)

See Also

See also ageadjust.direct

Examples

dth60 <- c(141, 926, 1253, 1080, 1869, 4891, 14956, 30888, 41725, 26501,
5928)
pop60 <- c(1784033, 7065148, 15658730, 10482916, 9939972, 10563872, 9114202,
6850263, 4702482, 1874619, 330915)
dth40 <- c(45, 201, 320, 670, 1126, 3160, 9723, 17935, 22179, 13461, 2238)
pop40 <- c(906897, 3794573, 10003544, 10629526, 9465330, 8249558, 7294330,
5022499, 2920220, 1019504, 142532)
rate60 <- dth60/pop60
rate40 <- dth40/pop40
tab <- array(c(dth60, pop60, round(rate60*100000,1), dth40, pop40,
round(rate40*100000,1)), c(11,3,2))
agelabs <- c("<1", "1-4", "5-14", "15-24", "25-34", "35-44", "45-54",
"55-64", "65-74", "75-84", "85+")
dimnames(tab) <- list(agelabs, c("Deaths", "Population", "Rate"),
c("1960", "1940"))
tab
ageadjust.indirect(count = dth40, pop = pop40, stdcount = dth60, stdpop
= pop60)