A note on profile likelihood for exponential tilt mixture models

Suppose that independent observations are drawn from multipledistributions, each of which is a mixture of two component distributionssuch that their log density ratio satisfies a linear model witha slope parameter and an intercept parameter. Inference forsuch models has been studied using empirical likelihood, andmixed results have been obtained. The profile empirical likelihoodof the slope and intercept has an irregularity at the null hypothesisso that the two component distributions are equal. We derivea profile empirical likelihood and maximum likelihood estimatorof the slope alone, and obtain the usual asymptotic propertiesfor the estimator and the likelihood ratio statistic regardlessof the null. Furthermore, we show the maximum likelihood estimatorof the slope and intercept jointly is consistent and asymptoticallynormal regardless of the null. At the null, the joint maximumlikelihood estimator falls along a straight line through theorigin with perfect correlation asymptotically to the firstorder.     

A note on semiparametric efficient inference for two-stage outcome-dependent sampling with a continuous outcome

Outcome-dependent sampling designs have been shown to be a cost-effectiveway to enhance study efficiency. We show that the outcome-dependentsampling design with a continuous outcome can be viewed as anextension of the two-stage case-control designs to the continuous-outcomecase. We further show that the two-stage outcome-dependent samplinghas a natural link with the missing-data and biased-samplingframeworks. Through the use of semiparametric inference andmissing-data techniques, we show that a certain semiparametricmaximum-likelihood estimator is computationally convenient andachieves the semiparametric efficient information bound. Wedemonstrate this both theoretically and through simulation.     

Probability distribution of molecular evolutionary trees: A new method of phylogenetic inference

A new method is presented for inferring evolutionary trees using nucleotide sequence data. The birth-death process is used as a model of speciation and extinction to specify the prior distribution of phylogenies and branching times. Nucleotide substitution is modeled by a continuous-time Markov process. Parameters of the branching model and the substitution model are estimated by maximum likelihood. The posterior probabilities of different phylogenies are calculated and the phylogeny with the highest posterior probability is chosen as the best estimate of the evolutionary relationship among species. We refer to this as the maximum posterior probability (MAP) tree. The posterior probability provides a natural measure of the reliability of the estimated phylogeny. Two example data sets are analyzed to infer the phylogenetic relationship of human, chimpanzee, gorilla, and orangutan. The best trees estimated by the new method are the same as those from the maximum likelihood analysis of separate topologies, but the posterior probabilities are quite different from the bootstrap proportions. The results of the method are found to be insensitive to changes in the rate parameter of the branching process. Correspondence to: Z. Yang     

A note on conditional AIC for linear mixed-effects models

The conventional model selection criterion, the Akaike informationcriterion, AIC, has been applied to choose candidate modelsin mixed-effects models by the consideration of marginal likelihood.Vaida & Blanchard (2005) demonstrated that such a marginalAIC and its small sample correction are inappropriate when theresearch focus is on clusters. Correspondingly, these authorssuggested the use of conditional AIC. Their conditional AICis derived under the assumption that the variance-covariancematrix or scaled variance-covariance matrix of random effectsis known. This note provides a general conditional AIC but withoutthese strong assumptions. Simulation studies show that the proposedmethod is promising.     

A Measure to Evaluate the Closeness of Satterthwaite's Approximation

Satterthwaite's approximation of the distribution of a nonnegative linear combination of independent mean squares is addressed in this article. A measure is developed to quantify the closeness of this approximation by using certain optimization results given by THIBAUDEAU and STYAN (1985). The main advantage of the proposed measure is to provide a theoretical framework for determining conditions under which Satterthwaite's approximation may be inadequate. This is demonstrated in three examples portraying frequently encountered problems in analysis of variance.