Title:	Mass Measurement Corrections
Version:	0.0.7.1
Date:	2021-09-09
Description:	Mass measurement corrections and uncertainties using calibration data, as recommended by EURAMET's guideline No. 18 (2015) ISBN:978-3-942992-40-4 . The package provides classes, functions, and methods for storing information contained in calibration certificates and converting balance readings to both conventional mass and real mass. For the latter, the Magnitude of the Air Buoyancy Correction factor employs models (such as the CIMP-2007 formula revised by Picard, Davis, Gläser, and Fujii (2008) <doi:10.1088/0026-1394/45/2/004>) to estimate the local air density using measured environmental conditions.
License:	GPL (≥ 3)
Encoding:	UTF-8
LazyData:	true
RoxygenNote:	7.1.1
Suggests:	testthat (≥ 3.0.0), knitr, rmarkdown
Config/testthat/edition:	3
Imports:	metRology
Depends:	R (≥ 3.5.0)
VignetteBuilder:	knitr
NeedsCompilation:	no
Packaged:	2021-09-09 15:43:26 UTC; cris
Author:	Cristhian Paredes [aut, cre]
Maintainer:	Cristhian Paredes <craparedesca@unal.edu.co>
Repository:	CRAN
Date/Publication:	2021-09-13 08:20:02 UTC

masscor: Mass Measurement Corrections and Uncertainties.

Description

The R package masscor provides functions, classes and methods to support mass measurements using non automatic balances as described in EURAMET's Calibration Guide No. 18 (2015). The new classes are objects that can store the calibration information for balances and mass standards. Those objects can be used to convert balance readings to both conventional mass and mass, and to perform routine balance verification (by using the normalized error function). Air buoyancy correction factors are calculated using local air density that can be calculated using environmental conditions and applying one of several models available in the package. The uncertainty of (corrected) mass measurements can also be evaluated allowing us to further assess the suitability of given mass measurement.

masscor functions

This package uses list objects of class 'calibCert' to store information of balance calibration certificates. The functions use the information of this object to convert balance reading indications to conventional mass and calculate mass uncertainties.

Several models for calculating air density are included and this information can be used to calculate the Magnitude of the Air Buoyancy Correction factor (MABC). Uncertainties calculations are made using Gauss Approximation according to the Guide to the Expression of Uncertainty in Measurement (GUM) implemented in R by the package metRology (Ellison, 2018).

Author(s)

Cristhian Paredes, craparedesca@unal.edu.co

References

EURAMET, Calibration Guide No. 18. 2015. Guidelines on the Calibration of Non-Automatic Weighing Instruments. https://www.euramet.org/Media/docs/Publications/calguides/I-CAL-GUI-018_Calibration_Guide_No._18_web.pdf.

Picard, A; Davis, R S; Gläser, M; Fujii, K (2008). Revised formula for the density of moist air (CIPM-2007). Metrologia, 45(2), 149–155. doi:10.1088/0026-1394/45/2/004

Harris, G. (2019). Selected Laboratory and Measurement Practices and Procedures to Support Basic Mass Calibrations. SOP 2 - Recommended Standard Operating Procedure for Applying Air Buoyancy Corrections. National Institute of Standards and Technology (NIST). doi:10.6028/NIST.IR.6969-2019

BIMP JCGM (2008) Evaluation of measurement data — Guide to the expression of uncertainty in measurement.

Stephen L R Ellison. (2018). metRology: Support for Metrological Applications. R package version 0.9-28-1. https://CRAN.R-project.org/package=metRology

Calibration data of a weights kit in the interval 1 mg to 1 kg

Description

An object of class "massStandardKit", with calibration information of mass standards with nominal masses in the interval 1 mg to 1 kg, performed by the Mass Laboratory at the Instituto Nacional de Metrologia de Colombia (2020-08-12).

Usage

Box.E2.MS.Kit

Format

A list with 25 objects of class massStandard()

