Creates an uncertainty estimation object using a measurand model and an uncertainty budget object
an uncertainty budget object
a list with the measurand description and selected estimation method, the measurand description includes: measurand_name, measurand_model, measurand_label, alpha (significance level), method and method parameters.
the valid methods are: GFO, GSO, MC, Kragten.
currently the only method parameter implemented is the number of simulated samples (B) for the method MC.
Creates an uncertainty estimation object. Uses an uncertainty budget object to estimate the expected value and uncertainty of a measurand by applying a selected estimation method.
An uncertainty estimation object with the structure:
method selected estimating method,
call current call invocation,
uncertaintyBudget an uncertainty budget object,
measurand name, label, model describing the measurand,
mean the estimated mean,
sd the estimated standard deviation,
u the estimated standard uncertainty,
alpha the significante level used in the estimation,
dof the estimated degrees of freedom,
U the estimated expanded uncertainty,
lcl the lower confidence interval,
ucl the upper confidence interval,
variables a vector with the input quantities,
contribution a vector with the uncertainty contributions,
cor.contribution the uncertainty contribution due to overall correlation,
partial a vector of the partial derivatives of the measurand.model with respect to each input quantity,
coeff a vector of the sensibility coefficients for each input quantity.
Maintainer: H. Gasca-Aragon <[email protected]>
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
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
# 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), dof=c(10,10), label=c("x", "x"), distribution=c("normal","normal")), 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", method="GFO", alpha=0.05)) GFO.res
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.