comparePlot {SEL}

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(s)

Bjoern Bornkamp

References

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

See Also

SEL, plot.SEL

Examples

# 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")