Plotting SAFI model coefficients

Share:

Description

Barplots of the functional sensitivities from a safimodel object.

Usage

1
2
## S3 method for class 'safimodel'
plot(x, ylim = NULL, ...)

Arguments

x

safimodel object

ylim

limits for the y-axis. If NULL, limits are computed internally.

...

optional graphical parameters

Value

One plot for each functional input.

References

Fruth, J.; Roustant, O.; Kuhnt, S. (2014) Sequential designs for sensitivity analysis of functional inputs in computer experiments, Reliability Engineering & System Safety, doi: 10.1016/j.ress.2014.07.018, preprint on HAL: http://hal.archives-ouvertes.fr/hal-00943509.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
### simple example

s.d <- createSafiDesign(d.f = 1)
s.d2 <- splitSafiDesign(s.d = s.d, new.split.points = list(c(0.25, 0.75)))

# artificial model output (rising influence)
x <- accessSafiDesign(s.d = s.d2, n.timepoints = 4)
y <- x$x1 %*% c(0,1,2,3)
s.m <- safiModel(s.d = s.d2, y = y)
plot(s.m)


### d.f = 3, mirrored

s.d <- createSafiDesign(d.f = 3, mirrored.runs.included = TRUE)
s.d2 <- splitSafiDesign(s.d, list(c(0.5), c(0.25, 0.75), c(0.25, 0.5, 0.75)))

# artificial model output (x1 without influence, x2 rising, x3 falling)
x <- accessSafiDesign(s.d = s.d2, n.timepoints = 4)
y <- x$x2 %*% c(0, 1, 2, 3) + x$x3 %*% c(0, -1, -2, -3)
s.m <- safiModel(s.d2, y = y)
plot(s.m)