Abstract
In this article, we revisit the problem of fitting a mixture model under the assumption that the mixture components are symmetric and log-concave. To this end, we first study the nonparametric maximum likelihood estimation (MLE) of a monotone log-concave probability density. To fit the mixture model, we propose a semiparametric EM (SEM) algorithm, which can be adapted to other semiparametric mixture models. In our numerical experiments, we compare our algorithm to that of Balabdaoui and Doss (Citation2018, Inference for a two-component mixture of symmetric distributions under log-concavity. Bernoulli 24 (2):1053–71) and other mixture models both on simulated and real-world datasets.
Acknowledgement
We are grateful to Kaspar Rufibach and Lutz Dümbgen for their feedback on the adaptation of their work in Section 2. We are also grateful to Günther Walther for answering questions of identifiability and for calling our attention to the work of Balabdaoui and Butucea (Citation2014).
Notes
1 The method is based on an active set implementation and is available in the R package logcondens.mode.
2 Following the definition in Rockafellar (Citation2015), a concave function f is said to be proper if for at least one x and for every x. A closed function is a function that maps closed sets to closed sets.