Limit of Detection A Closer Look at The IUPAC Definition
Limit of Detection A Closer Look at The IUPAC Definition
ft • ft % m # ft ft « « * ft ft ft » m 9 ft ft ft ft ft ft # ft ft ft ft ft ft
ft © ft ft ft ft ft ft ft ft ft ft ft ft • ft ft ft ft ft ft 9 ft ft ft 0
ft ft ft * ft ft ft « ft ft ft ft ft ft # • ft ft ft ft ft ft ft ft • ft ft ft ft ft «
ft ft ft ft # ft ft ft ft ft ft ft ft ft ft ft ft # ft ft ft ft ft # ft
<2 « ft ® # ft ft ft e? # ft ft ft ft * ft ft ft ft ft ft ft ft ft « ft ft ft ft ft
ft ft ft ft ft ft ft * ft ft ft » • ft ft ft ft ft ft ft jit ft ft ft ft ft ft ft
#v
>
ft
ft
ft
•
ft
•
$
m W * I*
T-
ft
ft
ft
ft
ft
# ft
ft
ft
•
ft
ft m
ftw
ft
ft
ft
ft
m
ft
ft
ft
$
»
ft
« # ft ft * ft
See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.
Downloaded via UNIV INDUSTRIAL DE SANTANDER on August 30, 2021 at 18:37:01 (UTC).
ft ft ft ft ft ft ft ft ft ft ft ft ft ft • s ft • ft ft ft ft
ft ft ft ft ft ft ft ft ft • ft ft ft ft ft # • • ft ft ft ft
ft ft ft ft • ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft # ft *
ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft * « ft ft ft ft
.
ft ft
ft ft ft ft # ft ft ft ft ft ft ft ft ft ft ft ft ft O ft ft ft ft ft ft ft ft
ft ft ft # ft ft ft • ft ft ft ft • ft * ft ft ft ft ft « ft ft ft ft # ft ft #
ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft * ft # ft ft ft ft ft ft ft ft
* ft ft ft ft ft # ft ft ft ft # ft ft ft ft ft m « « ft # ft ft ft ft ft ft -
« ft ft ft ft ft ft ft ft ft ft ft ft ft * ft ft ft ft ft ft # ft ft ft ft
The ability to quantify a trace ele- approaches (1-13). In spite of the fact examine the statistical meaning of
ment molecule in chemical and bio-
or that the International Union of Pure limit of detection values in a format
logical matrices using specific analyti- and Applied Chemistry (IUPAC) consistent with the IUPAC definition.
cal methods is often viewed in terms adopted a model for the limit of detec- It is not intended to be a comprehen-
of the limit of detection. This limit of tion calculations in 1975 (7), and the sive review of the various methods of
detection is a number, expressed in ACS Subcommittee on Environmental calculating limit of detection values.
units of concentration (or amount), Analytical Chemistry reaffirmed this For a more complete review, the read
that describes the lowest concentra- standard in 1980 (2), acceptance of er is referred to the excellent articles
tion level (or amount) of the element this model by the general analytical by Kaiser (3-6), Boumans (7-9), Cur-
that an analyst can determine to be community has been slow. The result rie (10), Glaser et al. (11), as well as
statistically different from an analyti- of this slow acceptance has led to a available textbooks (12,13). Rather,
cal blank (1). Although this definition great deal of uncertainty when limits this REPORT is intended to be a sim-
seems rather straightforward, signifi- of detection are used as a basis for ple and general discussion on methods
cant problems have been encountered comparison between various analytical for calculating limits of detection; it is
in expressing these values because of procedures, methods, or analytical in- geared to the analyst who does not
the various approaches to the term struments. Unless the limits of detec- have a rigorous knowledge of statis-
“statistically different.” The calculat- tion are calculated in a consistent tics. In this examination, the signifi-
ed limit of detection for an element manner, the comparison may be cance of the limit of detection values
can easily vary an order of magnitude meaningless. will be emphasized, and the possible
through the use of different statistical It is the purpose of this REPORT to problems encountered when using
ft#*#########*####**#*#*####
ft###*#############*####*## ft #
######## ^###**#*##»##«######»#####« ft 4
####«#«#<; • •ft**##*##***##########*## «##••*#•
####### It*############*#*****###### & * * # # « * # ft
ft#################**###### ft#######
ft ft ft * ft * ft » •##*•#•#**#*####*••#•#*•*# »########
########
####*#«# l#« ••#••••••### ##•«##•#•###
The limit of detection is the lowest concentration level that can be determined
to be statistically different from an analytical blank. Significant problems * ft *
have been encountered in expressing these values because of the various # * * *
approaches to the term “statistically different” ###••##«
######*#
3 * ft ft # * ft #**##* #
# # # * # # i
.##•#*
#
* # # « # ® # ##«#####
# # # # # # * # # * n £#######
m # * m # ft#######
# # # » # ######« p ##*#•« #
* # # * « E # # # * ft ft ft#######
# # # • # # ft ft # ft ft S
limit of detection values obtained tected with reasonable certainty for a generally plots of signal, x, vs. analyte
from non-IUPAC methods as a basis given analytical procedure” (J). This concentration, c, and are represented
for comparing methods and instru- concept is further clarified by the ACS as in Figure 1. The relationship be-
ments will be discussed. Also, two definition, which states “the limit of tween x and c can be obtained by per-
methods will be introduced that per- detection is the lowest concentration forming a linear regression analysis on
mit measurement errors in the analyt- of an analyte that an analytical pro- the data. This analytical calibration
ical sensitivity to be included in the cess cart reliably detect” (2). To un- relationship can be expressed as
limit of detection calculation. These derstand what a reasonably certain
two methods as well as the IUPAC measure or a reliable detection is, the x =
me + i (1)
method will be evaluated for their method of measurement as well as the
ability to incorporate these errors into errors (including noises) associated where m is the slope or analytical sen-
the numerical limit of detection. with the measurement must be well sitivity and i is the intercept. When an
understood. unknown sample containing the ana-
Definition lyte is subjected to the analytical pro-
cedure, a value, xy, can be measured.
The IUPAC definition, adopted in Measurements This value may then be inserted into
1975, states that “the limit of detec- Generally, most analytical methods Equation 1 to determine the concen-
tion, expressed as a concentration cl require the construction of analytical tration of the unknown, cu- However,
(or amount, qy), is derived from the calibration curves for the determina- the ability to solve accurately for cu is
smallest measure, xy, that can be de- tion of unknowns. These curves are dependent upon how well the m and i
calibration curve or a poor choice of level of 99.86% that xl & (xb + 3sb)
calibration curve ranges, the result of for a measurement based on the error
the unknown determination, cu, may of the blank signal following a normal
be subject to considerable error. distribution. It must be emphasized
that if xb does not follow a normal
Statistics distribution, then the probability
Figure 3. Analytical calibration curve that xl > (xg + 3sg) would be
The amount of error associated with of signal, x, vs. concentration, c, show- 100(1 1/k2), or 89% according to
—
a measurement of x can be statistical- ing the relationship of ksB to the limit of Tschebyscheff’s inequality (5). Hence,
ly estimated. Most measurements are detection, cL values of k < 3 should not be used for
subject to error that follows a normal limit of detection calculations.
distribution. If a sufficiently large
number of observations is made, plot- Other Approaches
ting the measured responses would The majority of the other approach-
nB
produce a curve similar to that shown es to calculating cl values are similar
in Figure 2. The mean value of the re- E (xgj -
that a measured x value would fall normal distribution, a plot of these re- values (3, 4). This value has also been
somewhere under the curve). sponses (frequency of occurrence vs.
agreed upon by other authors (9, 13),
The relationship between area and xb values) would resemble Figure 2. by IUPAC (1, 14), and by the ACS (2).
probability can be measured to esti- The probability that the smallest dis- A value of 2 for k had been initially
mate the chance that a newly mea- cernible analytical signal, xl, can be suggested (8) but this value corre-
sured x value, xg, would be a certain measured and not be a random fluctu- sponds to a 97.7% confidence level for
number of standard deviation units ation of the blank is dependent upon normal distribution and 75% for a
away from the mean response, yx. In how many standard deviation units xl nonnormal distribution of measure-
Figure 2, xE is shown to the right of yu is from xg. If xl is 3sb away from xB, ment error.
and can be measured to be k<r away the area to the right of xl is no less Although the use of k 3 instead of
=
from n. By dissecting the curve with a than 0.0013. Thus, there is a 0.13% k 2 slightly increases the cl value, it
=
line drawn at this x value, the area to chance that a signal measured at xg or is clear that cl values must differ by a
the right of the line, a, is the probabil- greater would be the result of a ran- factor of three for the values to be sig-
ity that xE > (m + kcr). This chance, dom fluctuation of the blank signal. nificantly different. Nevertheless, fac-
which is represented by the pink shad- This small chance of error can then tors of less than three have been com-
ed area, can be determined from fulfill the requirement of a reasonably monly used for comparison purposes.
certain signal. In order to minimize confusion,
area =
(2) In defining cl, IUPAC states that IUPAC suggested that xl values be re-
ibX~exp(::r)dk xl =
xB + ksB (5) ported in all literature with their k
where the xE value is kc away from n, value, xg(k—si* It would be extremely
(i.e., (xE n)/a k).
-
= where k is a numerical factor chosen in useful to go one step further and in-
This illustration can also be used to accordance with the confidence level clude the k values when cl values are
aid in the explanation of the smallest desired. The cl is a function of xl and reported, CL(k=3)- This change would
detectable signal, xl, in the IUPAC therefore be beneficial because cl values are
definition. When the determination of more commonly reported than xL
(xL -
xb)
a limit of detection is performed, cl =-
_
(6) values.
m
blank measurements, xg, are normally A problem encountered in the com-
taken. The question, however, is “how where m is the analytical sensitivity. parative use of cl values is the use of
well are these xg values known?’’ A Because the mean blank reading, xB, the standard deviation of the mean, sB
mean value of the blank responses, xg, is not always 0, the signal must be (11), the pooled standard deviation, sp
can be calculated as background corrected. By substituting (7, 12), or the relative standard devia-
"B
Equation 5 into Equation 6, Equation tion (RSD) (7, 9). Although each of
7 is obtained these standard deviation expressions
E xBj is important and has its place in ana-
xg = J—-- (3)
nB
(7) lytical chemistry, the use or misuse of
these expressions in cl calculations
and the standard deviation as This definition of cl can be illustrated may result in significant deviation
SB =
SB
(8) determination,” cb (7). Here, the mea- ysis procedure would be meaningless.
(nB)1/2 sure xd corresponding to cd is defined A well-based but seldom used con-
as cept in the calculation of detection
This value is calculated by dividing limits is the limit of identification, cB
the standard deviation of the blank ko sb xd =
(12) as described by Boumans (7) or the
measures, sB, by the square root of the where sB is the standard deviation of limit of guarantee for purity, eg, de-
number of blank observations, nB. the blank (not the RSD), and kp, is a scribed by Kaiser (3). These concepts
Usually when sp values are used, the k confidence factor that is defined by are essentially the same and are based
factor is replaced by a t distribution the reciprocal of the magnitude of the on the idea that the lowest statistical-
value. Even though the t values reflect
RSD value. For a maximum allowed ly discernible signal should be
the problem of defining a standard de- RSD of 5%, ko 20 and for a RSD of
=
viation for a finite number of observa- 10. xi =
xL + 3si (15)
10%, kD =
tions, 3sb > ta SB/(nB)1/2 for a = 0.005 Another limiting expression involv-
and five observations (four degrees of where si is the error associated with the
freedom). This inequality greatly in- ing RSD values involves the rearrang- measurements at the xi level, or
creases as nB increases. If a large num-
ing of Equation 7 (9). If the back-
ber of observations are taken on the ground signal is included in the nu- xi =
xl + 3sb (16)
merator and denominator of Equation
blank, say 30, the use of the t and sb 7,then if si =
sB, This concept is illustrated in
will reduce the value of cl by a factor
of six from the IUPAC model.
The use of the pooled standard de-
viation includes the number of blank
measurements as well as sample mea-
surements. Generally, the pooled stan- Table I. Statistical Expressions for the Slope and Intercept
dard deviation, sp, is expressed as
« n n 2
Sp -
SB (10)
ns nB \ l E Xj
/
For ng 1, sp is essentially the same
= Sxx ~
E x2
!j=i rr
as 3b- However, if ns and nB are both Sm
(Sec)1
large, Sp < sg. Again, the use of a large
number of observations as mentioned
above would cause a significant reduc-
tion in the cl value as compared to the
IUPAC model.
(Mfv
(ns^)1'2
Another common practice involves
the use of the RSD, which is defined
as
A B TARF
(20)
(21)
CIRCLE 33 ON READER SERVICE CARD
m —
350 Si
=
545 m 614=
Si
=
1200 value; this deviation is the result of
Sm *“10 sB = 3.5 sm = 24 sB = 97.6
the error associated with the intercept
i = -628 i = -980
value.
The values obtained in case B em-
phasize the problem of an ill-defined,
nonzero intercept. This problem is the
direct result of constructing calibra-
tion curves for detection limits when
the lowest point of the calibration
Table III. curve data is considerably removed
cMk=3) Values (ppm) from the CL(k ?,) value. Only the
=
Equation 7, which is the IUPAC defi- 0.03 ppm. Repeating the cp calcula-
cL =
ksc (22) nition of the limit of detection. tion using the graphical approach re-
sults in a value of 0.3 ppm. Finally al-
where the sB/m term describes the Evaluation of Approaches
lowing m and i errors to be included
error in terms of c (if c Cp). Equa- ~
The IUPAC, graphical, and propa- results in a cp(k » of 5 ppm for the
tion 21 could now be used to evaluate gation of errors models will be applied propagation of errors method. The
sc where x is the blank signal, xb, and to four different sets of experimental propagation of errors Cp value is 170
sB is substituted for sx. By measuring data to show the effect of certain ex- times the IUPAC cp<k m value.
»
xB, Sb, and calculating m, i, sm, and s„ perimental conditions on the estima- The problem of nonlinearity can be
the value of sc can be determined. tion of cp values (Table II). The data further emphasized by considering
In most determinations, the data in Table 11 have been taken from a re- case D. Here, the errors associated
are background corrected, that is, cent paper on ICP-excited ICP fluo- with m and i are greater; however, the
0- Substituting the above mea- detection limits (15). The
xb large value of sm results in a special
=
rescence
sured and calculated values in Equa- four sets of typical experimental con- problem with the graphical approach.
tions 21 and 22, the expression for cp ditions are: A, Ca(II) fluorescence If t„sm > m, then the concentration
is further simplified to data which have well-defined m and i value for the limit of detection can
1/2 values; B, Ca(II) fluorescence data even be negative. Such negative values
k SB* + Sp +
—-
Sm* where sB is essentially the same as in are the direct result of the graphical
(mj A, but the calibration curve data are model not being statistically valid. Al-
taken far away from eg resulting in a though the graphical model is easier to
system of compatible
instrumentation.
The catalogue TL-10 »lnstrumental
Thin-Layer Chromatography«
contains detailed information
on CAMAG's complete range of
TLC instrumentation.
Methodological explanations
together with an easy-to-use
guide of the instruments and
their compatibility will help you
to compose the COMPLETE
SYSTEM that suits your James D. Winefordner (/.) is a gradu- ment of new instrumental methods
requirements best. ate research professor in chemistry at based upon new optical and detection
the University of Florida. His re- approaches; and applications of ana-
Ask for this 56-page cataloque, interests include diagnostical mea- lytical techniques.
surements of flames and plasmas;
or even better, laser-excited luminescence, pho- Gary L. Long is a postdoctoral asso-
have the CAMAG product toacoustic and photoionization of ciate in the department of chemistry
specialist (of the CAMAG molecules in gases and liquids; laser- at the University of Florida. He re-
distributor in your territory) excited fluorescence of atoms in ceived his BS in chemistry from Wake
discuss your requirements with flames and plasmas; development of Forest University in 1978 and his
him. trace analytical methods for atoms PhD in analytical chemistry from
and molecules based upon lumines- North Carolina State University in
CAMAG has agents in almost all cence, chemiluminescence, pho- 1982. His research interests include
countries. toacoustic, and ionization methods; ICP-excited ICP fluorescence spec-
specialized methods for molecules trometry, laser-excited flame fluores-
based upon synchronous lumines- cence spectrometry, and the investi-
_£A/\A/\& cence, room-temperature phospho-
rimetry, energy transfer lumines-
gation of chemical and physical inter-
ferences in flame and plasma spectro-
Sonnenmattstr. 11 •
CH-4132 Muttenz/Swltzerland
Tel. 0 61-61 34 34 Telex 62 649
•
cence, and time resolution; develop- metric analyses.
CIRCLE 34 ON READER SERVICE CARD
724 A • ANALYTICAL CHEMISTRY, VOL. 55, NO. 7, JUNE 1983