Source

Mass Laboratory - Instituto Nacional de Metrologia de Colombia. CALIBRATION CERTIFICATE No. 4687.

Calibration data of a E2 class mass standard of nominal mass 20 g.

Description

An object of class "massStandard", with calibration information of a mass standards with nominal mass of 20 g, performed by the Mass Laboratory at the Instituto Nacional de Metrologia de Colombia (2020-08-12).

Usage

E2.MS.20g

Format

A list with 25 objects of class massStandard()

Source

Mass Laboratory - Instituto Nacional de Metrologia de Colombia. CALIBRATION CERTIFICATE No. 4687.

Magnitude of the Air Buoyancy Correction

Description

Calculates the Magnitude of the Air Buoyancy Correction (MABC). If no parameters are provided the function returns MABC for weighing water at standard conditions.

Usage

MABC(rho = 0.997, rho_w = 8, rho_air = airDensity())

Arguments

Details

Comparing masses (weighing) in air produces results that are influenced by the objects densities due to their buoyancy in air. This air buoyancy effects are usually small but must be taken in account when high accuracy is required. The effect can be corrected by using the densities of the object, the mass standard and the air filling the room where the measurement process takes place (Harris, 2019).

The uncertainty associated to MABC can be calculated by the function uncertMABC().

Value

Numeric value of the Magnitude of the Air Buoyancy Correction factor.

References

Harris, G. (2019). Selected Laboratory and Measurement Practices and Procedures to Support Basic Mass Calibrations. SOP 2 - Recommended Standard Operating Procedure for Applying Air Buoyancy Corrections. National Institute of Standards and Technology (NIST). doi:10.6028/NIST.IR.6969-2019

See Also

uncertMABC(), airDensity(), uncertMABC()

Examples

## Magnitude of the air buoyancy correction for some materials:
# Water
MABC()
# Zinc metal
MABC(rho = 7.133)
# Copper metal
MABC(rho = 8.96)

Calibration data for a balance Mettler Toledo XP 2002

Description

An object of class calibCert() with the calibration information of a balance Mettler Toledo XP 2002, performed by the Mass Laboratory at the Instituto Nacional de Metrologia de Colombia (2021-03-12).

Usage

MT.XP.2002

Format

An object of class calibCert()

Source

Mass Laboratory - Instituto Nacional de Metrologia de Colombia. CALIBRATION CERTIFICATE No. 5142.

Calibration data for a balance Mettler Toledo XPE 204

Description

An object of class calibCert() with the calibration information of a balance Mettler Toledo XPE 204, performed by the Mass Laboratory at the Instituto Nacional de Metrologia de Colombia (2021-03-18).

Usage

MT.XPE.204

Format

An object of class calibCert()

Source

Mass Laboratory - Instituto Nacional de Metrologia de Colombia. CALIBRATION CERTIFICATE No. 5143.

Models for calculating air density based on environmental conditions

Description

The function uses environmental conditions information (barometric pressure, temperature and relative humidity.) to calculate a local air density value. If no parameter is defined, the air density at 20^oC, 1013.25 hPa and 50% relative humidity is returned. The air density value value can later be used to calculate the Magnitude of Air Buoyancy Correction (MABC()). The uncertainty of the air density value can be calculated using the function uncertAirDensity().

Usage

airDensity(Temp = 20, p = 1013.25, h = 50, unitsENV = c("deg.C", "hPa",
  "%"), x_CO2 = 4e-04, model = "CIMP2007")

Arguments

Details

Local air density can be estimated using one of several methods. The most complete approach is the CIMP complete formula (the default, method = 'CIMP2007') as described in Picard et al (2008). The CIMP approximated exponential formula (method = 'CIMP.approx') and the method reported by Jones, (method = 'Jones1978') (Harris, 2019) are also included.

Value

Numeric value of air density in g~cm^{-3}, according to chosen model.

