LARisk (version 1.0.0)

LAR_group: Average Estimated Lifetime Attributable Risk by Group

Description

LAR_group is used to estimate lifetime attributable radiation-related cancer risk by group.

Usage

LAR_group(
  data,
  pid,
  group,
  basedata,
  sim = 300,
  seed = 99,
  current = as.numeric(substr(Sys.Date(), 1, 4)),
  ci = 0.9,
  weight = NULL,
  DDREF = TRUE,
  basepy = 1e+05
)

Arguments

data

data frame containing demographic information and exposure information. See 'Details'.

pid

a vector which distinguish each person.

group

a vector or list of vectors which distinguish each group.

basedata

a list of the data of lifetime table and incidence rate table. The first element is lifetime table and the second is incidence rate table.

sim

number of iteration of simulation.

seed

a random seed number.

current

a current year. default is year of the system time.

ci

confidence level of the confidence interval.

weight

a list containing the value between 0 and 1 which is a weight on ERR model. See 'Details'.

DDREF

logical. Whether to apply the dose and dose-rate effectiveness factor.

basepy

number of base person-years

Value

LAR_group returns an object of multiple classes "LAR_group", "LAR". An object of class LAR_group is a list of LAR class objects which names of elements are group of each groups.

References

Berrington de Gonzalez, A., Iulian Apostoaei, A., Veiga, L., Rajaraman, P., Thomas, B., Owen Hoffman, F., Gilbert, E. and Land, C. (2012). RadRAT: a radiation risk assessment tool for lifetime cancer risk projection. Journal of Radiological Protection, 32(3), pp.205-222.

National Research Council (NRC) and Committee to Assess Health Risks from Exposure to Low Levels of Ionizing Radiation (2005) Health Risks from Exposure to Low Levels of Ionizing Radiation: BEIR VII Phase 2 (Washington, DC: National Academy of Sciences)

Examples

# NOT RUN {
## example with lifetime and incidence rate table in 2010 Korea.
lar1 <- LAR_group(nuclear, pid=nuclear$ID, group=nuclear$distance,
                 basedata = list(life2010, incid2010))
summary(lar1)

lar2 <- LAR_group(nuclear, pid=nuclear$ID, group=list(nuclear$sex, nuclear$distance),
                 basedata = list(life2010, incid2010))
summary(lar2)


# }