A stochastic model related to the Richards-type growth curve. Estimation by means of simulated annealing and variable neighborhood search

For a Richards-type curve a diffusion process is constructed.The distribution and main characteristics of the process are analyzed.The estimation of the parameters is addressed by maximum likelihood.A procedure for bounding the parametric space is established.VND and SA algorithms are used for estimating the parameters. A stochastic diffusion model related to a reformulation of the Richards growth curve is proposed. The main characteristics of the process are studied, and the problem of maximum likelihood estimation for the parameters of the process is considered. Since a complex system of equations appears whose solution cannot be guaranteed via the classic numerical procedures, we suggest the use of metaheuristic optimization algorithms such as simulated annealing and variable neighborhood search. Given that the space of solutions is continuous and unbounded, some strategies are suggested for bounding it, and a description is provided for the application of the selected algorithms. In the case of the variable neighborhood search algorithm, a hybrid method is proposed in which it is combined with simulated annealing. Some examples based on simulated sample paths are developed in order to test the validity of the bounding method for the space of solutions, and a comparison is made between the application of both methods.

[1]  M. Omizo,et al.  Modeling , 1983, Encyclopedic Dictionary of Archaeology.

[2]  J. Pitchford,et al.  Stochastic von Bertalanffy models, with applications to fish recruitment. , 2007, Journal of theoretical biology.

[3]  Pierre Hansen,et al.  Variable Neighborhood Search , 2018, Handbook of Heuristics.

[4]  F. J. Richards.  A Flexible Growth Function for Empirical Use , 1959 .

[5]  Alexandre Souto Martinez,et al.  Continuous growth models in terms of generalized logarithm and exponential functions , 2008, 0803.2635.

[6]  H. Schurz.  MODELING, ANALYSIS AND DISCRETIZATION OF STOCHASTIC LOGISTIC EQUATIONS , 2007 .

[7]  Estimation of the Gentan probability using the price of logs in private plantation forests , 2007, Journal of Forest Research.

[8]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[9]  F. Busetti.  Simulated annealing overview , 2000 .

[10]  V. Giorno,et al.  Inferring the effect of therapy on tumors showing stochastic Gompertzian growth. , 2011, Journal of theoretical biology.

[11]  Seyed Taghi Akhavan Niaki,et al.  Estimating the four parameters of the Burr III distribution using a hybrid method of variable neighborhood search and iterated local search algorithms , 2012, Appl. Math. Comput..

[12]  J. Santos,et al.  A growth model for secondary forest in Central Amazonia , 2005 .

[13]  J. Brida,et al.  Re-formulation of the Solow economic growth model whit the Richards population growth law , 2005 .

[14]  J. Ludwig,et al.  grofit: Fitting Biological Growth Curves with R , 2010 .

[15]  Patrick Siarry,et al.  A survey on optimization metaheuristics , 2013, Inf. Sci..

[16]  Christos H. Skiadas,et al.  A Stochastic Logistic Innovation Diffusion Model Studying the Electricity Consumption in Greece and the United States , 1999 .

[17]  R. G. Jáimez,et al.  A note on the Volterra integral equation for the first-passage-time probability density , 1995, Journal of Applied Probability.

[18]  P. Román-Román,et al.  A diffusion process to model generalized von Bertalanffy growth patterns: fitting to real data. , 2010, Journal of theoretical biology.

[19]  H. Dette,et al.  Efficient experimental designs for sigmoidal growth models , 2005 .

[20]  H. Simianer,et al.  Modeling the growth of the Goettingen minipig. , 2007, Journal of animal science.

[21]  R. Fletcher,et al.  A general solution for the complete Richards function , 1975 .

[22]  P. Lánský,et al.  Diffusion approximation and first-passage-time problem for a model neuron III. Abirth-and-death process approach , 1988 .

[23]  Jianhong Wu,et al.  Richards model revisited: validation by and application to infection dynamics. , 2012, Journal of theoretical biology.