unitsENV

Temperature units (Temp) can be either 'deg.C' (for Celsius degrees) or 'K'. Pressure units (p) can be any of 'mmHg', 'Pa', 'hPa' or 'kPa'. Relative humidity (h) can be expressed as fraction ('frac') or as percentage ('%'). A typical arrangement for the parameter unitsENV would be c('deg.C', 'hPa', '%').

References

Picard, A; Davis, R S; Gläser, M; Fujii, K (2008). Revised formula for the density of moist air (CIPM-2007). Metrologia, 45(2), 149–155. doi:10.1088/0026-1394/45/2/004

Harris, G. (2019). Selected Laboratory and Measurement Practices and Procedures to Support Basic Mass Calibrations. SOP 2 - Recommended Standard Operating Procedure for Applying Air Buoyancy Corrections. National Institute of Standards and Technology (NIST). doi:10.6028/NIST.IR.6969-2019

Preguntar a andres referencia de la f'ormula simplificada exponencial

See Also

MABC() to calculate the Magnitude of Air Buoyancy Correction and uncertAirDensity() to estimate the uncertainty of the calculated air density.

Examples

airDensity(Temp = 23.4, p = 612.3, h = 23,
           unitsENV = c('deg.C', 'mmHg', '%')) # [g/cm^3]
airDensity(Temp = 23.4, p = 612.3, h = 23,
           unitsENV = c('deg.C', 'mmHg', '%'), model = 'CIMP.approx') # [g/cm^3]
airDensity(Temp = 23.4, p = 612.3, h = 23,
           unitsENV = c('deg.C', 'mmHg', '%'), model = 'Jones1978')   # [g/cm^3]

Air density estimation using only height above sea level.

Description

Calculates the approximated density of local air using a model that relies on height above sea level (HASL) information. More accurate alternatives are found in airDensity() but those require data form environmental conditions (temperature, barometric pressure and relative humidity).

Usage

airDensityHASL(HASL)

Arguments

Value

Numeric value of an approximated air density in g cm^{-3}.

See Also

airDensity() for better models to predict air density.

Examples

airDensityHASL(HASL = 0)    # [g/cm^3]
airDensityHASL(HASL = 1600) # [g/cm^3]

Information of balance calibration certificate

Description

Creates an object of class calibCert that contains the information of a balance calibration certificate. The object can later be used to correct mass readings and calculate mass uncertainties. Mandatory arguments for this function are the balance division scale (d), the results of the indication error test (indError), the results of repeatability test (rep), and the results of the eccentricity test (eccen).

Usage

calibCert(balanceID = "BalanceID", serial = NULL, certificate = NULL, d,
  d.units = "mg", indError, indError.units = c("g", "mg", "mg"),
  expanded = TRUE, k = 2, rep, rep.units = c("g", "mg"), eccen,
  eccen.units = "mg", classSTD = NULL, traceability = NULL,
  Temp = NULL, p = NULL, h = NULL, unitsENV = c("deg.C", "hPa", "%"),
  institution = NULL, accreditation = NULL, date = NULL,
  add.info = NULL)

Arguments

Details

The units of d, indError, rep and eccen shall be provided to the arguments d.units, indError.units, rep.units and eccen.units, respectively. The units can be any multiple or subdivision of the SI unit for mass, the kilogram. The greek letter \mu used to represent a millionth part, is replaced by the vocal u. Remember that both R and the SI prefixes are case sensitive.

Value

Object of class calibCert with information of the calibration certificate for a balance.

unitsENV

Temperature units (Temp) can be either 'deg.C' (for Celsius degrees) or 'K'. Pressure units (p) can be any of 'mmHg', 'Pa', 'hPa' or 'kPa'. Relative humidity (h) can be expressed as fraction ('frac') or as percentage ('%'). A typical arrangement for the parameter unitsENV would be c('deg.C', 'hPa', '%').

See Also

