Computation of nonparametric convex hazard estimators via profile methods. (English) Zbl 1161.62014

Summary: This paper proposes a profile likelihood algorithm to compute the nonparametric maximum likelihood estimator of a convex hazard function. The maximisation is performed in two steps: First the support reduction algorithm is used to maximise the likelihood over all hazard functions with a given point of minimum (or antimode). Then it is shown that the profile (or partially maximised) likelihood is quasi-concave as a function of the antimode, so that a bisection algorithm can be applied to find the maximum of the profile likelihood, and hence also the global maximum. The new algorithm is illustrated using both artificial and real data, including lifetime data for Canadian males and females.


62G05 Nonparametric estimation
62N02 Estimation in survival analysis and censored data
65C60 Computational problems in statistics (MSC2010)
62A09 Graphical methods in statistics
Full Text: DOI Link


[1] Grenander U., Skand. Aktuarietidskr. 39 pp 125– (1957)
[2] Bray, T., Crawford, G. and Proschan, F. 1967. ”Maximum likelihood estimation of a U-shaped failure rate function”. Seattle, WA: Mathematics Research Laboratory, Boeing Scientific Research Laboratories. Mathematical Note 534, Available athttp://www.stat.washington.edu/jaw/RESEARCH/OLD-PAPERS-OTHERS/UMLE.pdf
[3] Prakasa Rao B. L.S., Sankhyā Ser. A. 31 pp 23– (1969)
[4] DOI: 10.1214/009053605000000138 · Zbl 1072.62023
[5] Baraud Y., Probab. Theory Related Fields
[6] Brunel E., Math. Methods Statist. 15 pp 233– (2006)
[7] DOI: 10.1007/s00440-003-0259-1 · Zbl 1019.62079
[8] Borwein J. M., Convex Analysis and Nonlinear Optimisation (2000)
[9] Jankowski, H. and Wellner, J. A. 2007. ”Nonparametric estimation of a convex bathtub-shaped hazard function”. Department of Statistics, University of Washington. Tech. Rep. 521 · Zbl 1200.62025
[10] Jankowski, H., Wang, X., McCauge, H. and Wellner, J. 2008. ”convexHaz: R functions for convex hazard rate estimation”. software available athttp://www.r-project.org
[11] Jankowski, H. and Wellner, J. A. 2008. ”Computation of nonparametric convex hazard estimators via likelihood profile methods”. Department of Statistics, University of Washington. Tech. Rep. 542 · Zbl 1161.62014
[12] DOI: 10.1111/j.1467-9469.2007.00588.x · Zbl 1199.65017
[13] DOI: 10.1214/aos/1176345890 · Zbl 0523.62074
[14] Dümbgen, L., Hüsler, A. and Rufibach, K. 2007. ”Active set and EM algorithms for log-concave densities based on complete and censored data”. University of Bern.
[15] Bazaraa M. S., Nonlinear programming; Theory and Algorithms, 3. ed. (2006) · Zbl 1140.90040
[16] DOI: 10.1016/S0378-3758(97)89710-8 · Zbl 0900.62537
[17] DOI: 10.2307/1266340
[18] Cox D. R., The Statistical Analysis of Series of Events (1966) · Zbl 0148.14005
[19] Proschan, F. and Pyke, R. Tests for monotone failure rate. Proceedings of the Fifth Berkeley Symposiom of Mathematical Statistics and Probability. 1965/66, Berkeley, CA. Vol. III, pp.293–312. Berkeley, CA: Physical Sciences, University of California Press.
[20] Statistics Canada, Complete life table, Canada, 2000 to 2002, females and males, (2002). Available online:http://www.statcan.ca/english/freepub/84-537-XIE/tables/txttables/cam.txt(and caf.txt)
[21] DOI: 10.1016/j.jtbi.2006.11.011
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.