December 2022 Full Bayesian inference in hidden Markov models of plant growth
Gautier Viaud, Yuting Chen, Paul-Henry Cournède
Author Affiliations +
Ann. Appl. Stat. 16(4): 2352-2368 (December 2022). DOI: 10.1214/21-AOAS1594

Abstract

Accurately modeling the growth process of plants in interaction with their environment is important for predicting their biophysical characteristics, referred to as phenotype prediction. Most models are described by discrete dynamic systems in general state-space representation with important domain-specific characteristics: First, plant model parameters have usually clear functional meanings and may be of genetic origins, thus necessitating a precise estimation. Second, critical growth variables, specifically biomass production and dynamic allocation to organs, are hidden variables not accessible to measure. Finally, the difficulty to assess the local plant environment may imply the introduction of process noises in models. Therefore, a precise understanding of the system’s behavior requires the joint estimation of functional parameters, hidden states, and noise parameters. In this paper we describe how a full Bayesian method of estimation can accurately estimate all these key model variables using Markov chain Monte Carlo (MCMC) techniques. In the presence of both process and observation noises, it requires to use adequate particle MCMC (PMCMC) algorithms to efficiently sample the hidden states which, consequently, allows for a precise estimation of all noise parameters involved. Thanks to the Bayesian framework, appropriate choices of prior distributions for the noise parameters have enabled analytical posterior distributions and only simple updates are required. Furthermore, this estimation strategy can be easily generalized and adapted to different types of plant growth models, such as organ-scale or compartmental, provided that they are formulated as hidden Markov models. Our estimation method improves on those classically used in plant growth modeling in several aspects: First, by building upon a general probabilistic framework the estimation results allow proper statistical analyses. It is useful in prediction, no only for uncertainty and risk analysis (e.g., for crop yield prediction) but also to analyze the results of experimental trials, for example, to compare genotypes in breeding. Moreover, the care taken in the estimation of hidden variables opens new perspectives in the understanding of inner growth processes, notably the balance and interaction between biomass production and allocation (referred to as source-sink dynamics). Applications of this estimation procedure are demonstrated on the GreenLab model for Arabidopsis thaliana and the Log-Normal Allocation and Senescence (LNAS) model for sugar beet, on both synthetic and real data.

Citation

Download Citation

Gautier Viaud. Yuting Chen. Paul-Henry Cournède. "Full Bayesian inference in hidden Markov models of plant growth." Ann. Appl. Stat. 16 (4) 2352 - 2368, December 2022. https://doi.org/10.1214/21-AOAS1594

Information

Received: 1 December 2020; Revised: 1 December 2021; Published: December 2022
First available in Project Euclid: 26 September 2022

MathSciNet: MR4489214
zbMATH: 1498.62260
Digital Object Identifier: 10.1214/21-AOAS1594

Keywords: Bayesian inference , Hidden Markov model , particle MCMC , plant growth

Rights: Copyright © 2022 Institute of Mathematical Statistics

Vol.16 • No. 4 • December 2022
Back to Top