Statistical Science

Inference and Modeling with Log-concave Distributions

Guenther Walther

Full-text: Open access


Log-concave distributions are an attractive choice for modeling and inference, for several reasons: The class of log-concave distributions contains most of the commonly used parametric distributions and thus is a rich and flexible nonparametric class of distributions. Further, the MLE exists and can be computed with readily available algorithms. Thus, no tuning parameter, such as a bandwidth, is necessary for estimation. Due to these attractive properties, there has been considerable recent research activity concerning the theory and applications of log-concave distributions. This article gives a review of these results.

Article information

Statist. Sci., Volume 24, Number 3 (2009), 319-327.

First available in Project Euclid: 31 March 2010

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Nonparametric density estimation shape constraint log-concave density Polya frequency function strongly unimodal iterative convex minorant algorithm active set algorithm


Walther, Guenther. Inference and Modeling with Log-concave Distributions. Statist. Sci. 24 (2009), no. 3, 319--327. doi:10.1214/09-STS303.

Export citation


  • An, M. Y. (1995). Log-concave probability distributions: Theory and statistical testing. Technical report, Economics Dept., Duke Univ.
  • An, M. Y. (1998). Logconcavity versus logconvexity: A complete characterization. J. Econom. Theory 80 350–369.
  • Avriel, M. (1972). r-convex functions. Math. Program. 2 309–323.
  • Bagnoli, M. and Bergstrom, T. (2005). Log-concave probability and its applications. Econometric Theory 26 445–469.
  • Balabdaoui, R., Rufibach, K. and Wellner, J. A. (2009). Limit distribution theory for maximum likelihood estimation of a log-concave density. Ann. Statist. 37 1299–1331.
  • Barlow, R. E. and Proschan, F. (1975). Statistical Theory of Reliability and Life Testing. Holt, Reinhart and Winston, New York.
  • Birgé, L. (1997). Estimation of unimodal densities without smoothness assumptions. Ann. Statist. 25 970–981.
  • Birgé, L. and Massart, P. (1993). Rates of convergence for minimum contrast estimators. Probab. Theory Related Fields 97 113–150.
  • Borell, C. (1975). Convex set functions in d-space. Period. Math. Hungar. 6 111–136.
  • Brooks, S. P. (1998). MCMC convergence diagnosis via multivariate bounds on log-concave densities. Ann. Statist. 26 398–433.
  • Caplin, A. and Nalebuff, B. (1991). Aggregation and social choice: A mean voter theorem. Econometrica 59 1–23.
  • Chang, G. and Walther, G. (2007). Clustering with mixtures of log-concave distributions. Comput. Statist. Data Anal. 51 6242–6251.
  • Cule, M., Samworth, R. and Stewart, M. (2008). Maximum likelihood estimation of a multidimensional log-concave density. Preprint.
  • Cule, M., Gramacy, R. and Samworth, R. (2009). LogConcDEAD: An R package for maximum likelihood estimation of a multivariate log-concave density. J. Statist. Software 29.
  • Dellaportas, P. and Smith, A. F. M. (1993). Bayesian inference for generalized linear and proportional hazards models via Gibbs sampling. J. Roy. Statist. Soc. Ser. C 42 443–460.
  • Dharmadhikari, S. and Joag-Dev, K. (1988). Unimodality, Convexity, and Applications. Academic Press, Boston.
  • 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. Preprint.
  • Dümbgen, L. and Rufibach, K. (2009). Maximum likelihood estimation of a log-concave density and its distribution function: Basic properties and uniform consistency. Bernoulli 15 40–68.
  • Eilers, P. H. C. (2005). Unimodal smoothing. J. Chemometrics 19 317–328.
  • Eilers, P. H. C. and Borgdorff, M. W. (2007). Non-parametric log-concave mixtures. Comput. Statist. Data Anal. 51 5444–5451.
  • Fletcher, R. (1987). Practical Methods for Optimization, 2nd ed. Wiley, New York.
  • Fraley, C. F. and Raftery, A. E. (2002). Model-based clustering, discriminant analysis, and density estimation. J. Amer. Statist. Assoc. 97 611–631.
  • Gilks, W. R. and Wild, P. (1992). Adaptive rejection sampling for Gibbs sampling. J. Roy. Statist. Soc. Ser. C 41 337–348.
  • Grenander, U. (1956). On the theory of mortality measurement, part II. Skandinavisk Aktuarietidskrift 39 125–153.
  • Groeneboom, P. and Wellner, J. A. (1992). Information Bounds and Nonparametric Maximum Likelihood Estimation. Birkhäuser, Basel.
  • Hampel, F. R. (1987). Design, modelling and analysis of some biological datasets. In Design, Data and Analysis, By Some Friends of Cuthbert Daniel (C. L. Mallows, ed.) 93–128. Wiley, New York.
  • Ibragimov, I. A. (1956). On the composition of unimodal distributions. Theory Probab. Appl. 1 255–260.
  • Jongbloed, G. (1998). The iterative convex minorant algorithm for nonparametric estimation. J. Comput. Graph. Statist. 7 310–321.
  • Kappel, F. and Kuntsevich, A. (2000). An implementation of Shor’s r-algorithm. Comput. Optim. Appl. 15 193–205.
  • Karlin, S. (1968). Total Positivity, Vol. 1. Stanford Univ. Press, Stanford.
  • Koenker, R. and Mizera, I. (2008). Quasi-concave density estimation. Technical report.
  • Marshall, A. W. and Olkin, I. (1979). Inequalities: Theory of Majorization and Its Applications. Academic Press, Orlando.
  • McLachlan, G. J. and Krishnan, T. (1997). The EM Algorithm and Extensions. Wiley, New York.
  • McLachlan, G. J. and Peel, D. (2000). Finite Mixture Models. Wiley, New York.
  • Müller, S. and Rufibach, K. (2009). Smooth tail index estimation. J. Statist. Comput. Simul. 79 1155–1167.
  • Pal, J., Woodroofe, M. and Meyer, M. (2007). Estimating a Polya frequency function. In Complex Datasets and Inverse Problems: Tomography, Networks and Beyond (R. Liu, W. Straderman, C.-H. Zhang, eds.) 239–249. IMS, Beachwood, OH.
  • Rufibach, K. (2006). Log-concave density estimation and bump hunting for I.I.D. observations. Ph.D. thesis, Univ. Bern and Göttingen.
  • Rufibach, K. (2007). Computing maximum likelihood estimators of a log-concave density function. J. Statist. Comput. Simul. 77 561–574.
  • Seregin, A. (2008). Estimation of convex-transformed densities. General Exam Proposal. Statistics Dept., Univ. Washington.
  • Silverman, B. W. (1982). On the estimation of a probability density function by the maximum penalized likelihood method. Ann. Statist. 10 795–810.
  • Street, W. M., Wolberg, W. H. and Mangasarian, O. L. (1993). Nuclear feature extraction for breast tumor diagnosis. IS& T/SPIE 1993 International Symposium on Electronic Imaging: Science and Technology, San Jose, CA, 1905, 861–870.
  • Terlaky, T. and Vial, J. P. (1998). Computing maximum likelihood estimators of convex density functions. SIAM J. Sci. Comput. 19 675–694.
  • Walther, G. (2001). Multiscale maximum likelihood analysis of a semiparametric model, with applications. Ann. Statist. 29 1298–1319.
  • Walther, G. (2002). Detecting the presence of mixing with multiscale maximum likelihood. J. Amer. Statist. Assoc. 97 508–513.
  • Wang, X., Woodroofe, M., Walker, M., Mateo, M. and Olszewski, E. (2005). Estimating dark matter distributions. Astrophys. J. 626 145–158.
  • Woodroofe, M. and Sun, J. (1993). A penalized maximum likelihood estimate of f(0+) when f is non-increasing. Statist. Sinica 3 501–515.