Skip to contents

The longevity package proposes estimation routines for modeling excess lifetime. Core functionalities include maximum likelihood estimation for parametric models (exponential, Gompertz, Weibull, generalized Pareto, extended generalized Pareto, piecewise generalized Pareto), threshold selection plots for survival data, nonparametric maximum likelihood estimation, profile likelihood estimation for the endpoint of the distribution of exceedances.


You can install the development version of longevity from GitHub with:

# install.packages("devtools")


  • maximum likelihood estimation routines
  • simulation of left-truncated and right-truncated/right-censored data
  • hazard plots with profile-likelihood based confidence intervals
  • threshold selection diagnostics with profile and Wald pointwise confidence intervals
  • quantile-quantile plots
  • likelihood ratio tests for nested models
  • likelihood ratio tests for a categorical explanatory
  • score and likelihood ratio tests for piecewise generalized Pareto distribution, extending Northrop and Coleman (2014), with p-value paths
  • nonparametric maximum likelihood estimate of the distribution function with arbitrary truncation and censoring using the EM algorithm of Turnbull (1976) - C++ implementation.
  • profile likelihood for generalized Pareto (endpoint)
  • hazard functions for all parametric models, using helife