Sensitivity analysis
Sensitivity analysis is the study of how the uncertainty in the output of a mathematical model or system (numerical or otherwise) can be divided and allocated to different sources of uncertainty in its inputs.[1][2] A related practice is uncertainty analysis, which has a greater focus on uncertainty quantification and propagation of uncertainty; ideally, uncertainty and sensitivity analysis should be run in tandem.
The process of recalculating outcomes under alternative assumptions to determine the impact of a variable under sensitivity analysis can be useful for a range of purposes,[3] including:
- Testing the robustness of the results of a model or system in the presence of uncertainty.
- Increased understanding of the relationships between input and output variables in a system or model.
- Uncertainty reduction, through the identification of model input that cause significant uncertainty in the output and should therefore be the focus of attention in order to increase robustness (perhaps by further research).
- Searching for errors in the model (by encountering unexpected relationships between inputs and outputs).
- Model simplification – fixing model input that has no effect on the output, or identifying and removing redundant parts of the model structure.
- Enhancing communication from modelers to decision makers (e.g. by making recommendations more credible, understandable, compelling or persuasive).
- Finding regions in the space of input factors for which the model output is either maximum or minimum or meets some optimum criterion (see optimization and Monte Carlo filtering).
- In case of calibrating models with large number of parameters, a primary sensitivity test can ease the calibration stage by focusing on the sensitive parameters. Not knowing the sensitivity of parameters can result in time being uselessly spent on non-sensitive ones.[4]
- To seek to identify important connections between observations, model inputs, and predictions or forecasts, leading to the development of better models.[5][6]
Overview
A mathematical model (for example in biology, climate change, economics or engineering) can be highly complex, and as a result, its relationships between inputs and outputs may be poorly understood. In such cases, the model can be viewed as a black box, i.e. the output is an "opaque" function of its inputs.
Quite often, some or all of the model inputs are subject to sources of uncertainty, including errors of measurement, absence of information and poor or partial understanding of the driving forces and mechanisms. This uncertainty imposes a limit on our confidence in the response or output of the model. Further, models may have to cope with the natural intrinsic variability of the system (aleatory), such as the occurrence of stochastic events.[7]
Good modeling practice requires that the modeler provide an evaluation of the confidence in the model. This requires, first, a quantification of the uncertainty in any model results (uncertainty analysis); and second, an evaluation of how much each input is contributing to the output uncertainty. Sensitivity analysis addresses the second of these issues (although uncertainty analysis is usually a necessary precursor), performing the role of ordering by importance the strength and relevance of the inputs in determining the variation in the output.[2]
In models involving many input variables, sensitivity analysis is an essential ingredient of model building and quality assurance. National and international agencies involved in impact assessment studies have included sections devoted to sensitivity analysis in their guidelines. Examples are the European Commission (see e.g. the guidelines for impact assessment),[8] the White House Office of Management and Budget, the Intergovernmental Panel on Climate Change and US Environmental Protection Agency's modeling guidelines.[9] In a comment published in 2020 in the journal Nature 22 scholars take COVID-19 as the occasion for suggesting five ways to make models serve society better. One of the five recommendations, under the heading of 'Mind the assumptions' is to 'perform global uncertainty and sensitivity analyses [...] allowing all that is uncertain — variables, mathematical relationships and boundary conditions — to vary simultaneously as runs of the model produce its range of predictions.'[10]
Settings, constraints, and related issues
Settings and constraints
The choice of method of sensitivity analysis is typically dictated by a number of problem constraints or settings. Some of the most common are
- Computational expense: Sensitivity analysis is almost always performed by running the model a (possibly large) number of times, i.e. a sampling-based approach.[11] This can be a significant problem when,
- A single run of the model takes a significant amount of time (minutes, hours or longer). This is not unusual with very complex models.
- The model has a large number of uncertain inputs. Sensitivity analysis is essentially the exploration of the multidimensional input space, which grows exponentially in size with the number of inputs. See the curse of dimensionality.
- Computational expense is a problem in many practical sensitivity analyses. Some methods of reducing computational expense include the use of emulators (for large models), and screening methods (for reducing the dimensionality of the problem). Another method is to use an event-based sensitivity analysis method for variable selection for time-constrained applications.[12] This is an input variable selection (IVS) method that assembles together information about the trace of the changes in system inputs and outputs using sensitivity analysis to produce an input/output trigger/event matrix that is designed to map the relationships between input data as causes that trigger events and the output data that describes the actual events. The cause-effect relationship between the causes of state change i.e. input variables and the effect system output parameters determines which set of inputs have a genuine impact on a given output. The method has a clear advantage over analytical and computational IVS method since it tries to understand and interpret system state change in the shortest possible time with minimum computational overhead.[12][13]
- Correlated inputs: Most common sensitivity analysis methods assume independence between model inputs, but sometimes inputs can be strongly correlated. This is still an immature field of research and definitive methods have yet to be established.
- Nonlinearity: Some sensitivity analysis approaches, such as those based on linear regression, can inaccurately measure sensitivity when the model response is nonlinear with respect to its inputs. In such cases, variance-based measures are more appropriate.
- Model interactions: Interactions occur when the perturbation of two or more inputs simultaneously causes variation in the output greater than that of varying each of the inputs alone. Such interactions are present in any model that is non-additive, but will be neglected by methods such as scatterplots and one-at-a-time perturbations.[14] The effect of interactions can be measured by the total-order sensitivity index.
- Multiple outputs: Virtually all sensitivity analysis methods consider a single univariate model output, yet many models output a large number of possibly spatially or time-dependent data. Note that this does not preclude the possibility of performing different sensitivity analyses for each output of interest. However, for models in which the outputs are correlated, the sensitivity measures can be hard to interpret.
- Given data: While in many cases the practitioner has access to the model, in some instances a sensitivity analysis must be performed with "given data", i.e. where the sample points (the values of the model inputs for each run) cannot be chosen by the analyst. This may occur when a sensitivity analysis has to be performed retrospectively, perhaps using data from an optimisation or uncertainty analysis, or when data comes from a discrete source.[15]
Assumptions vs. inferences
In uncertainty and sensitivity analysis there is a crucial trade off between how scrupulous an analyst is in exploring the input assumptions and how wide the resulting inference may be. The point is well illustrated by the econometrician Edward E. Leamer:[16][17]
I have proposed a form of organized sensitivity analysis that I call 'global sensitivity analysis' in which a neighborhood of alternative assumptions is selected and the corresponding interval of inferences is identified. Conclusions are judged to be sturdy only if the neighborhood of assumptions is wide enough to be credible and the corresponding interval of inferences is narrow enough to be useful.
Note Leamer's emphasis is on the need for 'credibility' in the selection of assumptions. The easiest way to invalidate a model is to demonstrate that it is fragile with respect to the uncertainty in the assumptions or to show that its assumptions have not been taken 'wide enough'. The same concept is expressed by Jerome R. Ravetz, for whom bad modeling is when uncertainties in inputs must be suppressed lest outputs become indeterminate.[18]
Pitfalls and difficulties
Some common difficulties in sensitivity analysis include
- Too many model inputs to analyse. Screening can be used to reduce dimensionality. Another way to tackle the curse of dimensionality is to use sampling based on low discrepancy sequences[19]
- The model takes too long to run. Emulators (including HDMR) can reduce the number of model runs needed.
- There is not enough information to build probability distributions for the inputs. Probability distributions can be constructed from expert elicitation, although even then it may be hard to build distributions with great confidence. The subjectivity of the probability distributions or ranges will strongly affect the sensitivity analysis.
- Unclear purpose of the analysis. Different statistical tests and measures are applied to the problem and different factors rankings are obtained. The test should instead be tailored to the purpose of the analysis, e.g. one uses Monte Carlo filtering if one is interested in which factors are most responsible for generating high/low values of the output.
- Too many model outputs are considered. This may be acceptable for the quality assurance of sub-models but should be avoided when presenting the results of the overall analysis.
- Piecewise sensitivity. This is when one performs sensitivity analysis on one sub-model at a time. This approach is non conservative as it might overlook interactions among factors in different sub-models (Type II error).
- Commonly used OAT approach is not valid for nonlinear models. Global methods should be used instead.[20]
Sensitivity analysis methods
There are a large number of approaches to performing a sensitivity analysis, many of which have been developed to address one or more of the constraints discussed above.[2] They are also distinguished by the type of sensitivity measure, be it based on (for example) variance decompositions, partial derivatives or elementary effects. In general, however, most procedures adhere to the following outline:
- Quantify the uncertainty in each input (e.g. ranges, probability distributions). Note that this can be difficult and many methods exist to elicit uncertainty distributions from subjective data.[21]
- Identify the model output to be analysed (the target of interest should ideally have a direct relation to the problem tackled by the model).
- Run the model a number of times using some design of experiments,[22] dictated by the method of choice and the input uncertainty.
- Using the resulting model outputs, calculate the sensitivity measures of interest.
In some cases this procedure will be repeated, for example in high-dimensional problems where the user has to screen out unimportant variables before performing a full sensitivity analysis.
The various types of "core methods" (discussed below) are distinguished by the various sensitivity measures which are calculated. These categories can somehow overlap. Alternative ways of obtaining these measures, under the constraints of the problem, can be given.
One-at-a-time (OAT)
One of the simplest and most common approaches is that of changing one-factor-at-a-time (OAT), to see what effect this produces on the output.[23][24][25] OAT customarily involves
- moving one input variable, keeping others at their baseline (nominal) values, then,
- returning the variable to its nominal value, then repeating for each of the other inputs in the same way.
Sensitivity may then be measured by monitoring changes in the output, e.g. by partial derivatives or linear regression. This appears a logical approach as any change observed in the output will unambiguously be due to the single variable changed. Furthermore, by changing one variable at a time, one can keep all other variables fixed to their central or baseline values. This increases the comparability of the results (all 'effects' are computed with reference to the same central point in space) and minimizes the chances of computer programme crashes, more likely when several input factors are changed simultaneously. OAT is frequently preferred by modellers because of practical reasons. In case of model failure under OAT analysis the modeller immediately knows which is the input factor responsible for the failure.[14]
Despite its simplicity however, this approach does not fully explore the input space, since it does not take into account the simultaneous variation of input variables. This means that the OAT approach cannot detect the presence of interactions between input variables.[26]
Derivative-based local methods
Local derivative-based methods involve taking the partial derivative of the output Y with respect to an input factor Xi :
where the subscript x0 indicates that the derivative is taken at some fixed point in the space of the input (hence the 'local' in the name of the class). Adjoint modelling[27][28] and Automated Differentiation[29] are methods in this class. Similar to OAT, local methods do not attempt to fully explore the input space, since they examine small perturbations, typically one variable at a time. It is possible to select similar samples from derivative-based sensitivity through Neural Networks and perform uncertainty quantification. [30]
Regression analysis
Regression analysis, in the context of sensitivity analysis, involves fitting a linear regression to the model response and using standardized regression coefficients as direct measures of sensitivity. The regression is required to be linear with respect to the data (i.e. a hyperplane, hence with no quadratic terms, etc., as regressors) because otherwise it is difficult to interpret the standardised coefficients. This method is therefore most suitable when the model response is in fact linear; linearity can be confirmed, for instance, if the coefficient of determination is large. The advantages of regression analysis are that it is simple and has a low computational cost.
Variance-based methods
Variance-based methods[31][32][33] are a class of probabilistic approaches which quantify the input and output uncertainties as probability distributions, and decompose the output variance into parts attributable to input variables and combinations of variables. The sensitivity of the output to an input variable is therefore measured by the amount of variance in the output caused by that input. These can be expressed as conditional expectations, i.e., considering a model Y = f(X) for X = {X1, X2, ... Xk}, a measure of sensitivity of the ith variable Xi is given as,
where "Var" and "E" denote the variance and expected value operators respectively, and X~i denotes the set of all input variables except Xi. This expression essentially measures the contribution Xi alone to the uncertainty (variance) in Y (averaged over variations in other variables), and is known as the first-order sensitivity index or main effect index. Importantly, it does not measure the uncertainty caused by interactions with other variables. A further measure, known as the total effect index, gives the total variance in Y caused by Xi and its interactions with any of the other input variables. Both quantities are typically standardised by dividing by Var(Y).
Variance-based methods allow full exploration of the input space, accounting for interactions, and nonlinear responses. For these reasons they are widely used when it is feasible to calculate them. Typically this calculation involves the use of Monte Carlo methods, but since this can involve many thousands of model runs, other methods (such as emulators) can be used to reduce computational expense when necessary. Note that full variance decompositions are only meaningful when the input factors are independent from one another.[34]
Variogram analysis of response surfaces (VARS)
One of the major shortcomings of the previous sensitivity analysis methods is that none of them considers the spatially ordered structure of the response surface/output of the model Y=f(X) in the parameter space. By utilizing the concepts of directional variograms and covariograms, variogram analysis of response surfaces (VARS) addresses this weakness through recognizing a spatially continuous correlation structure to the values of Y, and hence also to the values of .[35][36]
Basically, the higher the variability the more heterogeneous is the response surface along a particular direction/parameter, at a specific perturbation scale. Accordingly, in the VARS framework, the values of directional variograms for a given perturbation scale can be considered as a comprehensive illustration of sensitivity information, through linking variogram analysis to both direction and perturbation scale concepts. As a result, the VARS framework accounts for the fact that sensitivity is a scale-dependent concept, and thus overcomes the scale issue of traditional sensitivity analysis methods.[37] More importantly, VARS is able to provide relatively stable and statistically robust estimates of parameter sensitivity with much lower computational cost than other strategies (about two orders of magnitude more efficient).[38] Noteworthy, it has been shown that there is a theoretical link between the VARS framework and the variance-based and derivative-based approaches.
Screening
Screening is a particular instance of a sampling-based method. The objective here is rather to identify which input variables are contributing significantly to the output uncertainty in high-dimensionality models, rather than exactly quantifying sensitivity (i.e. in terms of variance). Screening tends to have a relatively low computational cost when compared to other approaches, and can be used in a preliminary analysis to weed out uninfluential variables before applying a more informative analysis to the remaining set. One of the most commonly used screening method is the elementary effect method.[39][40]
Scatter plots
A simple but useful tool is to plot scatter plots of the output variable against individual input variables, after (randomly) sampling the model over its input distributions. The advantage of this approach is that it can also deal with "given data", i.e., a set of arbitrarily-placed data points, and gives a direct visual indication of sensitivity. Quantitative measures can also be drawn, for example by measuring the correlation between Y and Xi, or even by estimating variance-based measures by nonlinear regression.[15]
Alternative methods
A number of methods have been developed to overcome some of the constraints discussed above, which would otherwise make the estimation of sensitivity measures infeasible (most often due to computational expense). Generally, these methods focus on efficiently calculating variance-based measures of sensitivity.
Emulators
Emulators (also known as metamodels, surrogate models or response surfaces) are data-modeling/machine learning approaches that involve building a relatively simple mathematical function, known as an emulator, that approximates the input/output behavior of the model itself.[41] In other words, it is the concept of "modeling a model" (hence the name "metamodel"). The idea is that, although computer models may be a very complex series of equations that can take a long time to solve, they can always be regarded as a function of their inputs Y = f(X). By running the model at a number of points in the input space, it may be possible to fit a much simpler emulator η(X), such that η(X) ≈ f(X) to within an acceptable margin of error.[42] Then, sensitivity measures can be calculated from the emulator (either with Monte Carlo or analytically), which will have a negligible additional computational cost. Importantly, the number of model runs required to fit the emulator can be orders of magnitude less than the number of runs required to directly estimate the sensitivity measures from the model.[43]
Clearly, the crux of an emulator approach is to find an η (emulator) that is a sufficiently close approximation to the model f. This requires the following steps,
- Sampling (running) the model at a number of points in its input space. This requires a sample design.
- Selecting a type of emulator (mathematical function) to use.
- "Training" the emulator using the sample data from the model – this generally involves adjusting the emulator parameters until the emulator mimics the true model as well as possible.
Sampling the model can often be done with low-discrepancy sequences, such as the Sobol sequence – due to mathematician Ilya M. Sobol or Latin hypercube sampling, although random designs can also be used, at the loss of some efficiency. The selection of the emulator type and the training are intrinsically linked since the training method will be dependent on the class of emulator. Some types of emulators that have been used successfully for sensitivity analysis include,
- Gaussian processes[43] (also known as kriging), where any combination of output points is assumed to be distributed as a multivariate Gaussian distribution. Recently, "treed" Gaussian processes have been used to deal with heteroscedastic and discontinuous responses.[44][45]
- Random forests,[41] in which a large number of decision trees are trained, and the result averaged.
- Gradient boosting,[41] where a succession of simple regressions are used to weight data points to sequentially reduce error.
- Polynomial chaos expansions,[46] which use orthogonal polynomials to approximate the response surface.
- Smoothing splines,[47] normally used in conjunction with HDMR truncations (see below).
The use of an emulator introduces a machine learning problem, which can be difficult if the response of the model is highly nonlinear. In all cases, it is useful to check the accuracy of the emulator, for example using cross-validation.
High-dimensional model representations (HDMR)
A high-dimensional model representation (HDMR)[48][49] (the term is due to H. Rabitz[50]) is essentially an emulator approach, which involves decomposing the function output into a linear combination of input terms and interactions of increasing dimensionality. The HDMR approach exploits the fact that the model can usually be well-approximated by neglecting higher-order interactions (second or third-order and above). The terms in the truncated series can then each be approximated by e.g. polynomials or splines (REFS) and the response expressed as the sum of the main effects and interactions up to the truncation order. From this perspective, HDMRs can be seen as emulators which neglect high-order interactions; the advantage is that they are able to emulate models with higher dimensionality than full-order emulators.
Fourier amplitude sensitivity test (FAST)
The Fourier amplitude sensitivity test (FAST) uses the Fourier series to represent a multivariate function (the model) in the frequency domain, using a single frequency variable. Therefore, the integrals required to calculate sensitivity indices become univariate, resulting in computational savings.
Applications
Examples of sensitivity analyses can be found in various area of application, such as:
- Environmental sciences
- Business
- Social sciences
- Chemistry
- Engineering
- Epidemiology
- Meta-analysis
- Multi-criteria decision making
- Time-critical decision making
- Model calibration
- Uncertainty Quantification
Sensitivity auditing
It may happen that a sensitivity analysis of a model-based study is meant to underpin an inference, and to certify its robustness, in a context where the inference feeds into a policy or decision making process. In these cases the framing of the analysis itself, its institutional context, and the motivations of its author may become a matter of great importance, and a pure sensitivity analysis – with its emphasis on parametric uncertainty – may be seen as insufficient. The emphasis on the framing may derive inter-alia from the relevance of the policy study to different constituencies that are characterized by different norms and values, and hence by a different story about 'what the problem is' and foremost about 'who is telling the story'. Most often the framing includes more or less implicit assumptions, which could be political (e.g. which group needs to be protected) all the way to technical (e.g. which variable can be treated as a constant).
In order to take these concerns into due consideration the instruments of SA have been extended to provide an assessment of the entire knowledge and model generating process. This approach has been called 'sensitivity auditing'. It takes inspiration from NUSAP,[53] a method used to qualify the worth of quantitative information with the generation of `Pedigrees' of numbers. Likewise, sensitivity auditing has been developed to provide pedigrees of models and model-based inferences.[54] Sensitivity auditing has been especially designed for an adversarial context, where not only the nature of the evidence, but also the degree of certainty and uncertainty associated to the evidence, will be the subject of partisan interests.[55] Sensitivity auditing is recommended in the European Commission guidelines for impact assessment,[8] as well as in the report Science Advice for Policy by European Academies.[56]
Related concepts
Sensitivity analysis is closely related with uncertainty analysis;[2] while the latter studies the overall uncertainty in the conclusions of the study, sensitivity analysis tries to identify what source of uncertainty weighs more on the study's conclusions.
The problem setting in sensitivity analysis also has strong similarities with the field of design of experiments.[57] In a design of experiments, one studies the effect of some process or intervention (the 'treatment') on some objects (the 'experimental units'). In sensitivity analysis one looks at the effect of varying the inputs of a mathematical model on the output of the model itself. In both disciplines one strives to obtain information from the system with a minimum of physical or numerical experiments.
See also
- Causality
- Elementary effects method
- Experimental uncertainty analysis
- Fourier amplitude sensitivity testing
- Info-gap decision theory
- Interval FEM
- Perturbation analysis
- Probabilistic design
- Probability bounds analysis
- Robustification
- ROC curve
- Uncertainty quantification
- Variance-based sensitivity analysis
References
- Saltelli, A. (2002). "Sensitivity Analysis for Importance Assessment". Risk Analysis. 22 (3): 1–12. CiteSeerX 10.1.1.194.7359. doi:10.1111/0272-4332.00040. PMID 12088235.
- Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. (2008). Global Sensitivity Analysis: The Primer. John Wiley & Sons.
- Pannell, D. J. (1997). "Sensitivity Analysis of Normative Economic Models: Theoretical Framework and Practical Strategies" (PDF). Agricultural Economics. 16 (2): 139–152. doi:10.1016/S0169-5150(96)01217-0.
- Bahremand, A.; De Smedt, F. (2008). "Distributed Hydrological Modeling and Sensitivity Analysis in Torysa Watershed, Slovakia". Water Resources Management. 22 (3): 293–408. doi:10.1007/s11269-007-9168-x. S2CID 9710579.
- Hill, M.; Kavetski, D.; Clark, M.; Ye, M.; Arabi, M.; Lu, D.; Foglia, L.; Mehl, S. (2015). "Practical use of computationally frugal model analysis methods". Groundwater. 54 (2): 159–170. doi:10.1111/gwat.12330. OSTI 1286771. PMID 25810333.
- Hill, M.; Tiedeman, C. (2007). Effective Groundwater Model Calibration, with Analysis of Data, Sensitivities, Predictions, and Uncertainty. John Wiley & Sons.
- Der Kiureghian, A.; Ditlevsen, O. (2009). "Aleatory or epistemic? Does it matter?". Structural Safety. 31 (2): 105–112. doi:10.1016/j.strusafe.2008.06.020.
- http://ec.europa.eu/governance/impact/commission_guidelines/docs/iag_2009_en.pdf
- http://www.epa.gov/CREM/library/cred_guidance_0309.pdf
- A. Saltelli, G. Bammer, I. Bruno, E. Charters, M. Di Fiore, E. Didier, W. Nelson Espeland, J. Kay, S. Lo Piano, D. Mayo, R.J. Pielke, T. Portaluri, T.M. Porter, A. Puy, I. Rafols, J.R. Ravetz, E. Reinert, D. Sarewitz, P.B. Stark, A. Stirling, P. van der Sluijs, Jeroen P. Vineis, Five ways to ensure that models serve society: a manifesto, Nature 582 (2020) 482–484.
- Helton, J. C.; Johnson, J. D.; Salaberry, C. J.; Storlie, C. B. (2006). "Survey of sampling based methods for uncertainty and sensitivity analysis". Reliability Engineering and System Safety. 91 (10–11): 1175–1209. doi:10.1016/j.ress.2005.11.017.
- Tavakoli, Siamak; Mousavi, Alireza (2013). "Event tracking for real-time unaware sensitivity analysis (EventTracker)". IEEE Transactions on Knowledge and Data Engineering. 25 (2): 348–359. doi:10.1109/tkde.2011.240. S2CID 17551372.
- Tavakoli, Siamak; Mousavi, Alireza; Poslad, Stefan (2013). "Input variable selection in time-critical knowledge integration applications: A review, analysis, and recommendation paper". Advanced Engineering Informatics. 27 (4): 519–536. doi:10.1016/j.aei.2013.06.002.
- Saltelli, A.; Annoni, P. (2010). "How to avoid a perfunctory sensitivity analysis". Environmental Modeling and Software. 25 (12): 1508–1517. doi:10.1016/j.envsoft.2010.04.012.
- Paruolo, P.; Saisana, M.; Saltelli, A. (2013). "Ratings and Rankings: Voodoo or Science?". Journal of the Royal Statistical Society, Series A. 176 (3): 609–634. arXiv:1104.3009. doi:10.1111/j.1467-985X.2012.01059.x. S2CID 54074392.
- Leamer, Edward E. (1983). "Let's Take the Con Out of Econometrics". American Economic Review. 73 (1): 31–43. JSTOR 1803924.
- Leamer, Edward E. (1985). "Sensitivity Analyses Would Help". American Economic Review. 75 (3): 308–313. JSTOR 1814801.
- Ravetz, J.R., 2007, No-Nonsense Guide to Science, New Internationalist Publications Ltd.
- Tsvetkova, O.; Ouarda, T.B.M.J. (2019). "Quasi-Monte Carlo technique in global sensitivity analysis of wind resource assessment with a study on UAE". J. Renew. Sustain. Energy. 11 (5): 053303. doi:10.1063/1.5120035.
- Saltelli, A.; Aleksankina, K.; Becker, W.; Fennell, P.; Ferretti, F.; Holst, N.; Li, S.; Wu, Q. (2019). "Why so many published sensitivity analyses are false: A systematic review of sensitivity analysis practices". Environ. Model. Softw. 114: 29–39. doi:10.1016/J.ENVSOFT.2019.01.012.
- O'Hagan, A.; et al. (2006). Uncertain Judgements: Eliciting Experts' Probabilities. Chichester: Wiley. ISBN 9780470033302.
- Sacks, J.; Welch, W. J.; Mitchell, T. J.; Wynn, H. P. (1989). "Design and Analysis of Computer Experiments". Statistical Science. 4 (4): 409–435. doi:10.1214/ss/1177012413.
- Campbell, J.; et al. (2008). "Photosynthetic Control of Atmospheric Carbonyl Sulfide During the Growing Season". Science. 322 (5904): 1085–1088. Bibcode:2008Sci...322.1085C. doi:10.1126/science.1164015. PMID 19008442. S2CID 206515456.
- Bailis, R.; Ezzati, M.; Kammen, D. (2005). "Mortality and Greenhouse Gas Impacts of Biomass and Petroleum Energy Futures in Africa". Science. 308 (5718): 98–103. Bibcode:2005Sci...308...98B. doi:10.1126/science.1106881. PMID 15802601. S2CID 14404609.
- Murphy, J.; et al. (2004). "Quantification of modelling uncertainties in a large ensemble of climate change simulations". Nature. 430 (7001): 768–772. Bibcode:2004Natur.430..768M. doi:10.1038/nature02771. PMID 15306806. S2CID 980153.
- Czitrom, Veronica (1999). "One-Factor-at-a-Time Versus Designed Experiments". American Statistician. 53 (2): 126–131. doi:10.2307/2685731. JSTOR 2685731.
- Cacuci, Dan G. Sensitivity and Uncertainty Analysis: Theory. I. Chapman & Hall.
- Cacuci, Dan G.; Ionescu-Bujor, Mihaela; Navon, Michael (2005). Sensitivity and Uncertainty Analysis: Applications to Large-Scale Systems. II. Chapman & Hall.
- Griewank, A. (2000). Evaluating Derivatives, Principles and Techniques of Algorithmic Differentiation. SIAM.
- Kabir HD, Khosravi A, Nahavandi D, Nahavandi S. Uncertainty Quantification Neural Network from Similarity and Sensitivity. In2020 International Joint Conference on Neural Networks (IJCNN) 2020 Jul 19 (pp. 1-8). IEEE.
- Sobol', I (1990). "Sensitivity estimates for nonlinear mathematical models". Matematicheskoe Modelirovanie (in Russian). 2: 112–118.; translated in English in Sobol', I (1993). "Sensitivity analysis for non-linear mathematical models". Mathematical Modeling & Computational Experiment. 1: 407–414.
- Homma, T.; Saltelli, A. (1996). "Importance measures in global sensitivity analysis of nonlinear models". Reliability Engineering and System Safety. 52: 1–17. doi:10.1016/0951-8320(96)00002-6.
- Saltelli, A.; Chan, K.; and Scott, M. (eds.) (2000). Sensitivity Analysis. Wiley Series in Probability and Statistics. New York: John Wiley and Sons.
- Saltelli, A.; Tarantola, S. (2002). "On the relative importance of input factors in mathematical models: safety assessment for nuclear waste disposal". Journal of the American Statistical Association. 97 (459): 702–709. doi:10.1198/016214502388618447. S2CID 59463173.
- Razavi, Saman; Gupta, Hoshin V. (January 2016). "A new framework for comprehensive, robust, and efficient global sensitivity analysis: 1. Theory". Water Resources Research. 52 (1): 423–439. Bibcode:2016WRR....52..423R. doi:10.1002/2015WR017558. ISSN 1944-7973.
- Razavi, Saman; Gupta, Hoshin V. (January 2016). "A new framework for comprehensive, robust, and efficient global sensitivity analysis: 2. Application". Water Resources Research. 52 (1): 440–455. Bibcode:2016WRR....52..440R. doi:10.1002/2015WR017559. ISSN 1944-7973.
- Haghnegahdar, Amin; Razavi, Saman (September 2017). "Insights into sensitivity analysis of Earth and environmental systems models: On the impact of parameter perturbation scale". Environmental Modelling & Software. 95: 115–131. doi:10.1016/j.envsoft.2017.03.031.
- Gupta, H; Razavi, S (2016). "Challenges and Future Outlook of Sensitivity Analysis". In Petropoulos, George; Srivastava, Prashant (eds.). Sensitivity Analysis in Earth Observation Modelling (1st ed.). pp. 397–415. ISBN 9780128030318.
- Morris, M. D. (1991). "Factorial sampling plans for preliminary computational experiments". Technometrics. 33 (2): 161–174. CiteSeerX 10.1.1.584.521. doi:10.2307/1269043. JSTOR 1269043.
- Campolongo, F.; Cariboni, J.; Saltelli, A. (2007). "An effective screening design for sensitivity analysis of large models". Environmental Modelling and Software. 22 (10): 1509–1518. doi:10.1016/j.envsoft.2006.10.004.
- Storlie, C.B.; Swiler, L.P.; Helton, J.C.; Sallaberry, C.J. (2009). "Implementation and evaluation of nonparametric regression procedures for sensitivity analysis of computationally demanding models". Reliability Engineering & System Safety. 94 (11): 1735–1763. doi:10.1016/j.ress.2009.05.007.
- Wang, Shangying; Fan, Kai; Luo, Nan; Cao, Yangxiaolu; Wu, Feilun; Zhang, Carolyn; Heller, Katherine A.; You, Lingchong (2019-09-25). "Massive computational acceleration by using neural networks to emulate mechanism-based biological models". Nature Communications. 10 (1): 4354. doi:10.1038/s41467-019-12342-y. ISSN 2041-1723. PMC 6761138. PMID 31554788.
- Oakley, J.; O'Hagan, A. (2004). "Probabilistic sensitivity analysis of complex models: a Bayesian approach". J. Royal Stat. Soc. B. 66 (3): 751–769. CiteSeerX 10.1.1.6.9720. doi:10.1111/j.1467-9868.2004.05304.x.
- Gramacy, R. B.; Taddy, M. A. (2010). "Categorical Inputs, Sensitivity Analysis, Optimization and Importance Tempering with tgp Version 2, an R Package for Treed Gaussian Process Models" (PDF). Journal of Statistical Software. 33 (6). doi:10.18637/jss.v033.i06.
- Becker, W.; Worden, K.; Rowson, J. (2013). "Bayesian sensitivity analysis of bifurcating nonlinear models". Mechanical Systems and Signal Processing. 34 (1–2): 57–75. Bibcode:2013MSSP...34...57B. doi:10.1016/j.ymssp.2012.05.010.
- Sudret, B. (2008). "Global sensitivity analysis using polynomial chaos expansions". Reliability Engineering & System Safety. 93 (7): 964–979. doi:10.1016/j.ress.2007.04.002.
- Ratto, M.; Pagano, A. (2010). "Using recursive algorithms for the efficient identification of smoothing spline ANOVA models". AStA Advances in Statistical Analysis. 94 (4): 367–388. doi:10.1007/s10182-010-0148-8. S2CID 7678955.
- Li, G.; Hu, J.; Wang, S.-W.; Georgopoulos, P.; Schoendorf, J.; Rabitz, H. (2006). "Random Sampling-High Dimensional Model Representation (RS-HDMR) and orthogonality of its different order component functions". Journal of Physical Chemistry A. 110 (7): 2474–2485. Bibcode:2006JPCA..110.2474L. doi:10.1021/jp054148m. PMID 16480307.
- Li, G. (2002). "Practical approaches to construct RS-HDMR component functions". Journal of Physical Chemistry. 106 (37): 8721–8733. doi:10.1021/jp014567t.
- Rabitz, H (1989). "System analysis at molecular scale". Science. 246 (4927): 221–226. Bibcode:1989Sci...246..221R. doi:10.1126/science.246.4927.221. PMID 17839016. S2CID 23088466.
- Hornberger, G.; Spear, R. (1981). "An approach to the preliminary analysis of environmental systems". Journal of Environmental Management. 7: 7–18.
- Saltelli, A.; Tarantola, S.; Campolongo, F.; Ratto, M. (2004). Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models. John Wiley and Sons.
- Van der Sluijs, JP; Craye, M; Funtowicz, S; Kloprogge, P; Ravetz, J; Risbey, J (2005). "Combining quantitative and qualitative measures of uncertainty in model based environmental assessment: the NUSAP system". Risk Analysis. 25 (2): 481–492. doi:10.1111/j.1539-6924.2005.00604.x. hdl:1874/386039. PMID 15876219. S2CID 15988654.
- Saltelli, A.; van der Sluijs, J.; Guimarães Pereira, Â. (2013). "Funtowiz, S.O., What do I make of your Latinorum? Sensitivity auditing of mathematical modelling". International Journal Foresight and Innovation Policy. 9: 213–234. arXiv:1211.2668. doi:10.1504/ijfip.2013.058610. S2CID 55591748.
- Lo Piano, S; Robinson, M (2019). "Nutrition and public health economic evaluations under the lenses of post normal science". Futures. 112: 102436. doi:10.1016/j.futures.2019.06.008.
- Science Advice for Policy by European Academies, Making sense of science for policy under conditions of complexity and uncertainty, Berlin, 2019.
- Box GEP, Hunter WG, Hunter, J. Stuart. Statistics for experimenters [Internet]. New York: Wiley & Sons
Further reading
- Cannavó, F. (2012). "Sensitivity analysis for volcanic source modeling quality assessment and model selection". Computers & Geosciences. 44: 52–59. Bibcode:2012CG.....44...52C. doi:10.1016/j.cageo.2012.03.008.
- Fassò A. (2007) "Statistical sensitivity analysis and water quality". In Wymer L. Ed, Statistical Framework for Water Quality Criteria and Monitoring. Wiley, New York.
- Fassò A., Perri P.F. (2002) "Sensitivity Analysis". In Abdel H. El-Shaarawi and Walter W. Piegorsch (eds) Encyclopedia of Environmetrics, Volume 4, pp 1968–1982, Wiley.
- Fassò A., Esposito E., Porcu E., Reverberi A.P., Vegliò F. (2003) "Statistical Sensitivity Analysis of Packed Column Reactors for Contaminated Wastewater". Environmetrics. Vol. 14, n.8, 743–759.
- Haug, Edward J.; Choi, Kyung K.; Komkov, Vadim (1986) Design sensitivity analysis of structural systems. Mathematics in Science and Engineering, 177. Academic Press, Inc., Orlando, FL.
- Pianosi, F.; Beven, K.; Freer, J.; Hall, J.W.; Rougier, J.; Stephenson, D.B.; Wagener, T. (2016). "Sensitivity analysis of environmental models: A systematic review with practical workflow". Environmental Modeling and Software. 79: 214–232. doi:10.1016/j.envsoft.2016.02.008.
- Pilkey, O. H. and L. Pilkey-Jarvis (2007), Useless Arithmetic. Why Environmental Scientists Can't Predict the Future. New York: Columbia University Press.
- Santner, T. J.; Williams, B. J.; Notz, W.I. (2003) Design and Analysis of Computer Experiments; Springer-Verlag.
- Taleb, N. N., (2007) The Black Swan: The Impact of the Highly Improbable, Random House.
Special issues
- Reliability Engineering & System Safety, 2003, 79:121–2: SAMO 2001: Methodological advances and innovative applications of sensitivity analysis , edited by Tarantola S, Saltelli.
- Reliability Engineering and System Safety,Volume 91, 2006, Special issue on sensitivity analysis, edited by Helton JC, Cooke RM, McKay MD, Saltelli.
- International Journal of Chemical Kinetics 2008, Volume 40, Issue 11 – Special Issue on Sensitivity Analysis,edited by Turányi T.
- Reliability Engineering & System Safety, Volume 94, Issue 7,Pages 1133-1244 (July 2009), Special Issue on Sensitivity Analysis, edited by Andrea Saltelli.
- Reliability Engineering & System Safety, Volume 107, November 2012, Advances in sensitivity analysis, SAMO 2010, edited by Borgonovo E, Tarantola S.
- Journal of Statistical Computation and Simulation Volume 85, 2015 - Issue 7: Special Issue: Selected Papers from the 7th International Conference on Sensitivity Analysis of Model Output, July 2013, Nice, France, edited by David Ginsbourger, Bertrand Iooss & Luc Pronzato.
- Reliability Engineering & System Safety, Volume 134, February 2015, edited by Stefano Tarantola, and Nathalie Saint-Geours.
- Environmental Modelling & Software, Special issue: Sensitivity analysis for environmental modelling (2020), Edited by Saman Razavi, Andrea Saltelli, Tony Jakeman, Qiongli Wu.
External links
- Joseph Hart, Julie Bessac, Emil Constantinescu (2018), "Global sensitivity analysis for statistical model parameters", arXiv:1708.07441
- web-page on Sensitivity analysis – (Joint Research Centre of the European Commission)
- SimLab, the free software for global sensitivity analysis of the Joint Research Centre
- MUCM Project – Extensive resources for uncertainty and sensitivity analysis of computationally-demanding models.