plot.summary.uncertainty: Plots an output quantity's uncertainty summary. It shows the...

View source: R/plot.summary.uncertainty.R

plot.summary.uncertaintyR Documentation

Plots an output quantity's uncertainty summary. It shows the uncertainty contribution from each involved input quantity

Description

Builds a barplot with a bar for each source of uncertainty. If correlation is present then an additional entry is added. The current metric used to display is When correlation is present its contribution may be negative.

Usage

## S3 method for class 'summary.uncertainty'
plot(x, y = NULL, ...)

Arguments

x

an uncertainty summary object

y

not used.

...

additional parameters to customize the plot

Details

none

Value

None (invisible NULL)

Note

none

Author(s)

Hugo Gasca-Aragon

Maintainer: Hugo Gasca-Aragon <hugo_gasca_aragon@hotmail.com>

References

JCGM 200:2012. International vocabulary of metrology—Basic and general concepts and associated terms (VIM)

JCGM 100:2008. Guide to the expression of uncertainty of measurement

JCGM 100:2005. Supplement 1 Propagation of distributions usign a Monte Carlo method.

EURACHEM/CITAC Guide CG 4. Quantifying Uncertainty in Analytical Measurement

See Also

summary.uncertainty, plot

Examples

# create an uncertainty budget
cor.mat <- matrix(c(1, -0.7, -0.7, 1), 2, 2)

u.budget <- uncertaintyBudget(x = list(name = c("x0", "x1"), 
	mean = c(10, 20), u = c(1, 5), unit = c("kg", "kg"), dof = c(10, 10),
	description = c("measurand mass", "sample mass"),
	label = c("x[0]", "x[1]"), type = c("A", "A"), distribution = c("normal", "normal")), 
	y = cor.mat)

# estimate the measurand uncertainty using an uncertainty budget,
# a measurand definition and a selected estimating method.
GFO.res <- uncertainty(x = u.budget, 
y = list(measurand_name = "ratio.GFO",
         measurand_label = "ratio[GFO]", 
         measurand_model = "x0/x1",
	 measurand_description = "ratio of masses at 20 degrees celsius",
         method = "GFO", alpha = 0.05))

# create an uncertainty summary object
GFO.sum <- summary(GFO.res)

# display the chart
## Not run: plot(GFO.sum)


uncertainty documentation built on Dec. 15, 2025, 1:06 a.m.