Maximum likelihood inference and bootstrap methods for plant organ growth via multi-phase kinetic models and their application to maize

* BACKGROUND AND AIMS: Fitting the parameters of models of plant organ growth is a means to investigate how environmental conditions affect plant architecture. The aim of this article is to evaluate some non-linear methods for fitting the parameters of multi-phase models of the kinetics of extension of plant organs such as laminae, sheaths and internodes. * METHODS: A set of computational procedures was developed allowing parameter-fitting of multi-phase models, using the maximum likelihood criterion, in which phases are identified with reference to ontogenic processes. Two bootstrap methods were compared to assess the precision of the estimates of fitted parameters, and of functions of these parameters such as the final leaf length, and the duration and rate of the rapid extension phase. Methods were applied to an experimental dataset, representing the kinetics of laminae, sheaths and internodes along the maize shoot, for two contrasting densities. * KEY RESULTS: A set of multi-phase models was proposed to describe the extension of laminae, sheaths and internodes along the shoot. The distinguishable phases differed between laminae, sheaths and internodes. For sheaths and laminae, but not for internodes, the same model could be fitted to all phytomers along the shoot. The variation of parameters along the shoot and between density treatments, as well as derived functions such as the durations of the phases of extension, are presented for laminae. It was the duration of the fast extension period, rather than its rate, which determined the difference in final length between treatments. * CONCLUSIONS: Such methods permit a large degree of objectivity and facilitate the analysis of such rather complicated but co-ordinated datasets. The work also illustrates some natural limitations of maximum likelihood methods, and viable ways of overcoming them by including a priori knowledge in the model fitting method are discussed.     

Diagnostic measures for empirical likelihood of general estimating equations

We develop diagnostic measures for assessing the influence ofindividual observations when using empirical likelihood withgeneral estimating equations, and we use these measures to constructgoodness-of-fit statistics for testing possible misspecificationin the estimating equations. Our diagnostics include case-deletionmeasures, local influence measures and pseudo-residuals. Ourgoodness-of-fit statistics include the sum of local influencemeasures and the processes of pseudo-residuals. Simulation studiesare conducted to evaluate our methods, and real datasets areanalyzed to illustrate the use of our diagnostic measures andgoodness-of-fit statistics.     

Estimation of the Size Distribution of Spherical Grains in the Boolean Model

A method of strong consistent estimation of the radius' distribution of spherical grains in the Boolean model in R^d, d ≦3, is proposed. It is based on the evaluation of the empirical capacity functional and application of the Glivenko-Cantelli theorem for random closed sets.     

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.     

Efficiencies of maximum likelihood methods of phylogenetic inferences when different substitution models are used

Choice of a substitution model is a crucial step in the maximum likelihood (ML) method of phylogenetic inference, and investigators tend to prefer complex mathematical models to simple ones. However, when complex models with many parameters are used, the extent of noise in statistical inferences increases, and thus complex models may not produce the true topology with a higher probability than simple ones. This problem was studied using computer simulation. When the number of nucleotides used was relatively large (1000 bp), the HKY+Gamma model showed smaller d(T) topological distance between the inferred and the true trees) than the JC and Kimura models. In the cases of shorter sequences (300 bp) simpler model and search algorithm such as JC model and SA+NNI search were found to be as efficient as more complicated searches and models in terms of topological distances, although the topologies obtained under HKY+Gamma model had the highest likelihood values. The performance of relatively simple search algorithm SA+NNI was found to be essentially the same as that of more extensive SA+TBR search under all models studied. Similarly to the conclusions reached by Takahashi and Nei [Mol. Biol. Evol. 17 (2000) 1251], our results indicate that simple models can be as efficient as complex models, and that use of complex models does not necessarily give more reliable trees compared with simple models.     

Conditional likelihood score functions for mixed models in linkage analysis

In this paper, we develop a general strategy for linkage analysis, applicable for arbitrary pedigree structures and genetic models with one major gene, polygenes and shared environmental effects. Extending work of Whittemore (1996), McPeek (1999) and Hossjer (2003d), the efficient score statistic is computed from a conditional likelihood of marker data given phenotypes. The resulting semiparametric linkage analysis is very similar to nonparametric linkage based on affected individuals. The efficient score S depends not only on identical-by-descent sharing and phenotypes, but also on a few parameters chosen by the user. We focus on (1) weak penetrance models, where the major gene has a small effect and (2) rare disease models, where the major gene has a possibly strong effect but the disease causing allele is rare. We illustrate our results for a large class of genetic models with a multivariate Gaussian liability. This class incorporates one major gene, polygenes and shared environmental effects in the liability, and allows e.g. binary, Gaussian, Poisson distributed and life-length phenotypes. A detailed simulation study is conducted for Gaussian phenotypes. The performance of the two optimal score functions S(wpairs) and S(normdom) are investigated. The conclusion is that (i) inclusion of polygenic effects into the score function increases overall performance for a wide range of genetic models and (ii) score functions based on the rare disease assumption are slightly more powerful.     

Flexible maximum likelihood methods for bivariate proportional hazards models

This article presents methodology for multivariate proportional hazards (PH) regression models. The methods employ flexible piecewise constant or spline specifications for baseline hazard functions in either marginal or conditional PH models, along with assumptions about the association among lifetimes. Because the models are parametric, ordinary maximum likelihood can be applied; it is able to deal easily with such data features as interval censoring or sequentially observed lifetimes, unlike existing semiparametric methods. A bivariate Clayton model (1978, Biometrika 65, 141-151) is used to illustrate the approach taken. Because a parametric assumption about association is made, efficiency and robustness comparisons are made between estimation based on the bivariate Clayton model and "working independence" methods that specify only marginal distributions for each lifetime variable.