Nothing
## 2021-11-19: se incluyo la opcion de mostrar la incertidumbre expandida o la incertidumbre estandar
plot.comparison <- function(x, ...) {
# x is a data.frame with method:string, mean:numeric, u:numeric
display.order <- getOption("display.order")
display.shownames <- getOption("display.shownames")
display.orientation <- getOption("display.orientation")
display.length.out <- getOption("display.length.out")
display.tab.size <- getOption("display.tab.size")
display.signif.digits <- getOption("display.signif.digits")
display.expandedUncertainty <- getOption("display.expandedUncertainty")
if (is.null(display.order)) display.order <- "location"
if (is.null(display.shownames)) display.shownames <- FALSE
if (is.null(display.orientation)) display.orientation <- "horizontal"
if (is.null(display.length.out)) display.length.out <- 101
if (is.null(display.tab.size)) display.tab.size <- 12
if (is.null(display.signif.digits)) display.signif.digits <- 2
if (is.null(display.expandedUncertainty)) display.expandedUncertainty <- FALSE
n <- length(x$fit$value)
ci <- matrix(cbind(x$fit$value, x$fit$expandedUnc/x$fit$coverageFactor, x$fit$expandedUnc), n, 3,
byrow = FALSE)
if (display.order == "location") {
ss <- order(x$fit$value)
} else if (display.order == "dispersion") {
ss <- order(x$fit$expandedUnc)
} else {
ss <- order(x$fit$method)
}
zz <- c(1:n)
zlim <- c(1, n)
ww <- x$fit$value[ss]
wlim <- range(c(x$fit$value - x$fit$expandedUnc, x$fit$value + x$fit$expandedUnc))
ssxlab <- "method"
wlab <- paste0(x$gconsensus$exercise, ": ", x$gconsensus$measurand, " /(", x$fit$unit[1], ")")
if (display.orientation == "horizontal") {
xx <- zz
yy <- ww
xlim <- zlim
ylim <- wlim
xaxis <- 1
yaxis <- 2
xlab <- ssxlab
ylab <- wlab
} else {
xx <- ww
yy <- zz
xlim <- wlim
ylim <- zlim
xaxis <- 2
yaxis <- 1
xlab <- wlab
ylab <- ssxlab
}
plot(xx[ss],
yy[ss],
xlim = xlim,
ylim = ylim,
axes = FALSE,
xlab = xlab,
main = expression("Comparing Consensus Values " * (1 - alpha) * "% CI"),
ylab = ylab
)
if (display.orientation == "horizontal") {
if (display.expandedUncertainty) {
for (i in 1:n) {
lines(rep(i, 2), x$fit$value[ss][i] + c(-1, 1) * x$fit$expandedUnc[ss][i], lwd=1)
}
} else {
for (i in 1:n) {
lines(rep(i, 2), x$fit$value[ss][i] + c(-1, 1) * x$fit$expandedUnc[ss][i]/x$fit$coverageFactor[ss][i], lwd=1)
}
}
} else {
if (display.expandedUncertainty) {
for (i in 1:n) {
lines(x$fit$value[ss][i] + c(-1, 1) * x$fit$expandedUnc[ss][i], rep(i, 2), lwd=1)
}
} else {
for (i in 1:n) {
lines(x$fit$value[ss][i] + c(-1, 1) * x$fit$expandedUnc[ss][i]/x$fit$coverageFactor[ss][i], rep(i, 2), lwd=1)
}
}
}
box()
axis(yaxis)
axis(xaxis, at = c(1:n), labels = x$fit$method[ss])
}
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.