S3 methods print.calibCert() and plot.calibCert() are available. See convertMassUnitsSI() for information about mass units.

Examples

massSTD  <- c(0.01, 0.5, 1, 10, 20, 50, 100, 120, 150, 200, 220)  ## [g]
indError <- c(0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -0.2, -0.2) ## [mg]
uncert   <- c(0.1, 0.1, 0.1, 0.1, 0.1, 0.2, 0.2, 0.3, 0.3, 0.4, 0.5) ## [mg]
d <- 0.1 ## [mg]

Balance.D1 <- calibCert(balanceID = 'MT XPE 204', serial = 'B403223982',
                        d = d, d.units = 'mg',
                        indError = data.frame(massSTD, indError, uncert),
                        indError.units = c('g', 'mg', 'mg'),
                        rep = data.frame(load = c(0.1, 100, 220),
                                         sd = c(0.00, 0.04, 0.03)),
                        rep.units = c('g', 'mg'),
                        eccen = c(100, 0.1), eccen.units = c('g', 'mg'),
                        Temp = c(17.4, 17.9), ## [deg.C]
                        p = c(750.4, 751.0), ## [hPa]
                        h = c(70.5, 71.4), ## [%]
                        unitsENV = c('deg.C', 'hPa', '%'),
                        institution = 'Instituto Nacional de Metrologia de Colombia',
                        date = '2021-03-18')
print(Balance.D1)

Corrects a balance reading using balance calibration data.

Description

Given a balance reading indication and the calibration information of the balance, the function interpolates error correction for the reading using the errors of indication for the two closest calibration points. The output is generally the mass measurement result under the conditions of calibration. If densities from the object and the local air are provided the conventional mass of the object can be calculated. See Details.

Usage

convMass(calibCert, reading, units = NULL, rho = NULL, rho_air = NULL)

Arguments

Details

The conventional mass value of a body is equal to the mass m_c of a mass standard that balances this body under conventionally chosen conditions: at a temperature t_{ref} = 20^oC, with mass standards of density \rho_c=8000 kg m^{-3}, in normal air of density \rho_0=1.2 kg m^{-3} (OIML, 2004).

Value

Numeric value of conventional mass.

References

OIML, (2004). ORGANISATION INTERNATIONALE DE MÉTROLOGIE LÉGALE. International Document D 28: Conventional value of the result of weighing in air.

Harris, G. (2019). Selected Laboratory and Measurement Practices and Procedures to Support Basic Mass Calibrations. SOP 2 - Recommended Standard Operating Procedure for Applying Air Buoyancy Corrections. National Institute of Standards and Technology (NIST). doi:10.6028/NIST.IR.6969-2019

See Also

uncertConvMass()

Examples

data(minimalCert)
convMass(reading = 12.4835, calibCert = minimalCert)

Conversion between mass units of The International System of Units

Description

Mass values are converted from a SI unit to another SI unit according to the SI prefixes as shown in BIMP (2019). The greek letter \mu is replaced by the vocal u.

Usage

convertMassUnitsSI(value, from, to)

Arguments

Value

Numeric vector of mass values converted from a SI unit to another mass unit.

References

BIMP, 2019. Bureau International des Poids et Mesures. Brochure of The International System of Units. 9th Edition.

Examples

convertMassUnitsSI(value = c(0.2, 0.4), from = 'mg', to = 'g')

Creates an object of class "massStandard".

Description

The object of class "massStandard" contains the calibration information of a mass standard that is used in routine balance verification (e g. to calculate normalized error. See normalizedError()). A version to store information of several mass standards that belong to the same kit is massStandardKit().

Usage

massStandard(nominal, convMassCor, uncert, units = c("g", "mg", "mg"),
  serial = NULL, manufacturer = NULL, class = NULL, certificate = NULL,
  traceability = NULL, Temp = NULL, p = NULL, h = NULL,
  unitsENV = c("deg.C", "hPa", "%"), expanded = TRUE, k = 2, rho = 8,
  u_rho = 0.06, unitsrho = "g/cm^3", institution = NULL, date = NULL,
  add.info = NULL, partofakit = FALSE)

