`delsa`

implements Distributed Evaluation of
Local Sensitivity Analysis to calculate first order parameter
sensitivity at multiple locations in parameter space. The locations
in parameter space can either be obtained by a call to `parameterSets`

or by specifying `X0`

directly, in which case the prior variance
of each parameter `varprior`

also needs to be specified. Via `plot`

(which uses functions of the package `ggplot2`

and `reshape2`

),
the indices can be visualized.

1 2 3 4 5 6 7 8 9 10 11 12 13 |

`model` |
a function, or a model with a |

`perturb` |
Perturbation used to calculate sensitivity at each evaluation location |

`par.ranges` |
A named list of minimum and maximum parameter values |

`samples` |
Number of samples to generate. For the |

`method` |
Sampling scheme. See |

`X0` |
Parameter values at which to evaluate sensitivity indices.
Can be used instead of specifying sampling |

`varprior` |
Prior variance. If |

`...` |
any other arguments for |

`x` |
a list of class |

`y` |
a vector of model responses. |

`which` |
if a subset of the plots is required, specify a subset of the numbers 1:3 |

`ask` |
logical; if TRUE, the user is asked before each plot, see |

`print`

shows summary of the first order indices across parameter space.

`plot`

shows: (1) the cumulative distribution function of first order
sensitivity across parameter space, (2) variation of first order sensitivity
in relation to model response, and (3) sensitivity in relation to parameter value.

`delsa`

returns a list of class `"delsa"`

, containing all
the input arguments detailed before, plus the following components:

`call` |
the matched call. |

`X` |
a |

`y` |
a vector of model responses. |

`delsafirst` |
the first order indices for each location in |

Conversion for `sensitivity`

package by Joseph Guillaume,
based on original R code by Oldrich Rakovec

Rakovec, O., M. C. Hill, M. P. Clark, A. H. Weerts, A. J. Teuling, R. Uijlenhoet (2014), Distributed Evaluation of Local Sensitivity Analysis (DELSA), with application to hydrologic models, Water Resour. Res., 50, 1-18, doi:10.1002/2013WR014063

`parameterSets`

which is used to generate points, `sensitivity`

for other methods in the package

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ```
# Test case : the non-monotonic Sobol g-function
# (there are 8 factors, all following the uniform distribution on [0,1])
## Not run:
library(randtoolbox)
x <- delsa(model=sobol.fun,
par.ranges=replicate(8,c(0,1),simplify=FALSE),
samples=100,method="sobol")
# Summary of sensitivity indices of each parameter across parameter space
print(x)
library(ggplot2)
library(reshape2)
x11()
plot(x)
## End(Not run)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.