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

View source: R/compareGroupRIO.R

compareGroupRIOR Documentation

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 <j.oakley@sheffield.ac.uk>

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)

SHELF documentation built on Sept. 11, 2024, 6:54 p.m.