Using Monte Carlo techniques to judge model prediction accuracy: Validation of the pesticide root zone model 3.12

Individuals from the Federal Insecticide, Fungicide, and Rodenticide Act (FIFRA) Environmental Model Validation Task Force (FEMVTF) Statistics Committee periodically met to discuss the mechanism for conducting an uncertainty analysis of Version 3.12 of the pesticide root zone model (PRZM 3.12) and to identify those model input parameters that most contribute to model prediction error. This activity was part of a larger project evaluating PRZM 3.12. The goal of the uncertainty analysis was to compare site-specific model predictions and field measurements using the variability in each as a basis of comparison. Monte Carlo analysis was used as an integral tool for judging the model's ability to predict accurately. The model was judged on how well it predicts measured values, taking into account the uncertainty in the model predictions. Monte Carlo analysis provides the tool for inferring model prediction uncertainty. We argue that this is a fairer test of the model than a simple one-to-one comparison between predictions and measurements. Because models are known to be imperfect predictors prior to running the model, the inaccuracy in model predictions should be considered when models are judged for their predictive ability. Otherwise, complex models can easily fail a validation test. Few complex models, such as PRZM 3.12, would pass a typical model validation exercise. This paper describes the approaches to the validation of PRZM 3.12 used by the committee and discusses issues in sampling distribution selection and appropriate statistics for interpreting the model validation results.

[1]  Piyush Singh,et al.  Sensitivity analysis for validating expert opinion as to ideal data set criteria for transport modeling , 2002, Environmental toxicology and chemistry.

[2]  J. Hammersley,et al.  A new Monte Carlo technique: antithetic variates , 1956, Mathematical Proceedings of the Cambridge Philosophical Society.

[3]  T. Kloek,et al.  Bayesian estimates of equation system parameters, An application of integration by Monte Carlo , 1976 .

[4]  D E Burmaster,et al.  Principles of good practice for the use of Monte Carlo techniques in human health and ecological risk assessments. , 1994, Risk analysis : an official publication of the Society for Risk Analysis.

[5]  M. Stein Large sample properties of simulations using latin hypercube sampling , 1987 .

[6]  W. K. Robertson,et al.  Effects of Plant‐Water Stress on Root Distribution of Corn, Soybeans, and Peanuts in Sandy Soil1 , 1980 .

[7]  W. Warren‐hicks,et al.  Development of a Monte Carlo sampling shell for the pesticide root zone model and its application by the Federal Insecticide, Fungicide, and Rodenticide Act Environmental modeling validation Task Force , 2002, Environmental Toxicology and Chemistry.

[8]  James R. Wilson Variance Reduction Techniques for Digital Simulation , 1984 .

[9]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[10]  L. Stone,et al.  Top and Root Growth of Irrigated and Nonirrigated Soybeans 1 , 1976 .

[11]  William Warren-Hicks,et al.  Monte Carlo analysis: Classical and Bayesian applications , 1996 .

[12]  M. D. McKay,et al.  Monte Carlo estimation under different distributions using the same simulation , 1987 .

[13]  John P Carbone,et al.  Validation of pesticide root zone model 3.12: Employing uncertainty analysis , 2002, Environmental toxicology and chemistry.

[14]  Boxin Tang Orthogonal Array-Based Latin Hypercubes , 1993 .