compareGroupRIO: Compare individual elicited distributions with linear pool...

In OakleyJ/SHELF: Tools to Support the Sheffield Elicitation Framework

Compare individual elicited distributions with linear pool and RIO distribution

Description

Produce one of three plots to compare the individual elicited judgements with the final elicited distribution, chosen to represent the views of a "Rational Impartatial Observer" (RIO) as part of the SHELF process. A linear pool of fitted distributions from the individually elicited judgements is also obtained. The plot choices are a display of the quartiles, a display of the tertiles, and a plot of the various density functions.

Usage

compareGroupRIO(
  groupFit,
  RIOFit,
  type = "density",
  dLP = "best",
  dRIO = "best",
  xlab = "x",
  ylab = expression(f[X](x)),
  fs = 12
)

Arguments

groupFit	either an object of class elicitation, or the file path for a .csv file exported from the elicitMultiple() app. This should contain the individually elicited judgements from the experts
RIOFit	an object of class elicitation containing a single set of of probability judgements corresponding to the "Rational Impartial Observer (RIO)".
type	the plot used to show the comparison: one of "quartiles", "tertiles" or "density".
dLP	the distribution fitted to each expert's judgements and to the linear pool. Options are Options are "normal", "t", "gamma", "lognormal", "logt","beta", "mirrorgamma", "mirrorlognormal", "mirrorlogt" "hist" (for a histogram fit), and "best" (for best fitting).
dRIO	the distribution fitted to RIO's judgements. Options are the same as for dLP.
xlab	x-axis label in plot
ylab	y-axis label in plot
fs	font size used in plot

Author(s)

Jeremy Oakley <[e-mail]>

Examples

## Not run: 
l <- c(2, 1, 5, 1)
u <- c(95, 90, 65, 40)
v <- matrix(c(15, 25, 40,
 10, 20, 40,
 10, 15, 25,
 5, 10, 20),
 3, 4)
p <- c(0.25, 0.5, 0.75)
group <- fitdist(vals = v, probs = p, lower = l, upper = u)
rio <- fitdist(vals = c(12, 20, 25), probs = p, lower = 1, upper = 100)
compareGroupRIO(groupFit = group, RIOFit = rio, dRIO = "gamma")

## End(Not run)

OakleyJ/SHELF documentation built on March 17, 2024, 8:13 p.m.