Nothing
## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----define v-----------------------------------------------------------------
v <- matrix(c(25, 30, 35, 30, 35, 50), nrow = 3, ncol = 2)
## ----show v, echo = F---------------------------------------------------------
v
## ----define p-----------------------------------------------------------------
p <- c(0.25, 0.5, 0.75)
## ----use myfit----------------------------------------------------------------
library(SHELF)
myfit <- fitdist(vals = v, probs = p, lower = 0, upper = 100)
## ----show myfit---------------------------------------------------------------
names(myfit)
## ----show beta----------------------------------------------------------------
myfit$Beta
## ----show ssq-----------------------------------------------------------------
myfit$ssq
## ----plot-fitted-distributions, fig.height = 4, fig.width = 5, fig.pos="h", fig.align="center", fig.cap = "The two fitted distributions and an equal-weighted linear pool."----
plotfit(myfit, lp = TRUE)
## ----define single set--------------------------------------------------------
v <-c(25, 30, 40)
p <-c(0.25, 0.5, 0.75)
consensus <- fitdist(vals = v, probs = p, lower = 0, upper = 100)
## ----plot-RIO, fig.height = 4, fig.width = 5, fig.pos="h", fig.align="center", fig.cap = "The fitted consensus distribution, with the lower and upper 5\\% tail areas shown as feedback."----
plotfit(consensus, ql = 0.05, qu = 0.95, d = "beta")
## ----show feedback------------------------------------------------------------
feedback(consensus, quantiles = c(0.05, 0.95))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.