10.6084/m9.figshare.3493526.v2
Pengfei Li
Pengfei
Li
Yukun Liu
Yukun
Liu
Jing Qin
Jing
Qin
Semiparametric Inference in a Genetic Mixture Model
Taylor & Francis Group
2017
Cluster mixture data
Degenerate Fisher information
Exponential tilting model
Genetic mixture model
2017-04-25 20:50:08
Journal contribution
https://tandf.figshare.com/articles/journal_contribution/Semiparametric_inference_in_a_genetic_mixture_model/3493526
<p>In genetic backcross studies, data are often collected from complex mixtures of distributions with known mixing proportions. Previous approaches to the inference of these genetic mixture models involve parameterizing the component distributions. However, model misspecification of any form is expected to have detrimental effects. We propose a semiparametric likelihood method for genetic mixture models: the empirical likelihood under the exponential tilting model assumption, in which the log ratio of the probability (density) functions from the components is linear in the observations. An application to mice cancer genetics involves random numbers of offspring within a litter. In other words, the cluster size is a random variable. We wish to test the null hypothesis that there is no difference between the two components in the mixture model, but unfortunately we find that the Fisher information is degenerate. As a consequence, the conventional two-term expansion in the likelihood ratio statistic does not work. By using a higher-order expansion, we are able to establish a nonstandard convergence rate <i>N</i><sup>− 1/4</sup> for the odds ratio parameter estimator <math><mi>β</mi><mo>^</mo></math>. Moreover, the limiting distribution of the empirical likelihood ratio statistic is derived. The underlying distribution function of each component can also be estimated semiparametrically. Analogously to the full parametric approach, we develop an expectation and maximization algorithm for finding the semiparametric maximum likelihood estimator. Simulation results and a real cancer application indicate that the proposed semiparametric method works much better than parametric methods. Supplementary materials for this article are available online.</p>