-
-
Notifications
You must be signed in to change notification settings - Fork 0
/
calibration.Rd
161 lines (135 loc) · 5.35 KB
/
calibration.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/calibration.R
\name{calibration}
\alias{calibration}
\alias{print.calibration}
\alias{summary.calibration}
\alias{plot.calibration}
\alias{lod}
\alias{lod.calibration}
\alias{loq}
\alias{loq.calibration}
\title{Analytical calibration functions}
\usage{
calibration(
formula,
data = NULL,
blanks = NULL,
weights = NULL,
model = "lm",
check_assumptions = TRUE,
...
)
\method{print}{calibration}(x, ...)
\method{summary}{calibration}(object, ...)
\method{plot}{calibration}(x, interval = "conf", level = 0.95, ...)
lod(object, ...)
\method{lod}{calibration}(object, blanks = NULL, alpha = 0.01, level = 0.05, ...)
loq(object, ...)
\method{loq}{calibration}(
object,
blanks = NULL,
alpha = 0.01,
k = 3,
level = 0.05,
maxiter = 10,
...
)
}
\arguments{
\item{formula}{model formula providing the recorded signal intensities with
respect to the nominal analyte concentrations in the form of
\code{signal ~ concentration} or \code{signal ~ concentration - 1}; model
formulas are currently restricted to those forms, however, the possibility
to use \code{log} or \code{sqrt} transformed data will be implemented in the
future.}
\item{data}{an optional data frame containing the variables in the model.}
\item{blanks}{a vector of numeric blank values overriding those automatically
retrieved from calibration data.}
\item{weights}{an optional character string containing one or more model
variables, for example, in the form of "\code{1/concentration^0.5}" or
"\code{1/signal}" which is internally converted to a numeric vector and
passed to the fitting process of the selected model.}
\item{model}{model class to be used for fitting; currently,
\code{\link[stats]{lm}()} and \code{\link[MASS]{rlm}()} are supported.}
\item{check_assumptions}{automatically check for normality and
homoscedasticity of model residuals using \code{\link[stats]{shapiro.test}()}
and \code{\link[lmtest]{bptest}()}, respectively; only executed if
\code{weights == NULL}.}
\item{\dots}{further arguments passed to the submethod, namely the
respective model environment such as \code{lm}(), \code{plot}(), or
\code{print}().}
\item{x, object}{an object of class '\code{calibration}' with a model formula
as shown above.}
\item{interval}{Type of interval calculation (can be abbreviated); see
\code{\link[stats]{predict}()} for details.}
\item{level}{tolerance/confidence level; see \code{\link[stats]{predict}()}
and \code{\link[stats]{confint}()} for details.}
\item{alpha}{error tolerance for the detection limit (critical value).}
\item{k}{relative uncertainty for the limit of quantification
(\code{1/beta}).}
\item{maxiter}{a positive integer specifying the maximum number of iterations
to calculate the LOQ.}
}
\value{
\code{calibration} returns an object of \code{\link[base]{class}}
'\code{calibration}'. \code{print()} calls the function parameters together
with the respective LOD and LOQ. \code{plot()} plots the respective
calibration curve together with the measurement values. \code{summary()} may
be used to retrieve the summary of the underlying model.
}
\description{
Defines a \code{calibration} object for the calculation of concentrations
from measurement signals including estimations for the limit of detection
(LOD) and limit of quantification (LOQ) in accordance with DIN 32645 (2008).
The LOD is defined as the lowest quantity of a substance that can be
distinguished from the absence of that substance (blank value) within a given
confidence level (\code{alpha}). The LOQ is defined as the lowest quantity of
a substance that can be quantified/distinguished from another sample given
with respect to a defined confidence level (\code{k}).
}
\details{
If the \code{data} supplied to \code{calibration} contain more than one blank
value, namely measurements with a nominal concentration of or close to zero,
the LOD and LOQ are calculated from the deviation of the blank samples. This
method is called "blank method" according to DIN 32645 (2008) and supposed
to be more accurate than the so-called "calibration method" which will be
used for the estimation of LOD and LOQ when \code{data} does not contain zero
concentration measurements.
}
\examples{
data(din32645)
din <- calibration(Area ~ Conc, data = din32645)
din
plot(din)
summary(din)
lod(din)
loq(din)
}
\references{
Almeida, A. M. D., Castel-Branco, M. M., & Falcao, A. C. (2002). Linear
regression for calibration lines revisited: weighting schemes for
bioanalytical methods. \emph{Journal of Chromatography B}, \bold{774}(2),
215-222. \doi{10.1016/S1570-0232(02)00244-1}.
Currie, L.A. (1999). Nomenclature in evaluation of analytical methods
including detection and quantification capabilities: (IUPAC Recommendations
1995). \emph{Analytica Chimica Acta} \bold{391}, 105-126.
DIN 32645 (2008). \emph{Chemical analysis - Decision limit, detection limit
and determination limit under repeatability conditions - Terms, methods,
evaluation}. Technical standard. Deutsches Institut für Normung, Berlin.
Massart, D.L., Vandeginste, B.G., Buydens, L.M.C., Lewi, P.J., &
Smeyers-Verbeke, J. (1997). \emph{Handbook of chemometrics and qualimetrics:
Part A}. Elsevier Science Inc.
}
\seealso{
Other calibration:
\code{\link{din32645}},
\code{\link{icp}},
\code{\link{matrix_effect}()},
\code{\link{neitzel2003}},
\code{\link{weight_select}()}
}
\author{
Zacharias Steinmetz
}
\concept{calibration}