View source: R/compareGroupRIO.R
compareGroupRIO | R Documentation |
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.
compareGroupRIO(
groupFit,
RIOFit,
type = "density",
dLP = "best",
dRIO = "best",
xlab = "x",
ylab = expression(f[X](x)),
fs = 12
)
groupFit |
either an object of class |
RIOFit |
an object of class |
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 |
xlab |
x-axis label in plot |
ylab |
y-axis label in plot |
fs |
font size used in plot |
Jeremy Oakley <j.oakley@sheffield.ac.uk>
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.