map_maker: Create two dimensional hazard rate from...

Description Usage Arguments Value Examples

Description

Create two dimensional hazard rate from Age-Cohort-Interaction model

Usage

1
2
map_maker(mu, age_coef, cohort_coef, age_seq, cohort_seq,
  islet_ls = NULL)

Arguments

mu

Constant part of the log-hazard

age_coef

Age-effect parameter sequence

cohort_coef

Cohort-effect parameter sequence

age_seq

Age discretizatioin vector

cohort_seq

Cohort discretizatioin vector

islet_ls

List of parameters of the interaction terms.

Value

The interactions terms are densities of two dimensional normal distributions. islet_ls gives the mean, variance-covariance matrix and multiplicative coefficients of each interaction term.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## Not run: 
library(lattice)
mu <- log(1e-2)
age_seq <- seq(0, 100, length.out = 100)
cohort_seq <- seq(1900, 1900 + age_seq[101], length.out = age_length + 1)
age_coef <- seq(0, 2, length.out = age_length - 1)
cohort_coef <- seq(0, 0.8, length.out = age_length - 1)
islet_ls = list(list(mean = c(1945, 45), sigma = 100 * diag(2), coef = 40))
hazard <- map_maker(mu, age_coef, cohort_coef, age_seq, cohort_seq, islet_ls = islet_ls)
persp(hazard$sheet)

## End(Not run)

goepp/hazreg documentation built on June 7, 2019, 4:03 a.m.