quarterly_variables: Data frame of quarterly variables

View source: R/quarterly_variables.R

quarterly_variablesR Documentation

Data frame of quarterly variables

Description

Computes punctual risk coordinates in the Lexis diagram and quarterly biometric variables of a population.

Usage

quarterly_variables(
  date.birth,
  date.event,
  random.b = TRUE,
  random.e = TRUE,
  constant.age.year = FALSE
)

Arguments

date.birth

A character vector with the dates of birth in format either "yyyy-mm-dd" or "yyyy-mm-dd hour:min:secs" (for instance, "2016-01-20 12:00:00") of a population. If "hour:min:secs" is omitted the function imputes either "12:00:00", if 'random.b = FALSE', or a random hour by default.

date.event

A character vector with the dates of events in format either "yyyy-mm-dd" or "yyyy-mm-dd hour:min:secs" (for instance, "2016-01-20 12:00:00") of a population. If "hour:min:secs" is omitted the function imputes either "12:00:00", if 'random.e = FALSE', or a random hour, by default. This vector must have either length 1, when the aim is to compute the exact age or the (1x1-Lexis) age coordinate of all the members of the population in the same temporal point or the same length as 'date.birth' when the aim is to compute for each member of the population the exact age or the (1x1-Lexis) age coordinate in the moment of the event (e.g., death).

random.b

A 'TRUE/FALSE' argument indicating whether the exact moment ("hour:min:secs") when the birth occurs within the day is randomly selected. If TRUE, this overwrites "hour:min:secs" in 'date.birth' even if those have been declared. By default, TRUE.

random.e

A 'TRUE/FALSE' argument indicating whether the exact moment ("hour:min:secs") when the event occurs within the day is randomly selected. If TRUE, this overwrites "hour:min:secs" in 'date.event' even if those have been declared. By default, TRUE.

constant.age.year

A 'TRUE/FALSE' argument indicating whether the length of the year should be constant, 365.25 days, or variable, depending on the time lived for the person in each year since her/his dates of birth and event. By default, FALSE. The advantage of using a non-constant (person-dependent) length of year is congruence when estimating time exposed at risk: in each year the time exposed along the time and age axes will coincide.

Value

A data.frame with the following components:

coord.age

Time elapsed, measure in years, between the last birthday and the date when the event happens.

coord.time

Time coordinate: time elapsed, measure in years, between the begining of the year and the date when the event happens.

age.last.birthday

The integer age at last birthday.

exact.age.at.event

Time elapsed, measure in years, between the dates of birth and event.

quarter.age

Age quarter when the event happens.

quarter.calendar

Calendar (time, season) quarter to which the time exposed at risk corresponds.

year

Year when the event happens.

Note

In the age axis, the length of the years are assumed either constant 365.25 days ('constant.age.year = TRUE') or variable ('constant.age.year = FALSE'), depending on the person. In the time axis, the length of the year is either 365 in non-leap years and 366 in leap years. The advantage of using a non-constant (person-dependent) length of year in the age axis is that in each year the lengths of the years when computing 'coord.age' and 'coord.time' in both axis are equal.

Author(s)

Jose M. Pavia pavia@uv.es

Josep Lledo josep.lledo@uv.es

References

Pavia, JM and Lledo, J (2022). Estimation of the Combined Effects of Ageing and Seasonality on Mortality Risk. An application to Spain. *Journal of the Royal Statistical Society, Series A (Statistics in Society)*, 185(2), 471-497. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/rssa.12769")}

Examples

dates.b <- c("1920-05-13", "1999-04-12", "2019-01-01")
dates.e <- c("2002-03-23", "2009-04-12", "2019-01-01")
quarterly_variables(dates.b, dates.e)

qlifetable documentation built on April 19, 2023, 5:09 p.m.