Description Usage Arguments Details Warning Author(s) Source References Examples
Generates the Latin Hypercube sampling for uncertainty and sensitivity analyses.
1 2 3 4 5 6 7 8 9 10 11 |
model |
The function to be run, representing the model or simulation.
If |
factors |
The names of the input variables (used for naming the 'data' data.frame and in plotting) Either a vector of strings or a single number representing the number of factors |
N |
The size of the hypercube, i.e., how many samples are generated. Must be at least the number of factors plus 2. |
q |
The quantile functions to be used. If only one is provided, it will be used for all parameters. Defaults to "qunif". |
q.arg |
A list containing the arguments for the 'q' functions. Each parameter must be specified by a named list, containing all of the arguments for the quantile distribution. If unsupplied, default values for the parameters are used. |
res.names |
Optional: what are the names of the model results? (Used mainly for plotting) |
method |
Currently, two methods are supported. "random" generates a simple LH, with no modifications. "HL" (the default) generates a random LH, and subsequently corrects the correlation matrix using the Huntington & Lyrintzis method. |
opts |
Further options for the method used. The method HL supports the following options: ‘COR’ The desired correlation matrix between the model variables. If none is provided, the function will generate a zero-correlation Latin Hypercube. ‘eps’ The tolerance between the prescribed correlation and the actual correlation present in the generated Latin Hypercube. ‘maxIt’ The maximum number of iterations to be run for each factor. The default is set by a heuristic, but it might need some adjustments. |
nboot |
Number of bootstrap replicates for calculating the PRCC. |
repetitions |
The number of model repetitions to be run for a single data point. See the vignette on stochastic models for details |
cl |
Cluster generated with the “parallel” library. May be of any type supported. If a cluster is provided, the model will be run in parallel or distributed across the cluster via clusterApply. No load balancing is provided, so the model results are reproducible. NOTE: You should manually export ALL objects required for the model to run, including the model
function itself. See the help on |
x |
An LHS/PLUE object. For "tell", an incomplete LHS object (created with model=NULL) |
... |
Currently ignored |
y |
A data.frame containing the model responses |
A Latin Hypercube of size N
is generated from the desired quantile distribution functions
The following methods are currently supported for generating the LHS: random LHS and
Huntington & Lyrintzis method for correcting the correlation matrix to be similar
to the prescribed by the option COR
(see the arguments for description).
The specified model
is run with the data from the LHS. If repetitions
is set to more than one, the model will be run several times for each data point.
Partial rank correlation coefficients are estimated using code based on the prcc
function from the "sensitivity" package.
When the LHS function is called with no model (i.e., with argument
model=NULL
), it generates an incomplete object storing the Latin
Hypercube samples, and allowing the user to run the simulation
model independently. The method tell
allows to pass the simulation
results to the incomplete object.
tell
and ask
are S3 generic methods for decoupling
simulations and sensitivity measures estimations in the package
‘sensitivity’. In general, they are
not used by the end-user for a simple R model, but rather for an
external computational code. The LHS object implements only the
tell
method. For help on the other methods, see the help
pages on the ‘sensitivity’ package.
NOTE: the tell method from sensitivity objects (like 'fast99') modifies the object passed as argument as a side effect. This is NOT the case with the LHS tell method.
Andre Chalom
Uses internal code originally published on package sensitivity, by Gilles Pujol, Bertrand Iooss, Alexandre Janon
McKay, M.D. and Beckman, R.J. 1979. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code, Technometrics 21: 239-244
Chalom, A. and Prado, P.I.K.L. 2012. Parameter space exploration of ecological models arXiv:1210.6278 [q-bio.QM]
1 2 3 4 5 6 7 8 9 10 | completeLHS <- LHS(model=function(x) x[,1]+x[,2]*x[,3], factors=3, N=20)
incompleteLHS <- LHS(factors=5, N=30)
incompleteLHS <- tell(incompleteLHS, seq(1,30))
## Not run:
new.cluster <- parallel::makePSOCKcluster(c("localhost", "localhost"))
clusterLHS <- LHS(model=function(x) x[,1]/x[,2], factors=2, N=100, cl = new.cluster)
stopCluster(new.cluster)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.