Markov and Semi‐Markov Switching Linear Mixed Models Used to Identify Forest Tree Growth Components

Tree growth is assumed to be mainly the result of three components: (i) an endogenous component assumed to be structured as a succession of roughly stationary phases separated by marked change points that are asynchronous among individuals, (ii) a time-varying environmental component assumed to take the form of synchronous fluctuations among individuals, and (iii) an individual component corresponding mainly to the local environment of each tree. To identify and characterize these three components, we propose to use semi-Markov switching linear mixed models, i.e., models that combine linear mixed models in a semi-Markovian manner. The underlying semi-Markov chain represents the succession of growth phases and their lengths (endogenous component) whereas the linear mixed models attached to each state of the underlying semi-Markov chain represent-in the corresponding growth phase-both the influence of time-varying climatic covariates (environmental component) as fixed effects, and interindividual heterogeneity (individual component) as random effects. In this article, we address the estimation of Markov and semi-Markov switching linear mixed models in a general framework. We propose a Monte Carlo expectation-maximization like algorithm whose iterations decompose into three steps: (i) sampling of state sequences given random effects, (ii) prediction of random effects given state sequences, and (iii) maximization. The proposed statistical modeling approach is illustrated by the analysis of successive annual shoots along Corsican pine trunks influenced by climatic covariates.

[1]  Paul Fearnhead,et al.  Computational methods for complex stochastic systems: a review of some alternatives to MCMC , 2008, Stat. Comput..

[2]  G. C. Wei,et al.  A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms , 1990 .

[3]  G. McLachlan,et al.  The EM Algorithm and Extensions: Second Edition , 2008 .

[4]  C. McCulloch.  Maximum Likelihood Variance Components Estimation for Binary Data , 1994 .

[5]  Galin L. Jones,et al.  Ascent‐based Monte Carlo expectation– maximization , 2005 .

[6]  S Y Lee,et al.  Latent variable models with mixed continuous and polytomous data , 2001, Biometrics.

[7]  J. Foulley,et al.  Heterogeneous variances in Gaussian linear mixed models , 1995, Genetics Selection Evolution.

[8]  W. Qian,et al.  Estimation of parameters in hidden Markov models , 1991, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[9]  G. McLachlan,et al.  The EM algorithm and extensions , 1996 .

[10]  Haikady N. Nagaraja,et al.  Inference in Hidden Markov Models , 2006, Technometrics.

[11]  Peter Bajorski,et al.  Wiley Series in Probability and Statistics , 2010 .

[12]  Padhraic Smyth,et al.  Segmental Hidden Markov Models with Random Effects for Waveform Modeling , 2006, J. Mach. Learn. Res..

[13]  Y. Guédon,et al.  Identifying ontogenetic, environmental and individual components of forest tree growth. , 2009, Annals of botany.

[14]  D. Harville,et al.  Some Bayesian and Non-Bayesian Procedures for the Analysis of Comparative Experiments and for Small-Area Estimation: Computational Aspects, Frequentist Properties, and Relationships , 1991 .

[15]  G. Celeux,et al.  Mixture of linear mixed models for clustering gene expression profiles from repeated microarray experiments , 2005 .

[16]  Y. Guédon.  Estimating Hidden Semi-Markov Chains From Discrete Sequences , 2003 .

[17]  S. Chib.  Calculating posterior distributions and modal estimates in Markov mixture models , 1996 .

[18]  J-B Durand,et al.  F F S S P P M M 0 0 4 4 Oral Presentations -session 2 Analysis of the Plant Architecture via Tree-structured Statistical Models: the Hidden Markov Trees , 2022 .

[19]  John M. Olin.  Calculating posterior distributions and modal estimates in Markov mixture models , 1996 .

[20]  Neri Merhav,et al.  Hidden Markov processes , 2002, IEEE Trans. Inf. Theory.

[21]  R. Altman.  Mixed Hidden Markov Models , 2007 .

[22]  C. McCulloch.  Maximum Likelihood Algorithms for Generalized Linear Mixed Models , 1997 .

[23]  Stephen Rapp,et al.  Qualitative longitudinal analysis of symptoms in patients with primary and metastatic brain tumours , 2008, Journal of the Royal Statistical Society. Series A,.

[24]  Yves Caraglio,et al.  Analyzing growth components in trees. , 2007, Journal of theoretical biology.

[25]  Sylvia Frühwirth-Schnatter,et al.  Finite Mixture and Markov Switching Models , 2006 .

[26]  Yann Guédon,et al.  Exploring the state sequence space for hidden Markov and semi-Markov chains , 2007, Comput. Stat. Data Anal..