age_adjusted: Returns age-adjusted proportions with its confidence interval

In mechkar: Useful Tools for Scientific Research

Description

This function is intended to calculate the age-adjusted rates given a dataset that includes age and an binary outcome. It returns the crude rates with their confidence interval and the age-adjusted rates with its confidence interval.

Usage

age_adjusted(dataset,outcome,age,agemin=0,agemax=130,source="who",alpha=0.05)

Arguments

dataset	the name of the dataframe
outcome	the name of the outcome variable
age	the name of the age variable
agemin	the minimum age to be used to calculate the weights (default=0)
agemax	the maximum age to be used to calculate the weights (default=130)
source	the name of the standard to be used for the weights (see the details section). Default: 'who'
alpha	the probability of a type I error. Default: 0.05 for the 95% CI (1-alpha)

Details

This function return a list containing nine values (see values section). It includes three types of standards: the European (euro), the American (us) and the World Health Organization (who). The confidence interval (CI) is calculated using the binomial probabilities using the binconf function of the Hmisc package. The default is to calculate the 95% CI (using an alpha of 0.05).

Value

This function returns a list containing the following fields:

• outcome - the name of the outcome variable

• population - the total population on the dataset

• events - the number of events for the outcome

• crude - a list with the crude rates. This contains a list we three elements:

  • rate - the value of the crude rate

  • CImin - the value of the lower CI

  • CImax - the value of the higher CI

• weighted - a list with the weighted rates. This contains a list we three elements:

  • rate - the value of the crude rate

  • CImin - the value of the lower CI

  • CImax - the value of the higher CI

Author(s)

Tomas Karpati M.D.

References

"The WHO Standard Population is defined to reflect the average age structure of the world's population over the next generation, from the year 2000 to 2025." https://www.nrscotland.gov.uk/files//statistics/age-standardised-death-rates-esp/age-standard-death-rates-background.pdf

American standard (2000): https://www.cdc.gov/nchs/data/statnt/statnt20.pdf

World (WHO 2000-2025) Standard: https://seer.cancer.gov/stdpopulations/world.who.html

Examples

age <- floor(runif(150,min=30,max=65))
outcome <- round(runif(150,min=0,max=1),0)
df <- data.frame(age=age,disease=outcome)
aa <- age_adjusted(df,"disease","age",agemin=30,agemax=65,source="who",alpha=0.0)

mechkar documentation built on March 13, 2020, 2:30 a.m.