Arguments

Value

Object of class "massStandard" with the information of a calibrated mass standard.

unitsENV

Temperature units (Temp) can be either 'deg.C' (for Celsius degrees) or 'K'. Pressure units (p) can be any of 'mmHg', 'Pa', 'hPa' or 'kPa'. Relative humidity (h) can be expressed as fraction ('frac') or as percentage ('%'). A typical arrangement for the parameter unitsENV would be c('deg.C', 'hPa', '%').

References

OIML, (2004). ORGANISATION INTERNATIONALE DE MÉTROLOGIE LÉGALE. International Document OIML R 111: Weights of classes E1, E2, F1, F2, M1, M1–2, M2, M2–3 and M3.

See Also

massStandardKit(), normalizedError()

Examples

singleMS.E2.10g <- massStandard(nominal = 10, convMassCor = 0.001,
                                uncert = 0.1,
                                units = c('g', 'mg', 'mg'))
print(singleMS.E2.10g)

Creates an object of class "massStandardKit".

Description

The object of class "massStandardKit" is a wrapper for several objects of class "massStandard". The object of class "massStandard" contains the calibration information of a mass standard that is used in routine balance verification (e g. to calculate normalized error. See normalizedError()). When several mass standards are part of a kit their information can be in conveniently be stored together in a "massStandardKit" class object.

Usage

massStandardKit(nominal, convMassCor, uncert, units = c("g", "mg", "mg"),
  serial = NULL, manufacturer = NULL, class = NULL, certificate = NULL,
  traceability = NULL, Temp = NULL, p = NULL, h = NULL,
  unitsENV = c("deg.C", "hPa", "%"), expanded = TRUE, k = 2,
  rho = NULL, u_rho = NULL, unitsrho = "g/cm^3", institution = NULL,
  date = NULL, add.info = NULL)

Arguments

Details

When two mass standards of equal nominal mass are present in the same kit (typically those of nominal mass in the form 2\times10^n), it is necessary to make a distinction between them because their real mass are not likely to be exactly the same. Physically this distinction is achieved by marking one of them, for example, with a dot in the head of the knot weights, or by bending the final part of the lifted extreme in wire weights. To differentiate those duplicated mass standards in the "massStandardKit" class object, its nominal mass value must be entered as a character including an asterisk after the value (ie. '200*' instead of just entering 200). The function returns error if the kit contains duplicated mass standards and no differentiation is indicated.

Value

Object of class "massStandard" with calibration information of a mass standards kit.

unitsENV

Temperature units (Temp) can be either 'deg.C' (for Celsius degrees) or 'K'. Pressure units (p) can be any of 'mmHg', 'Pa', 'hPa' or 'kPa'. Relative humidity (h) can be expressed as fraction ('frac') or as percentage ('%'). A typical arrangement for the parameter unitsENV would be c('deg.C', 'hPa', '%').

See Also

massStandard(), normalizedError()

Examples

nominal     <- c(1000, 500, 200, '200*', 100)   # [g]
convMassCor <- c(0.0, -0.03, 0.03, 0.06, 0.00)  # [mg]
uncert      <- c(0.50, 0.25, 0.10, 0.10, 0.05)  # [mg]
units       <- c('g', 'mg', 'mg')

rho         <- c(8012.217, 8008.640, 8011.126, 8010.722, 8010.935)# [kg/m^3]
u_rho       <- c(0.096, 0.090, 0.160, 0.160, 0.321)# [kg/m^3]
unitsrho    <- 'kg/m^3'

MS.Kit1 <- massStandardKit(nominal = nominal, convMassCor = convMassCor, uncert = uncert,
                           units = units, rho = rho, u_rho = u_rho, unitsrho = unitsrho)
