Combined use of local and ANOVA-based global sensitivity analyses for the investigation of a stochastic dynamic model: Application to the case study of an individual-based model of a fish population

被引:65
作者
Ginot, V
Gaba, S
Beaudouin, R
Aries, F
Monod, H
机构
[1] INRA, Unite Biometrie, F-84814 Avignon 9, France
[2] INRA, SCRIBE, F-35042 Rennes, France
[3] INRA, Unite Math & Informat Appl, F-78352 Jouy En Josas, France
关键词
sensitivity analysis; ANOVA; simulation design; computer experiments; stochastic individual-based model;
D O I
10.1016/j.ecolmodel.2005.08.025
中图分类号
Q14 [生态学(生物生态学)];
学科分类号
071012 ; 0713 ;
摘要
Global methods based on variance decomposition are increasingly being used for sensitivity analysis (SA). Of these, analysis of variance (ANOVA) is surprisingly rarely employed. Yet, it is a viable alternative to other model-free methods, as it gives comparable results and it readily available in most statistical packages. Furthermore, decomposing the input factors of ANOVA into orthogonal polynomial effects gives additional insights into the way a parameter impacts on model output (linear, quadratic and cubic). However, using global methods should not lead modellers to forego local methods, which provide additional information, as, for example, time course analysis of local sensitivity coefficients. We illustrate the use of these techniques, particularly ANOVA, on a stochastic individual-based model of a mosquitofish (Gambusia holbrooki) population in experimental tanks. Local SA led to unexpected and counter-intuitive results, indicating that the model output (population size) was much more sensitive to the fecundity threshold (length at first parturition) than to the fecundity parameter (brood size). Time course analysis of local coefficients suggested that, as far as calibration is concerned, it would probably be impossible to determine more than two parameters on the sole records of the population size in time. Global SA (ANOVA) was targeted to assess which processes had an impact on the model outcome in our experimental conditions, by exploring the parameter space over the entire biological range of all parameters. It showed that parameters had mainly linear and additive effects (few interactions) on the output in a logarithmic scale, and that the main processes involved in population growth were individual growth and adult survival, followed by the breeding process. Juvenile survival had a lesser impact. (C) 2005 Elsevier B.V. All rights reserved.
引用
收藏
页码:479 / 491
页数:13
相关论文
共 30 条
[1]  
Atkinson A.C., 1992, OPTIMUM EXPT DESIGNS
[2]  
Banks J., 1998, HDB SIMULATION PRINC
[3]  
Bard Y., 1974, Nonlinear Parameter Estimation
[4]  
Box G, 1987, EMPIRICAL MODEL BUIL
[5]  
Box G. E. P., 1964, JR STAT SOC B, V26, P211
[6]   Simplifying Sirius: sensitivity analysis and development of a meta-model for wheat yield prediction [J].
Brooks, RJ ;
Semenov, MA ;
Jamieson, PD .
EUROPEAN JOURNAL OF AGRONOMY, 2001, 14 (01) :43-60
[7]  
Campolongo F., 2000, SENSITIVITY ANAL, P51
[8]  
Christensen R., 1987, PLANE ANSWERS COMPLE
[9]   Understanding 3-PG using a sensitivity analysis [J].
Esprey, LJ ;
Sands, PJ ;
Smith, CW .
FOREST ECOLOGY AND MANAGEMENT, 2004, 193 (1-2) :235-250
[10]   Assessing uncertainty in ecological systems using global sensitivity analyses: a case example of simulated wolf reintroduction effects on elk [J].
Fieberg, J ;
Jenkins, KJ .
ECOLOGICAL MODELLING, 2005, 187 (2-3) :259-280