Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
# Render html vignetes by using devtools::document(roclets = "vignette")
# Render also pdf vignetes by using rmarkdown::render("vignettes/masscor.Rmd", "all")
## ---- eval = FALSE-----------------------------------------------------------
# install.packages("masscor")
## ---- eval = FALSE-----------------------------------------------------------
# devtools::install_github("Crparedes/masscor", build_vignettes = TRUE)
## ----setup--------------------------------------------------------------------
library(masscor)
## ----calibCert1---------------------------------------------------------------
Bal.1.Lab.317 <- calibCert(balanceID = '(Brand) Analytic balance lab 317',
d = 0.1, d.units = 'mg',
indError = data.frame(nominal = c(5, 20), # grams
error = c(-0.2, -0.1), # miligrams
uncert = c(0.200, 0.201)), # miligrams
indError.units = c('g', 'mg', 'mg'),
rep = c(20, 0.02), rep.units = c('g', 'mg'),
eccen = c(100, 0.1), eccen.units = c('g','mg'))
## ----calibCert2---------------------------------------------------------------
print(Bal.1.Lab.317)
plot(Bal.1.Lab.317)
## ----calibCert3, eval = FALSE------------------------------------------------
# save(Bal.1.Lab.317, file = 'certificateBalance.1_Lab.317.RData')
# load(file = 'certificateBalance.1_Lab.317.RData')
## ----calibCert4---------------------------------------------------------------
data(MT.XPE.204)
print(MT.XPE.204, complete = TRUE)
plot(MT.XPE.204)
## ----convMass1----------------------------------------------------------------
data(MT.XPE.204)
convMass(calibCert = MT.XPE.204, reading = 211.1342, units = 'g')
## ----convMass2----------------------------------------------------------------
convMass(calibCert = MT.XPE.204, reading = 211.1342, units = 'g', rho = 7.113,
rho_air = airDensity(Temp = 21.1, p = 990, h = 46.6,
unitsENV = c("deg.C", "hPa", "%")))
## ----convMass1a---------------------------------------------------------------
(u_err <- uncertErrorCorr(calibCert = MT.XPE.204, reading = 211.1342, units = 'g'))
## Result rounded to three significant figures:
signif(u_err, 3)
## ----convMass1b---------------------------------------------------------------
(u_read <- uncertReading(calibCert = MT.XPE.204, reading = 211.1342, units = 'g'))
## ----convMass1c---------------------------------------------------------------
(u_w <- sqrt(u_err ^ 2 + u_read ^ 2))
(u_w <- uncertConvMass(calibCert = MT.XPE.204, reading = 211.1342, units = 'g'))
## ----MABC1--------------------------------------------------------------------
MABC(rho = 2.115, rho_w = 8, rho_air = 0.001199314) # g/cm^3
## ----airDens------------------------------------------------------------------
Temp <- 21.5 # Celcius degrees
p <- 751 # Hectopascals
h <- 58.1 # Percentaje
unitsENV = c("deg.C", "hPa", "%")
(rho_a.1 <- airDensity(Temp = Temp, p = p, h = h, unitsENV = unitsENV, model = 'CIMP2007'))
(rho_a.2 <- airDensity(Temp = Temp, p = p, h = h, unitsENV = unitsENV, model = 'CIMP.approx'))
(rho_a.3 <- airDensity(Temp = Temp, p = p, h = h, unitsENV = unitsENV, model = 'Jones1978'))
(rho_a.4 <- airDensityHASL(HASL = 2640)) # Height above sea level for Bogotá-Colombia
## ----uncertAirDens------------------------------------------------------------
(u_rho_a.1 <- uncertAirDensity(Temp = Temp, p = p, h = h, unitsENV = unitsENV,
u_Temp = 2.9, u_p = 10.10, u_h = 11.3, model = 'CIMP2007',
plot = TRUE))
## ----MABC2--------------------------------------------------------------------
uncertMABC(rho = 2.115, rho_w = 8, rho_air = rho_a.1,
u_rho = 0.005/sqrt(3), u_rho_w = 0.060, u_rho_air = u_rho_a.1,
plot = TRUE)
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.