Usage
morris(model = NULL, factors, r, design, binf = 0, bsup = 1,
scale = TRUE, ...)
## S3 method for class 'morris':
tell(x, y = NULL, ...)
## S3 method for class 'morris':
print(x, ...)
## S3 method for class 'morris':
plot(x, identify = FALSE, atpen = FALSE, y_col = NULL,
y_dim3 = NULL, \dots)
## S3 method for class 'morris':
plot3d(x, alpha = c(0.2, 0), sphere.size = 1, y_col = NULL,
y_dim3 = NULL)Arguments
model
a function, or a model with a predict method,
defining the model to analyze.
factors
an integer giving the number of factors, or a vector of
character strings giving their names.
r
either an integer giving the number of repetitions of the design,
i.e. the number of elementary effect computed per factor, or a
vector of two integers c(r1, r2) for the space-filling
improvement (Campolongo et al.). In this ca
design
a list specifying the design type and its
parameters:
type = "oat"for Morris's OAT design (Morris 1991),
with the parameters:levels: either an integer specifying the number of
levels of the desi
binf
either an integer, specifying the minimum value for the
factors, or a vector for different values for each factor.
bsup
either an integer, specifying the maximum value for the
factors, or a vector for different values for each factor.
scale
logical. If TRUE, the input design of experiments is
scaled after building the design and before computing the elementary
effects so that all factors vary within the range [0,1]. For each factor,
the scaling is done relativel
x
a list of class "morris" storing the state of the
screening study (parameters, data, estimates).
y
a vector of model responses.
identify
logical. If TRUE, the user selects with the
mouse the factors to label on the $(\mu^*,\sigma)$
graph (see identify).
atpen
logical. If TRUE (and identify = TRUE), the
user-identified labels (more precisely: their lower-left corners) of the
factors are plotted at the place where the user had clicked (if near enough
to one of the fact
y_col
an integer defining the index of the column of x$y to be
used for plotting the corresponding Morris statistics $\mu^*$
and $\sigma$ (only applies if x$y is a matrix or an
array).
If set to NULL (
y_dim3
an integer defining the index in the third dimension of
x$y to be used for plotting the corresponding Morris statistics
$\mu^*$ and $\sigma$ (only applies if x$y is an
array).
If set to NULL (a
alpha
a vector of three values between 0.0 (fully transparent) and 1.0
(opaque) (see rgl.material). The first value is for the
cone, the second for the planes.
sphere.size
a numeric value, the scale factor for displaying the
spheres.
...
for morris: any other arguments for model which
are passed unchanged each time it is called. For plot.morris:
arguments to be passed to plot.default.
eqn
$\sigma$code
scale = TRUEpreformatted
mu <- sapply(1:dim(x$ee)[4], function(i){
apply(x$ee[, , , i, drop = FALSE], 3, function(M){
apply(M, 2, mean)
})
}, simplify = "array")
mu.star <- sapply(1:dim(x$ee)[4], function(i){
apply(abs(x$ee)[, , , i, drop = FALSE], 3, function(M){
apply(M, 2, mean)
})
}, simplify = "array")
sigma <- sapply(1:dim(x$ee)[4], function(i){
apply(x$ee[, , , i, drop = FALSE], 3, function(M){
apply(M, 2, sd)
})
}, simplify = "array")