1. What Is Sensitivity Analysis?
Sensitivity analysis examines how changes in model inputs affect model outputs. In QSP, PK/PD, and other mechanistic models, the inputs may include parameters such as clearance, receptor abundance, binding affinity, production rates, degradation rates, or pharmacodynamic potency.
Suppose a model produces an output \(Y\) as a function of parameters \(\boldsymbol{\theta}\):
Sensitivity analysis asks questions such as:
- Which parameters have the largest influence on the model output?
- Which parameters have little influence over the region being studied?
- Does parameter influence depend on the values of other parameters?
- Are there interactions between parameters?
- Which measurements or experiments could most reduce uncertainty in an important prediction?
2. Why Sensitivity Analysis Matters in QSP
QSP models can contain dozens, hundreds, or even thousands of parameters. Not every parameter contributes equally to every model output.
Sensitivity analysis provides a systematic way to understand this complexity. It can help distinguish parameters that strongly affect a specific prediction from parameters whose influence is negligible within the scenario being studied.
| Use | Question | Potential application |
|---|---|---|
| Model understanding | Which parameters control an output? | Interpret the mechanisms driving a prediction |
| Experimental design | Which uncertain quantities matter most? | Prioritize measurements |
| Model reduction | Which parameters have little influence? | Identify opportunities for simplification |
| Uncertainty analysis | Which inputs contribute most to output uncertainty? | Prioritize uncertainty reduction |
| Drug development | Which biological processes affect treatment response? | Explore mechanistic hypotheses |
| Robustness assessment | Does the conclusion change across plausible inputs? | Assess confidence in model-based conclusions |
The distinction between local and global sensitivity analysis is particularly important because they explore parameter influence in fundamentally different ways.
3. What Is Local Sensitivity Analysis?
Local sensitivity analysis evaluates how an output changes when parameters are perturbed around a particular nominal parameter set.
For a scalar output \(Y=f(\boldsymbol{\theta})\), the local sensitivity of \(Y\) to parameter \(\theta_i\) can be represented by the partial derivative:
The derivative describes the instantaneous change in the output associated with a small change in the parameter, while the other parameters are held fixed.
For dynamic models, sensitivity can be evaluated at each time point. If \(Y(t)\) is the model output:
This produces a sensitivity trajectory rather than a single sensitivity value.
4. Normalized Local Sensitivity
Raw derivatives can be difficult to compare because different parameters may have different units and numerical scales.
A commonly useful normalized sensitivity is:
Approximately, this quantity represents the fractional change in output associated with a fractional change in the parameter.
For example, a normalized sensitivity of \(0.8\) means that, locally, a 1% increase in the parameter corresponds to approximately a 0.8% increase in the output, under the local linear approximation.
A negative value indicates an inverse local relationship: increasing the parameter tends to decrease the output near the nominal parameter set.
| Normalized sensitivity | Local interpretation |
|---|---|
| Large positive | Output increases relatively strongly when the parameter increases |
| Large negative | Output decreases relatively strongly when the parameter increases |
| Near zero | Output is locally insensitive to the parameter |
The magnitude should be interpreted in the context of the model, output, parameterization, and perturbation region.
5. How Is Local Sensitivity Calculated?
The derivative can be obtained analytically, through sensitivity equations, automatic differentiation, or numerical perturbation.
A simple finite-difference approximation is:
Central finite differences often provide a useful numerical approximation when analytic derivatives are unavailable.
For a dynamic system represented by:
one can define a sensitivity state:
and propagate that sensitivity through the system alongside the model states.
6. What Is Global Sensitivity Analysis?
Global sensitivity analysis (GSA) evaluates parameter influence across a specified region of the parameter space rather than at only one nominal parameter set.
Instead of asking:
global analysis asks:
This distinction is important for nonlinear QSP models because parameter effects can change substantially across the parameter space.
Local analysis examines behavior near a nominal parameter set; global analysis explores influence across a specified region of plausible parameter values.
7. Common Global Sensitivity Methods
Morris screening
The Morris method, also called the elementary-effects method, evaluates changes in the output produced by perturbations at multiple locations in the parameter space.
It is particularly useful as a screening method when the model contains many parameters and a full variance-based analysis would be computationally expensive.
Sobol sensitivity indices
Sobol analysis decomposes output variance into contributions associated with individual parameters and their interactions.
The first-order Sobol index for parameter \(i\) can be written conceptually as:
The total-order index is commonly represented as:
Here \(\theta_{-i}\) denotes all parameters except \(\theta_i\).
The first-order index measures the contribution of parameter \(i\) acting alone, whereas the total-order index includes its contribution through interactions with other parameters.
Other approaches
Other global approaches include Latin hypercube sampling followed by regression or correlation-based analysis, FAST and extended FAST, variance decomposition, and distribution-based methods.
8. Parameter Interactions
One of the major advantages of global sensitivity analysis is that it can reveal interactions between parameters.
Suppose two parameters have little influence individually but have a substantial joint effect. A purely one-at-a-time local analysis may fail to reveal that relationship.
For Sobol analysis, the difference between total-order and first-order indices provides useful information:
A relatively large difference indicates that parameter \(i\)'s influence involves interactions with other parameters.
9. Local vs. Global Sensitivity Analysis
| Feature | Local sensitivity | Global sensitivity |
|---|---|---|
| Parameter region | Near a nominal point | Across a specified parameter space |
| Typical perturbation | Small | Broad, according to defined distributions or ranges |
| Common output | Derivative or normalized sensitivity | Variance contribution, elementary effects, or related measures |
| Nonlinearity | May be poorly represented away from the nominal point | Can capture nonlinear effects across the explored region |
| Interactions | Usually not characterized directly | Can be quantified by suitable methods |
| Computational cost | Often relatively low | Can be substantial for complex models |
| Best suited for | Local interpretation and mechanistic diagnostics | Parameter screening, uncertainty attribution, and robustness analysis |
Neither method is universally preferable. They answer different questions and can be complementary.
10. Worked Example: A Simple Nonlinear Model
Consider a hypothetical response model:
Suppose \(C=10\), \(\theta_1=100\), and \(\theta_2=10\).
Step 1: Calculate the output
Step 2: Local sensitivity to \(\theta_1\)
Because \(Y\) is linear in \(\theta_1\):
The normalized sensitivity is:
Thus, locally, a 1% increase in \(\theta_1\) produces approximately a 1% increase in \(Y\).
Step 3: Local sensitivity to \(\theta_2\)
At the nominal values:
The normalized sensitivity is:
Locally, a 1% increase in \(\theta_2\) produces approximately a 0.5% decrease in \(Y\).
Step 4: Why a global analysis can tell us more
The effect of \(\theta_2\) depends strongly on the relationship between \(\theta_2\) and \(C\). If \(\theta_2\) becomes much larger than \(C\), the model behaves differently than when \(\theta_2\) is much smaller than \(C\).
A global analysis would therefore explore many plausible values of \(\theta_1\) and \(\theta_2\), rather than relying only on the sensitivities at \(100\) and \(10\).
11. Sensitivity Analysis in QSP Models
QSP models are especially well suited to sensitivity analysis because they connect biological mechanisms across multiple scales.
For example, a QSP model might include:
- Target expression and turnover.
- Drug concentration and tissue distribution.
- Receptor binding and target engagement.
- Signal transduction.
- Biomarker production and degradation.
- Disease progression.
- Immune-cell dynamics.
- Clinical or translational endpoints.
A parameter can be highly influential for one output but almost irrelevant for another.
| Output | Potential sensitive mechanisms |
|---|---|
| Plasma concentration | Clearance, bioavailability, distribution parameters |
| Target engagement | Binding affinity, target abundance, drug concentration |
| Biomarker response | Production, turnover, pathway parameters, target engagement |
| Disease burden | Growth/decay rates, treatment effects, disease-state parameters |
| Clinical response | Multiple upstream mechanisms and parameter interactions |
This is why sensitivity analysis should generally be performed with respect to a clearly defined model output and scenario.
12. Defining the Global Parameter Space
A global sensitivity analysis is only as meaningful as the parameter space that it explores.
For each uncertain parameter, the analyst needs to specify a plausible range or probability distribution.
For example:
| Parameter | Example range | Possible rationale |
|---|---|---|
| Clearance | 5–15 L/h | Observed or plausible PK variability |
| Binding affinity | 1–100 nM | Experimental uncertainty |
| Production rate | 0.5–2× nominal | Biological variability |
| Turnover rate | 0.25–4× nominal | Uncertain kinetic parameter |
Probability distributions may be preferable when the uncertainty has a meaningful probabilistic interpretation. Uniform ranges can be useful for bounded exploration but should not automatically be interpreted as probability statements.
13. Sensitivity Depends on the Output
There is rarely one universal sensitivity ranking for a complex model.
Consider a QSP model with outputs such as:
- Maximum biomarker response.
- Time to maximum response.
- Area under the biomarker-time curve.
- Minimum disease burden.
- Probability of achieving a response threshold.
Each output can have a different sensitivity profile.
For example, one parameter may strongly influence the magnitude of response while having little effect on the timing of response.
For dynamic models, sensitivity can therefore be assessed at specific time points, over integrated summaries, or using clinically meaningful endpoint definitions.
14. Sensitivity Analysis Is Not the Same as Uncertainty Analysis
The concepts are closely related but distinct.
Sensitivity asks how strongly the output responds to changes in an input.
Uncertainty analysis asks how uncertainty in the inputs propagates into uncertainty in the output.
A parameter can be highly sensitive but precisely known. In that case, it may contribute little to output uncertainty.
Conversely, a parameter can have moderate sensitivity but very large uncertainty and therefore contribute substantially to predictive uncertainty.
| Parameter | Sensitivity | Parameter uncertainty | Potential impact on output uncertainty |
|---|---|---|---|
| A | High | Low | Potentially modest |
| B | Moderate | High | Potentially substantial |
| C | Low | High | Often limited |
This distinction is especially important when using sensitivity analysis to prioritize additional experiments.
15. Screening Before Detailed Global Analysis
Large QSP models can make exhaustive global analysis computationally expensive. A practical strategy is often to begin with a screening analysis.
- Define the scientific outputs of interest.
- Specify plausible parameter ranges.
- Run a computationally efficient screening method.
- Identify parameters with potentially important effects.
- Perform more detailed analysis on the influential subset.
- Investigate interactions and uncertainty contributions where appropriate.
The Morris method can be useful for this type of screening because it can identify parameters with strong, nonlinear, or interaction-related effects without requiring the same computational burden as some variance-based approaches.
16. A Practical Sensitivity Analysis Workflow
Step 1: Define the question
Specify exactly what you want to understand. For example: “Which parameters most influence predicted biomarker suppression at week 12?”
Step 2: Define the output
Choose a time point, summary statistic, threshold, or other scientifically meaningful endpoint.
Step 3: Define plausible parameter ranges
Use experimental data, literature, prior distributions, calibrated estimates, or other scientifically defensible information.
Step 4: Decide whether local or global analysis is appropriate
Local analysis is useful for behavior near a calibrated parameter set. Global analysis is more appropriate when broad parameter uncertainty, nonlinear behavior, or interactions are central to the question.
Step 5: Run the analysis
For local analysis, calculate derivatives or perturbation-based sensitivities. For global analysis, use an appropriate sampling and sensitivity method.
Step 6: Check numerical robustness
Assess whether the results are stable with respect to perturbation size, sampling density, solver tolerances, and other computational choices.
Step 7: Interpret biologically
Translate statistical or mathematical sensitivity into the context of the mechanism represented by the model.
Step 8: Use the results
Potential uses include prioritizing measurements, identifying influential mechanisms, assessing robustness, simplifying models, or selecting parameters for further investigation.
17. Common Mistakes in Sensitivity Analysis
1. Treating a local result as globally valid
A derivative at one parameter set does not establish that the same sensitivity holds across the entire plausible parameter space.
2. Ignoring parameter interactions
One-at-a-time perturbations can miss important joint effects in nonlinear models.
3. Comparing raw derivatives without considering scale
Parameters can have different units and magnitudes. Normalization can make comparisons more interpretable.
4. Using unrealistic parameter ranges
A mathematically broad range may produce a sensitivity analysis that is disconnected from biological reality.
5. Treating sensitivity as causal proof
A model parameter can be mathematically influential without establishing that the corresponding biological mechanism is uniquely responsible for an observed phenomenon.
6. Ignoring model failures
Sensitivity analysis cannot compensate for an inappropriate structural model, incorrect equations, or invalid assumptions.
7. Reporting rankings without the underlying context
A sensitivity ranking should specify the model version, output, parameter ranges, scenario, method, and relevant computational settings.
18. How Should Sensitivity Results Be Interpreted?
Sensitivity results are most useful when interpreted as conditional statements about a particular model and analysis setup.
For example:
This is more informative than simply stating that “A and B are the most important biological parameters.”
The distinction matters because a parameter's influence can change with:
- the model output;
- the treatment scenario;
- the time point;
- the parameter range;
- the parameter distribution;
- the values of other parameters;
- the model structure itself.
19. Using Local and Global Sensitivity Together
Local and global methods do not have to be competing approaches. In many modeling workflows, they provide complementary information.
A useful sequence is:
For example, local sensitivities can quickly reveal which parameters influence a calibrated model around its nominal solution. A global analysis can then determine whether those conclusions remain consistent across plausible parameter combinations.
When the two approaches agree, confidence in the robustness of the sensitivity pattern can increase. When they disagree, the disagreement can itself be informative, particularly when it reflects nonlinearities or parameter interactions.
20. Worked QSP Example: Biomarker Suppression
Consider a simplified QSP model in which drug exposure suppresses the production of a biomarker:
Suppose the model contains uncertain parameters \(k_{\mathrm{in}}\), \(k_{\mathrm{out}}\), \(I_{\max}\), and \(IC_{50}\).
Local question
At the calibrated parameter set, which parameter produces the largest change in week-12 biomarker concentration for a small perturbation?
A local sensitivity analysis answers this by evaluating quantities such as:
Global question
Across the plausible ranges of these four parameters, which parameters contribute most to variation in \(B(12)\)?
A global analysis could sample the four-dimensional parameter space and calculate first-order and total-order sensitivity measures.
Interpretation
The two analyses can produce different patterns. For example, \(IC_{50}\) may have modest local sensitivity at the nominal exposure but become highly influential when exposure varies across a broader range.
21. Using Sensitivity Analysis to Guide Experiments
Sensitivity analysis can help prioritize experiments, but the interpretation should distinguish model sensitivity from value of information.
If an uncertain parameter strongly influences an important prediction, measuring that parameter more accurately may be useful. But the practical value of the measurement also depends on:
- how uncertain the parameter currently is;
- how much the proposed experiment could reduce that uncertainty;
- how reliable the measurement would be;
- the cost and feasibility of the experiment;
- whether reducing uncertainty would materially change the decision or prediction.
Thus, sensitivity analysis can be an important input into experimental prioritization, but it is not by itself a complete decision-analysis framework.
22. Implementing Sensitivity Analysis Computationally
Sensitivity analysis can be implemented in many modeling environments, including R, Python, MATLAB, and specialized QSP or systems-biology platforms.
A basic computational workflow is:
- Define the model equations.
- Define the nominal parameter vector.
- Define the model output or outputs.
- For local analysis, calculate derivatives or perturb parameters around the nominal values.
- For global analysis, define parameter ranges or distributions and generate samples.
- Run the model for each sampled parameter set.
- Calculate the selected sensitivity measures.
- Visualize and interpret the results.
For computationally intensive QSP models, surrogate models, parallel computing, screening methods, or reduced-order approaches may be useful.
23. What Sensitivity Analysis Does Not Tell You Automatically
- It does not establish biological truth. Sensitivity is conditional on the model.
- It does not establish parameter identifiability. A parameter can strongly influence an output while remaining difficult to estimate from available data.
- It does not establish causality in the real system. The result describes the mathematical model.
- It does not eliminate structural uncertainty. Alternative model structures can produce different sensitivity patterns.
- It does not automatically identify the best experiment. Experimental value depends on uncertainty, feasibility, cost, and the scientific objective.
- It does not guarantee robust predictions. Sensitivity is one component of model evaluation.
24. Key Takeaways
- Local sensitivity analysis evaluates parameter influence near a specified nominal parameter set.
- Local sensitivity is commonly expressed using derivatives such as \(\partial Y/\partial\theta_i\), often with normalization to facilitate comparison.
- Global sensitivity analysis explores parameter influence across a defined region or probability distribution of parameter values.
- Global methods can capture nonlinear effects and, depending on the method, parameter interactions.
- Morris analysis is useful for screening, while Sobol methods provide variance-based measures of first-order and total-order influence.
- Sensitivity depends on the model output, scenario, parameter ranges, and model structure.
- Sensitivity and uncertainty are related but different: a highly sensitive parameter is not necessarily a major contributor to output uncertainty if that parameter is already well known.
- QSP models can exhibit strong nonlinearities and interactions, making global analysis particularly useful for understanding behavior across plausible biological states.
- Local and global sensitivity analyses are complementary rather than mutually exclusive.
- Sensitivity analysis can help prioritize experiments and uncertainty reduction, but it does not by itself establish biological causality, identifiability, or experimental value.
Where to Go Next
A natural next step is to study parameter identifiability and its relationship to sensitivity analysis. Sensitivity describes how model outputs respond to parameters, whereas identifiability asks whether those parameters can be uniquely and reliably inferred from available data.
For QSP models, the next progression can include structural versus practical identifiability, parameter estimation, uncertainty propagation, uncertainty decomposition, parameter estimation under sparse data, and global sensitivity analysis for complex nonlinear systems.