print(MS.Kit1)

Minimal calibCert() object

Description

An object of class calibCert() with the minimal calibration information of a hypothetical balance calibrated using only two points.

Usage

minimalCert

Format

An object of class calibCert()

Calculates normalized in balance verification using a mass standard

Description

Calculates normalized in balance verification using a mass standard

Usage

normalizedError(reading, standard, calibCert, u_massStandard = NULL)

Arguments

Value

Numeric value of normalized error for balance verification using a mass standard.

See Also

massStandard().

Examples

data(E2.MS.20g)
data(MT.XPE.204)
normalizedError(reading = 20.0000, standard = E2.MS.20g, calibCert = MT.XPE.204)

S3 method for plotting objects of class "calibCert"

Description

The function plots the indication error or the conventional mass correction for a balance whose calibration data is in a object of class "calibCert".

Usage

## S3 method for class 'calibCert'
plot(x, error = TRUE, y0line = TRUE, ...)

Arguments

Value

A base plot with calibration data of indication error or correction.

See Also

calibCert(), print.calibCert()

Examples

data(MT.XPE.204)
plot(MT.XPE.204)

S3 method for printing objects of class "calibCert"

Description

The function prints objects of class "calibCert".

Usage

## S3 method for class 'calibCert'
print(x, complete = FALSE, nudeCertificate = FALSE, ...)

Arguments

Value

No return value, called for side effects.

See Also

calibCert(), plot.calibCert()

Examples

data(MT.XPE.204)
print(MT.XPE.204)

S3 method for printing objects of class "massStandard"

Description

The function prints objects of class "massStandard".

Usage

## S3 method for class 'massStandard'
print(x, minimal = TRUE, description = TRUE,
  institution = TRUE, density = FALSE, envConditions = TRUE,
  addInfo = TRUE, ...)

Arguments

Value

No return value, called for side effects.

See Also

massStandard(), print.massStandardKit()

Examples

data(E2.MS.20g)
print(E2.MS.20g)
print(E2.MS.20g, minimal = FALSE)

S3 method for printing objects of class "massStandardKit"

Description

The function prints objects of class "massStandardKit".

Usage

## S3 method for class 'massStandardKit'
print(x, minimal = FALSE, description = TRUE,
  institution = TRUE, density = FALSE, envConditions = TRUE,
  addInfo = TRUE, ...)

Arguments

Value

No return value, called for side effects.

See Also

massStandardKit(), print.massStandard()

Examples

data(Box.E2.MS.Kit)
print(Box.E2.MS.Kit, minimal = TRUE)
# We can print individual information of a single mass standard:
print(Box.E2.MS.Kit[['20']])

Uncertainties in air density calculations

Description

Propagates the uncertainty of environmental conditions measurements to estimated values of air densities calculated using the function airDensity() with models 'CIMP2007' y 'CIMP.approx' function. Uncertainty arising from to the chosen model itself is considered. Calculations are made according to the Guide to the Guide to the expression of uncertainty in measurement (GUM, JCGM, 2008) as implemented by the package metRology (Ellison, 2018).

Usage

uncertAirDensity(model = "CIMP2007", Temp = 20, p = 1013.25, h = 50,
  u_Temp = 2.9, u_p = 10.1, u_h = 11.3, unitsENV = c("deg.C", "hPa",
  "%"), plot = TRUE, printRelSD = TRUE)

Arguments

Value

A numeric value of standard uncertainty of calculated air density in g~cm^{-3}.

unitsENV

Temperature units (Temp) can be either 'deg.C' (for Celsius degrees) or 'K'. Pressure units (p) can be any of 'mmHg', 'Pa', 'hPa' or 'kPa'. Relative humidity (h) can be expressed as fraction ('frac') or as percentage ('%'). A typical arrangement for the parameter unitsENV would be c('deg.C', 'hPa', '%').

References