[24]  Patricia Román-Román,et al.  Estimating the parameters of a Gompertz-type diffusion process by means of Simulated Annealing , 2012, Appl. Math. Comput..

[25]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[26]  Nuria Rico,et al.  APPROXIMATE AND GENERALIZED CONFIDENCE BANDS FOR SOME PARAMETRIC FUNCTIONS OF THE LOGNORMAL DIFFUSION PROCESS WITH EXOGENOUS FACTORS , 2006 .

[27]  Patricia Román-Román,et al.  FORECASTING FOR THE UNIVARIATE LOGNORMAL DIFFUSION PROCESS WITH EXOGENOUS FACTORS , 2003, Cybern. Syst..

[28]  M. Artzrouni,et al.  Stochastic differential equations in mathematical demography: a review , 1990 .

[29]  E. Tjørve,et al.  A unified approach to the Richards-model family for use in growth analyses: why we need only two model forms. , 2010, Journal of theoretical biology.

[30]  L. Pienaar,et al.  The Chapman-Richards Generalization of Von Bertalanffy's Growth Model for Basal Area Growth and Yield in Even - Aged Stands , 1973 .

[31]  Jianyang Xia,et al.  Independent effects of warming and nitrogen addition on plant phenology in the Inner Mongolian steppe. , 2013, Annals of botany.

[32]  Bruce D. McCullough,et al.  Statistical Analysis of Stochastic Processes in Time , 2005, Technometrics.

[33]  D. R. Causton.  A Computer Program for Fitting the Richards Function , 1969 .

[34]  Zeljko Bajzer,et al.  Mathematical modeling of cancer radiovirotherapy. , 2006, Mathematical biosciences.

[35]  Seyed Taghi Akhavan Niaki,et al.  A hybrid variable neighborhood search and simulated annealing algorithm to estimate the three parameters of the Weibull distribution , 2011, Expert Syst. Appl..

[36]  S. Aggrey,et al.  Growth characteristics of pearl gray guinea fowl as predicted by the Richards, Gompertz, and logistic models. , 2006, Poultry science.

[37]  R. Gutiérrez-Jáimez,et al.  A new Gompertz-type diffusion process with application to random growth. , 2007, Mathematical biosciences.

[38]  Virginia Giorno,et al.  On the effect of a therapy able to modify both the growth rates in a Gompertz stochastic model. , 2013, Mathematical biosciences.

[39]  Paolo Baldi,et al.  Lévy processes and stochastic von Bertalanffy models of growth, with application to fish population analysis. , 2009, Journal of theoretical biology.

[40]  Forecasting Fruit Size and Caliber by Means of Diffusion Processes. Application to “Valencia Late” Oranges , 2014 .

[41]  A reformulation of the Ramsey model , 2006 .

[42]  Shunsuke Sato,et al.  Diffusion approximation and first passage time problem for a model neuron. II. Outline of a computation method , 1983 .

[43]  José A. Díaz-García,et al.  A global simulated annealing heuristic for the three-parameter lognormal maximum likelihood estimation , 2008, Computational Statistics & Data Analysis.

[44]  R. Capocelli,et al.  A diffusion model for population growth in random environment. , 1974, Theoretical population biology.

[45]  Growth with regulation in fluctuating environments , 1984, Biological cybernetics.

[46]  J. D. Val,et al.  Inference for the Richards growth model using Box and Cox transformation and bootstrap techniques , 2006 .

[47]  Patricia Román-Román,et al.  Modelling logistic growth by a new diffusion process: Application to biological systems , 2012, Biosyst..

[48]  Virginia Giorno,et al.  AN OUTLINE OF THEORETICAL AND ALGORITHMIC APPROACHES TO FIRST PASSAGE TIME PROBLEMS WITH APPLICATIONS TO BIOLOGICAL MODELING , 1999 .