dfeehan / mortfit

Fit mortality hazard models, with special attention to humans at old ages.

mortfit

Fit mortality hazard models, with special attention to humans at old ages.

NOTES

• 20160526 - this depends on the version of RcppFaddeeva that is on github and not the CRAN version. (The github version exports the Rcpp headers properly.) So it's necessary to

  install_github("baptiste/RcppFaddeeva@5fb5c32")

  in order to build this version

(NB: The issue opened by Tim Riffe on github has a solution)

To add a hazard function

Example steps used to create the Beard hazard object:

• create mortalityhazard-beard.R

  • include code to create new mortalityHazard object

• calculate

  1. hazard to probability function
  2. gradients

• Using Sage, I produced C++ code to calculate (1) and (2) above. I put the results in mortalityhazard-beard.cpp. The function names are:

  • mortalityhazard_beard_cpp - calculates the mortality hazard, given params
  • mortalityhazard_to_prob_beard_cpp - given params and ages, calculates cond prob of death
  • mortalityhazard_beard_binomial_grad_cpp - calculates the gradient of the binomial likelihood
  • beard_partial_alpha_partI - calculates first part of gradient wrt alpha
  • beard_partial_alpha_partII - calculates second part of gradient wrt alpha
  • similar functions to calculate gradient wrt other params

mortfit's Issues

add Remotes to DESCRIPTION

Hi! Looking forward to checking this out.

You can handle the github package dependency issue by adding to your DESCRIPTION file:
Remotes: baptiste/RcppFaddeeva
or similar (maybe the commit too if necessary).