Picard, A; Davis, R S; Gläser, M; Fujii, K (2008). Revised formula for the density of moist air (CIPM-2007). Metrologia, 45(2), 149–155. doi:10.1088/0026-1394/45/2/004

Harris, G. (2019). Selected Laboratory and Measurement Practices and Procedures to Support Basic Mass Calibrations. SOP 2 - Recommended Standard Operating Procedure for Applying Air Buoyancy Corrections. National Institute of Standards and Technology (NIST). doi:10.6028/NIST.IR.6969-2019

BIMP JCGM (2008) Evaluation of measurement data — Guide to the expression of uncertainty in measurement.

Stephen L R Ellison. (2018). metRology: Support for Metrological Applications. R package version 0.9-28-1. https://CRAN.R-project.org/package=metRology

See Also

airDensity()

Examples

 uncertAirDensity(model = 'CIMP2007',
                  Temp = 20, p = 1013.25, h = 50,
                  u_Temp = 0.29, u_p = 1.01, u_h = 11.3)

Uncertainty in conventional mass value

Description

The function combines the uncertainty of the conventional mass correction (as obtained by uncertErrorCorr()) and the uncertainty in the balance reading (as obtained by uncertReading()), to produce the uncertainty of a conventional mass value.

Usage

uncertConvMass(calibCert, reading, units, sd, sd.units, d, d.units)

Arguments

Value

A numeric value of uncertainty for a conventional mass value.

See Also

convMass(), uncertReading(), uncertErrorCorr()

Examples

data(minimalCert)
uncertConvMass(reading = 12.4835, calibCert = minimalCert)

Uncertainty due to mass correction using calibration certificate

Description

Given a balance reading indication and the calibration information of the balance, the function uses the conventional mass correction uncertainties of the two closest calibration points to the balance reading to interpolate the uncertainty due to the conventional mass correction.

Usage

uncertErrorCorr(calibCert, reading, units = NULL)

Arguments

Value

A numeric value of uncertainty for a conventional mass correction.

See Also

convMass(), uncertReading(), uncertConvMass()

Examples

  data(minimalCert)
  uncertErrorCorr(reading = 12.4835, calibCert = minimalCert)

Uncertainty of the Magnitude of the Air Buoyancy Correction factor

Description

Propagates density uncertainties in the calculation of the Magnitude of Air Buoyancy Correction (See MABC()).

Usage

uncertMABC(rho = 0.998, rho_w = 8, rho_air = NULL, u_rho = 1e-04,
  u_rho_w = 0.006, u_rho_air = NULL, plot = TRUE, printRelSD = TRUE)

Arguments

Details

Calculations are made according to the Guide to the Guide to the expression of uncertainty in measurement (GUM, JCGM, 2008) as implemented by the package metRology (Ellison, 2018). If air density and associated uncertainty are not provided the default output values of the functions airDensity() and uncertAirDensity(), respectively, are used.

Value

Numeric value of uncertainty for the Magnitude of the Air Buoyancy Correction factor.

References

BIMP JCGM (2008) Evaluation of measurement data — Guide to the expression of uncertainty in measurement.

Andrej-Nikolai Spiess (2018). propagate: Propagation of Uncertainty. R package version 1.0-6. https://CRAN.R-project.org/package=propagate

Uncertainty of balance readings

Description

Uncertainty in a given balance reading considering the effects of rounding error, lack of repeatability, eccentricity and balance taring.

Usage

uncertReading(calibCert, reading, units = NULL, sd = NULL,
  sd.units = NULL, d = NULL, d.units = NULL)

Arguments

Value

A numeric value of uncertainty for a balance reading.

See Also

uncertErrorCorr(), uncertConvMass()

Examples

data(minimalCert)
uncertReading(calibCert = minimalCert, reading = 12.4835)
uncertReading(calibCert = minimalCert, reading = 12.484, d = 1, d.units = 'mg')