# 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), units = c("kg", "kg"), u = c(1, 5), dof = c(10, 10),
label = c("x[0]", "x[1]"), distribution = c("normal", "normal"),
description = c("measurand mass", "sample mass"),
type = c("A", "A")),
y = cor.mat)
u.budget
# 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",
method = "GFO", alpha = 0.05))
GFO.res
Run the code above in your browser using DataLab