comparePlot: Compare different elicitated densities.

Description

Compare different elicitated distributions in a trellis display.

Usage

comparePlot(..., type = c("density", "cdf"), deriv,
            points = TRUE, superpose = FALSE, n = 101,
            xlab= "", ylab, addArgs = NULL)

Arguments

...: Fitted SEL objects separated by a comma.
type: Determines whether to plot densities or cdfs (ignored if deriv is not missing).
deriv: Specify derivative of the expert's density to be plotted.
points: Logical indicating whether elicited quantiles should be displayed, when displaying the cdf.
superpose: Logical indicate if plots should be superposed.
n: Integer, number of points used for plotting.
xlab, ylab: Labels for x-axis and y-axis. If ylab is missing the value of type is used.
addArgs: List specifying additional arguments for xyplot function, list elements should match the corresponding arguments of the xyplot function.

Author

Bjoern Bornkamp

References

Bornkamp, B. and Ickstadt, K. (2009). A Note on B-Splines for Semiparametric Elicitation. The American Statistician, 63, 373--377

Examples

Run this code

# example from O'Hagan et al. (2006)
x <- c(177.5, 183.75, 190, 205, 220)
y <- c(0.175, 0.33, 0.5, 0.75, 0.95)

default   <- SEL(x, y, Delta = 0.05, bounds = c(165, 250))
bernst    <- SEL(x, y, d = 10, N = 0, Delta = 0.05, bounds = c(165, 250))
unifknots <- SEL(x, y, d = 3, N = 5, Delta = 0.05, bounds = c(165, 250))
lin       <- SEL(x, y, d = 1, inknts = x, Delta = 0.05, bounds = c(165, 250))
comparePlot(default, bernst, unifknots, lin, type = "cdf")
comparePlot(default, bernst, unifknots, lin, type = "density")