comparePlot {SEL}R Documentation

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, to appear

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

[Package SEL version 1.0-1 Index]