Uncertainty in measurement - Introduction and

1 downloads 0 Views 289KB Size Report
In the examples in this document we will use a worksheet from. Excel 5.0 that ... sample with a pipette with a relative uncertainty of 5 %. Thus, we now have each ...
How to Cite this artic le: Unc ertainty in Measurement, Introduc tion and Examples, Kallner A, eJIFCC vol 1 3 no1 : http://www.ifc c .org/ ej ifc c /vo l1 3 no 1 /1 3 0 1 2 0 0 1 0 3 .htm Artic le ID 1 3 0 1 2 0 0 1 0 3

Uncertainty in measurement Introduction and examples from laboratory medicine

All modern c omputers and c alc ulators have predefined routines for the c alc ulation of the mean and SD. However, the c alc ulated properties are only estimates of the true standard deviation and mean of the distribution. If these estimates are made from too small a number of observations the estimate will be non - reliable i.e. the unc ertainty of the standard deviation and mean will be

Author

large.

Anders Kallner, Dept. of Clinical Chemistry, Karolinska H ospital, S - 171 76 Stockholm, Sweden

When c harac terizing a result we do not only want to desc ribe how prec ise we have managed to perform the measurement ( repeatable or reproduc ible) but also how c orrec t it is i.e. the deviation of the result from the “true value”. We never know the “true

T he concept uncertainty in measurement

value” and it is therefore not possible to desc ribe the trueness of the result. This dilemma is resolved by assigning a value and c all it the “true value”. This is normally c alled “c onventional true

It is unavoidable that all dec isions, all ac tions and therefore all measurements harbour an inherent unc ertainty. For the lay man/woman the term unc ertainty is used in many c onnec tions and is a frequently used c onc ept with no sc ientific meaning.

value” ( c onvention is here used in the meaning of agreement) or ‘assigned value’. The differenc e between the mean of many measurements of the same quantity and the c onventional true value is the “ b ia s “. The statistic “ tru e n e ss “ will therefore describe the systematic error.

However, metrologists have defined unc ertainty and given it a sc ientific c ontent that is useful for all measurements. The definitio n o f unce r ta inty in me a sur e me nt ac c ording to ISO is:

The prec ision of a result ( standard error of the mean) c an be improved by inc reasing the number of observations whereas the trueness of the result c annot be influenc ed this way. In our role as analysts we are interested to find not only the random error but

‘parameter, associated with a result of a measurement that characterises the dispersion of the values that could be reasonably attributed to the measurand’.

also the systematic error. However, as a c onsumer of results of measurements we would prefer a c onc ept that would inc lude both the random and the systematic error for a single measurement i.e. both the prec ision and the trueness. This c onc ept is c alled a c c u ra c y and in laboratory medic ine often equivalent to

This definition may be diffic ult to apply to prac tic al work and we

total error . The c o n c e pt o f u n c e rta i n ty as defined above will

shall therefore expand on the definition and explain how it c an

offer a simple and tangible alternative to total error and will also

be used for desc ribing how well a measurement proc edure

give us a possibility to avoid the use of “error” to desc ribe the

performs’.

variation inherent in measurements.

The result of a measurement shall always give information about

W hy uncertainty?

an interval within whic h the results c an be expec ted with a given probability. This c an be reported as an interval, as a standard deviation ( SD) , c oeffic ient of variation or standard error of the mean ( SEM) . Most c ommonly the distribution of results is given as one SD. The SD is a desc ription of the prec ision and c an always be c alc ulated from 2 or more results but the interpretation of this value requires knowledge of the properties of the distribution of values. The most probable value within a Gaussian distribution is represented by the mean. Together with the mean, the SD fully defines the shape of a Gaussian distribution. In these distributions the mean defines the position of the distribution on the X axis and the SD its width and height. The SD is therefore a measure of the random error ( the width of the distribution) . In a Gaussian distribution about 2 /3 of all results will be found

down into smaller entities and study eac h of them. This approac h c an be used also to estimate the unc ertainty in measurements and is usually not too c omplic ated. However, as in all sc ienc es the synthesis of results that should be performed is the diffic ult part when sc ientists may make mistakes and end up with grossly erroneous c onc lusions. The c onc ept of unc ertainty as defined by ISO requires that all steps involved in a measurement are defined and evaluated with regard to their unc ertainty. This is a great help in c harac terizing and optimising the performanc e of a measurement proc edure bec ause the sourc es of unc ertainties will be revealed. The c hemist c an review the result and selec t the steps that require improvement in an orderly and reproduc ible

between -1 SD and + 1 SD c ounted from the mean of the distribution. Further, about 9 5 % will be found within the interval -2 SD to + 2 SD and about 9 9 % between -3 SD and + 3 SD.

An ac c epted sc ientific approac h to solve a problem is to break it

manner. All steps should be inc luded, preanalytic al, e.g. sampling and sample preparation, as well as analytic al, e.g. c alibrations,

Page 16 eJIFCC2001Vol13No1pp016-024

dilutions and postanalytic al e.g. transformations and c orrec tions.

rec ording is made one must c onsider the position of the eye in

The ac c umulated data shall then be c ombined. The following text

relation to the sc ale, i.e. one sourc e of unc ertainty that c an be

will desc ribe how these steps c an be systematized and performed

referred to the operator. When a c alibration func tion is estab-

in a pragmatic manner.

lished one must estimate the unc ertainty of the value assigned to the c alibrator. Possible sourc es that might influenc e the c alibration func tion are for instanc e the assigned value of the c alibrator,

Formal ac c reditation of laboratories and measurement proc edures or methods ac c ording to the ISO standards 1 5 1 8 9 ( 1 ) and 1 7 0 2 5 ( 2 ) requires that the unc ertainty in measurements is estimated. The preferred method for estimation of unc ertainties is desc ribed in ‘Guide to the Expression of unc ertainties n measure-

the fitting of the c urve, the dilution of the c alibrators etc . When an extrac tion is inc luded one must estimate the yield. Even if all sourc es of bias were eliminated an unc ertainty of the suc c ess of that proc edure remains and should be estimated.

ments’ ( GUM) ( 3 ) . All these unc ertainty sourc es must be evaluated and brought into the unc ertainty budget. An experienc ed c hemist may direc tly

Besides estimating the unc ertainty in measurements to identify

disregard some of the unc ertainty sourc es bec ause the experienc e

areas in whic h improvements should be foc used, there is a point

is that they are very small or bec ause those remaining are muc h

in estimating the unc ertainty in all measurement that produc e results for the diagnosis and management of diseases. The reason

larger.

is that the unc ertainty in measurements will inc lude a c ontribution from c alibration etc that affec t the bias as well as estimates of

The proc edure of measurement shall then be desc ribed in

the imprec ision and any pre - and postanalytic al unc ertainties.

mathematic al terms to illustrate how the different variables affec t eac h other and how their unc ertainties c an be joined in a c ombined unc ertainty. This proc edure is tric ky and requires

U ncertainty of measurement

experienc e from establishing stoc hiometric c alc ulations i.e. the qualitative and quantitative desc ription of a c hemic al reac tion

The c onc ept of unc ertainty is not a statistic al c onc ept in tradi-

and c alc ulations of c onc entrations and amounts involved.

tional sense i.e. it needs not be assoc iated with a known distribution of the data. The unc ertainty shall rather be understood as an interval within whic h the result c an be found with a given

The estimation of the size of the sourc es of unc ertainty c an be done in two different ways c alled type A and type B . The value and

probability. Thus, the result will be within the interval but all

results of both approac hes are treated similarly but the type A is

values within the interval have the same probability to represent the result. This is quite different from other distributions, e.g. the Gaussian distribution where the mean is the most probable value. Therefore, there is not nec essarily any symmetry in the c onc ept of unc ertainty. There is no value that is more probable or c ommon than any other. It is c onvenient, though, to let a number represent the distribution and one c an c hoose any representative value from

rec ommended whenever possible. In both c ases the ‘standard unc ertainty’ is estimated. The standard unc ertainty is abbreviated

u( x i ) .

Estimation of the uncertainty in measurement, type A

the interval. One is as good as the other but the c entral value has advantages from a presentation point of view. B y these statements

In this model the estimation is based on the standard deviations

we have in fac t desc ribed the essential features of a ‘rec tangular

derived from repeated measurements. The SD and the standard

distribution’

unc ertainty will therefore have the same size. A word of warning: sometimes the variation is given as a c oeffic ient of variation ( CV)

The c alc ulation or estimation of the unc ertainty in measurements

or c onfidenc e interval rather than SD. That information c an also

is princ ipally very simple. This is done by estimating the unc er-

be used but it must then be c orrec ted for what the numbers

tainty in eac h single step of a proc edure ( input variables, X i ) and

represent e.g. multiply with the result of measurement and

c o mbining them in an ‘unc ertainty budget’. The c o mbinatio n o f

divided with 1 0 0 if the CV is given in perc ent or multiplied with

the individual unc ertainties follows similar rules as the propaga-

the square root of the number of observations if the c onfidenc e

tion of errors in normal statistic s giving the ‘c ombined unc er-

interval or SEM ( standard error of the mean) is given. One must

tainty’ for the proc edure. The c ombined unc ertainty is the

also c larify how many SD the information refers to. As a rule the

c onc ept that is c losest to ‘the total error’.

number is one but sometimes multiples are given, for instanc e two SD.

In the examples in this doc ument we will use a worksheet from Exc el 5 .0 that simplifies these c alc ulations and do not require any knowledge of the mathematic s behind. The worksheet c an be used for evaluation of all unc ertainty budgets that deal with

Estimation of the uncertainty in measurement, type B

independent observations. Sometimes, and partic ularly in c omplex measurement proc edures, we do not have ac c ess to or we are not able to estimate the

Sources of uncertainties

variation from repeated experiments, whic h is a prerequisite for the estimation of the SD. The professional experienc e, informa-

It is essential when establishing a realistic unc ertainty budget to

tion in the literature or spec ific ations from a manufac turer will

identify the variables that give rise to the unc ertainty and their

usually allow the demarc ation of an interval within whic h a result

sizes. Therefore one must have detailed knowledge about the

c an reasonably be expec ted. For instanc e, it is fairly safe to

proc edure of measurement to allow identific ation and quantific a-

assume that the European woman is between 1 3 0 c m and 2 1 0

tion of all reasonable sourc es. Examples of sourc es of unc ertainty

c m tall, that a litre of milk c osts between USD 0 ,1 and 2 and that

in our field are the measurement of volumes, weighing, reac tion

one litre of water has a mass between 9 9 4 g and 1 0 0 4 g. The

temperature, purity of reagents, and value assigned to the

more one knows about the proc edure of measurement and the

c alibrator. Also properties of the instrument used e.g. the size and

items measured the better will the estimate of interval be. In

toleranc e of the c uvette, if the instrument uses more than one

laboratory medic ine it is often important to find a reasonable

c uvette they may be different, the wavelength etc . If a visual

interval for the volume, the mass or the reading of an instrument.

Page 17 eJIFCC2001Vol13No1pp016-024

For instanc e a 1 mL pipette might deliver between 0 .9 0 mL and

now have eac h volume defined to its size and ac c ompanying

1 .0 5 mL, the body mass of a grown up man without exc essive fat

unc ertainty. Observe that the unc ertainty of the measuring c ylinder

might be between 6 5 and 8 5 kg; based on your experienc e your

is given as an interval whereas the unc ertainty for the pipette is

might even assume that it is between 6 9 kg and 7 6 kg etc .

given as a relative standard unc ertainty.

The interval that is defined in this way will not identify any result

The volumes V 1 and V 2 should be added and thus the ‘reac tion’

within the interval as being more reasonable than any other in the

c ould be written

same interval. Thus, the interval represents a rec tangular distribution of possible results.

Throwing a single dic e will give results belonging to a rec tangular distribution. One c annot foresee the number of dots that will c ome up but a reasonable ( the largest possible) interval is 1 - 6 .

The standard unc ertainties are added but not as suc h but as their

In fac t no other results c an be obtained. The probability for eac h

squares:

result within the interval is the same, provided the dic e is fair. The mean of the interval ( 3 .5 ) , however, c an never be obtained!

Given an interval there is a somewhat smaller interval within whic h the probability that a result will oc c ur is inc reased. This new interval is c reated from the half - width of the initial interval divided by the square root of 3 ( about 1 .7 3 ) in both direc tions from the middle of the initial interval. It is about twic e as c ommon that results will be within this interval as outside these limits. Using the example with the dic e, the inner interval will be

The c ombined unc ertainty is then obtained by taking the square root from the sum, in this example 8 .7 mL. Perform the c alc ulations yourself and find that the c ontribution from the pipette is of minor importanc e in the c ombined unc ertainty!

Coverage factor

from 2 to and inc luding 5 . For the dic e it means that between 2 and 5 there are four alternatives ( 2 , 3 , 4 and 5 ) , whereas outside

The standard unc ertainty will thus demarc ate an interval where we

there are only two, 1 and 6 . It is therefore more probable that any

c an estimate that 2 /3 of all the results will be found. This is also

of the four dots 2 to 5 will show up than any of the two outside

the usual way to give the variation ( usually the standard deviation)

the smaller interval. The probability that any given number of dots

of results in sc ientific literature. If we want to find the interval

within the new interval will be obtained is the same as outside

that is large enough to give a 9 5 % probability to c over the results

and one c annot predic t whic h. I.e. it is more probable than one

in a Gaussian distribution the SD should be multiplied with about

of four given alternatives is obtained than one of two.

2 for a two - tailed distribution ( the size of the fac tor is depending on the number of observations) . To reac h a c orresponding

The width of this inner interval is twic e the standard unc ertainty 2

probability for the c ombined unc ertainty it shall be multiplied by

u ( x i ) . The standard unc ertainty will c orrespond approximately

2 . If we wish a larger probability then a larger fac tor, for instanc e

to the probability within the interval mean ± 1 SD in a Gaussian

3 for about 9 9 % probability, should be used. This fac tor is c alled

distribution where about 2 /3 of all observations will be found.

the co ve ra ge fa cto r ( k ) and the result e xpa n de d u n c e rta i n ty

You c an c onvinc e yourself about this by playing with a good dic e

U y . The c overage fac tor shall always be given in the answer

and you would expec t 2 /3 of all answers to be found between and

together with the unc ertainty. In c ase no c overage fac tor is given

inc lude 2 and 5 dots.

then the c ombined standard unc ertainty c overing about 2 /3 of the result is given ( k = 1 ) .

An example from our own profession is the estimation of the unc ertainty of a measured volume using a two - litre measure-

U ncertainty budget

ment c ylinder. Suppose we want to measure 5 0 0 mL, and assume a reasonable interval to be ± 3 % or ( 4 8 5 -5 1 5 ) mL. The

When estimating the c ombined unc ertainty ( u

standard unc ertainty is then 1 5 ( half the interval) divided by the

point is to define how the different parts of the proc edure

square root of 3 i.e. 8 .7 mL.

interac t. In our model we assume that the unc ertainty sourc es are

c

( y) ) the starting

independent. In the example above the two volumes were added

Combined standard uncertainty (uc-(y))

to reac h the total volume. For additions ( subtrac tions) , the c ombined unc ertainty is the square root of the sum of the squares of the ingoing standard unc ertainties. In c ase the variables shall

Onc e the sourc es of unc ertainty have been identified and the sizes

be multiplied ( divided) the squares of the ingoing relative

of the unc ertainties estimated then we fac e the problem to

standard unc ertainties shall be added. When the square root is

c ombine the c ontributions to an unc ertainty for the entire

drawn the relative c ombined unc ertainty is ac hieved.

proc edure i.e. the c ombined unc ertainty in the result y. This is abbreviated u c ( y) .

The mathematic al expressions that desc ribe how the standard unc ertainties shall be c ombined bec ome rather c omplic ated

Let us begin with an example:

when the unc ertainty budget for a proc edure c ontains all four operators and a variable may partic ipate as a logarithm or

We want to make a dilution of a sample and are only interested in

exponent. The generic rule for the c ombination is based on

the final volume. Ac c ording to the proc edure we shall take 1 0 mL

partial derivatives but we will use an approximate numeric al

of the sample solution and dilute that to 5 0 0 mL. We use the

method for solving partial derivatives. This c an c onveniently be

measuring c ylinder from the example above and measure the

ac hieved by using a spreadsheet program like Exc el® .

sample with a pipette with a relative unc ertainty of 5 % . Thus, we

Page 18 eJIFCC2001Vol13No1pp016-024

Kragten ( 4 ) originally desc ribed the numeric al approximation of

Examples

partial derivatives. The doc ument from Eurac hem ( 5 ) desc ribes how this method c an be used in a worksheet.

The following examples are c hosen to illustrate an inc reasing c omplexity in the c alc ulation and originate in routine laboratory

Estimation of the uncertainty in measurements using an Excel® worksheet

work. The unc ertainties and other numbers do not nec essarily represent reality but have been c hosen to illustrate their influenc es. Eac h example is solved on the attac hed diskette where also you will also find a template for the c alc ulations.

General

1)

All c ells in the worksheet that will be filled by c alc ulations are

Let us repeat the example above dealing with adding two volumes

protec ted to avoid unintentional c hanges in the program. Cells that c an be c hanged are marked with a blue border. The number

and use the template.

of dec imals is fixed to three that will give a suffic ient presentation of results within laboratory medic ine even if it sometimes

Enter the nominations “Sample” and “Dilution” in the c orrec t

exaggerates the prec ision. You c an perform c alc ulations direc tly

c ells in row 3 and the volumes 1 0 and 4 9 0 respec tively in row 4 .

in that c ell where the value finally will be plac ed. Never c opy the

The standard unc ertainties are entered in row 5 and 6 ( absolute

c ontents of a c ell into another, delete and input it again or

and relative standard unc ertainty, respec tively) . Move the c ursor

unforeseeable errors may oc c ur.

to c ell C2 1 under the label “Nominal” and enter the formula ‘= C9 + C1 0 ’. Copy this formula two c olumns to the right. The c ombined unc ertainty will be shown in c ell C2 5 , the relative in

It is suggested that you enter different data into the given

E2 5 and the interval c alc ulated with a c overage fac tor 2 in the

examples to see how the outc ome c hanges.

c ells K2 5 - L2 5 . The expanded unc ertainty is given in c ell I2 5 . c ompare the result with our manually worked example.

D ata entry and calculations

2) Enter the names of the input variables in row 3 of the worksheet. The input variables c an be entered in any order. The c orresponding names will appear in c olumn B , rows 1 0 to 1 8 . Enter ac tual

We want to weigh NaCl on an old fashioned sc ale. The empty vessel weighs 1 2 g. We add NaCl until 1 2 7 g. Weighing within the

or representative values in row 4 and the standard unc ertainty of the results in row 5 or 6 depending on if the standard unc ertainty is given in absolute ( row 5 ) or relative terms ( row 6 ) . Relative

given interval c an be made with an unc ertainty of 1 ,5 % . Calc ulate the mass of the NaCl.

unc ertainties shall be given as parts of 1 , i.e. 1 .5 % should be In c ell C2 1 the formula will be:’ = final weight-weight of empty

0 .0 1 5 .

vessel ( tare) ’. Test the result if the unc ertainty is 0 .2 g! Then the interrelation between the various c omponents of the budget shall be entered in c ell C2 1 ( ‘Nominal’) , i.e. how the final

3)

result shall be c alc ulated from the input variables. Often this formula presents itself easily but it c an sometimes be rather

We shall c alc ulate the volume of a water bath whic h has a length

c omplic ated and require deep thoughts. It may be helpful to

of 3 5 c m, width 2 5 c m and height 2 0 c m. The edges are

establish and solve an equation that desc ribes how to c alc ulate

measured using a ruler with an unc ertainty of 2 % .

the result. In the dilution example above the expression will be simply the volume of the original sample plus the dilution

The formula in C2 5 will be: ‘= the length x width x height’. Chec k

volume. Compare example 1 .

the result if the unc ertainty of measurement in the interval is 0 .5 c m!

The expression whic h goes into c ell C2 1 will be a mathematic al expression and shall therefore, with the Exc el® nomenc lature,

4)

begin with “= ” ( “equal sign”) . Then, enter the algorithm and drag the c ontents of C2 1 as far to the right as variables have been entered. Finally press “enter” and all c alc ulations will be c arried

Let us add a slight c omplic ation to example 3 . Let us assume that

out. Note that the c ells to the right of the last used c olumn in row

the water bath is half - filled ( 1 0 c m) and you want to fill it up to

2 1 shall be empty. If not, it is safe to mark them and press delete.

3 c m ( the margin) from the upper rim. How muc h water should be added? We make the same assumptions regarding the unc ertainty of the lengths as above, but sinc e the vessel walls are

Depending on how your c omputer has been c onfigured you may

not quite even we must add the unc ertainty this c ontributes. Let

need to press the F9 key to trigger the c alc ulations. If nec essary

us assume that this unc ertainty is multiplic ative and optimally 1

you c an c hange the c onfiguration to ‘automatic ’ under ‘Tools-

but with the unc ertainty 3 % .

Options-Calc ulation’ and c hec k ‘Automatic ’.

The formula in C2 1 will be ‘“length x width x ( the height -

The c ombined unc ertainty will be given on the last row in both absolute and relative terms. Also the c overage fac tor ( default = 2 ) and the unc ertainty interval will be given. Contributions from the

marginal - the height at half - filled water bath) x unevenness fac tor’. Test the results if the unc ertainty in all lengths is 0 .2 c m!

different sourc es of unc ertainty are shown on row 2 3 and graphic ally in the inserted diagram. The diagram c an be freely

5)

moved within the surfac e if you want to study the underlying c ontents of the table. The sc ale of the Y - axis c an also be c hanged to improve the presentation.

Page 19 eJIFCC2001Vol13No1pp016-024

Let us inc rease the c omplexity of example 1 above. Assume that

spec ific absorbanc e at 5 6 0 nm will therefore be ( 0 .1 2 5 ± u a -2 0

the sample c onc entration is 3 2 0 mmol/L with a standard

x 0 .0 0 1 ± u l ) . We c an then formulate the absorbanc e at 5 8 0 nm

unc ertainty of 3 % . What is the c onc entration in the final solution

± 5 nm. It will be ( 0 .1 2 5 ± u a -2 0 x 0 .0 0 1 ± u l ) + 5 8 0 nm ±

and what is its c ombined unc ertainty? Enter the numbers to the

5 nm-5 6 0 ) x ( 0 .0 0 1 ± u l ) .

template as in example 1 but add the c onc entration of the sample solution in a c olumn of its own.

Although the wavelength has a small unc ertainty this will be the dominating sourc e of unc ertainty. Disc uss what happens if the

The formula will be: ‘= sample x c onc entration/( sample +

slope inc reases or if the unc ertainty in the volume measurements

dilutio n) ’.

is c hanged.

6)

10)

We want to estimate the unc ertainty in an HPLC - method. The

The c onc ept of unc ertainty c an also be used for other purposes

sample is diluted 1 + 9 with a so lutio n c o ntaining 0 .2 5 mmo l/L

than measurements and it c an inc lude for instanc e preanalytic al

internal standard ( IS) . We injec t 2 0 µL of the dilution. The

and post analytic al sourc es of unc ertainty. Let us examine this

sample peak is 2 4 7 mm, the IS - peak is 2 3 5 mm. How muc h

example that is a simplific ation of an artic le in “Clinic al Chemis-

substanc e did the sample c ontain and whic h is the c ombined

try and Laboratory Medic ine” ( 6 ) .

unc ertainty of the result? The standard unc ertainty in the sample vo lume is ± 3 % and fo r the dilutio n vo lume ± 2 % . The standard unc ertainty in weighing of IS is ± 1 % . Measuring of the peaks is assoc iated with an interval of ± 2 % . The same amounts of IS and sample give the same response on the printer with an

B - Gluc ose is used in the primary health c are for diagnosis and c ontrol of diabetes. Whic h is the smallest differenc e in results that, with a given probability c an be assumed to be different?

unc ertainty of ± 3 % . The patient shall be fasting. The patient shall be c alm and relaxed to avoid that c atec holamines and other hormones will give a

7)

falsely elevated gluc ose c onc entration. The c onc entration of gluc ose in erythroc ytes is less than in plasma and therefore one

The relation between signal and c onc entration ( i.e. the c alibra-

also has to c ontrol the intake of fluids to avoid falsely low results

tion func tion) is: signal = c onc entration x 2 -1 ,5 ; ( Y= 2 x X-1 ,5 )

and inc reased hematoc rit. The c apillary sampling is diffic ult. One

in the c onc entration interval 0 ,5 -5 units. The standard unc er-

c annot avoid a varying addition of interstitial fluid etc . There are

tainty in the slope ( b) is 2 % in the interc ept ( a) 5 % and in the

instruments on the market that measure P - Gluc ose and use an

signal 0 ,0 1 units. Estimate the c ombined unc ertainty in the

algorithm to transfer the value to B - Gluc ose. The algorithm is

middle of the interval, i.e. 2 ,7 5 units.

based on the assumption that P - Gluc ose is 1 5 % higher than B Gluc ose, whic h however is depending on the hematoc rit of the patient. Finally we also have unc ertainty c ontributions from the

8)

measurement itself. All these unc ertainties shall be c onsidered in

We c ompare the results from measuring the same samples on two different instruments. The results are assumed to be identic al,

estimating the c ombined unc ertainty that shall be the basis for estimating the smallest signific ant differenc e between two results.

that is the regression func tion is assumed to be Y = X. Whic h deviation in Y results c an we expec t if the standard unc ertainty in

Let us selec t 6 .0 mmol/L as an interesting c onc entration.

the estimation of the c oeffic ient is 5 % and the interval of the interc ept is -0 ,5 - + 0 ,8 at the c ritic al limit 1 7 mmo l/L? The

Fa sti n g a n d f lu i d i n ta k e . Just standing up c an c ause to up to

standard unc ertainty in the referenc e method is 3 % .

1 5 % c hanges of the plasma volume. B etween day variation of the plasma volume is given in the literature ( 7 ) to ± 6 .5 % based on repeated measurements. Let us assume that this inc ludes the

9)

variations in stress, fluid intake and fasting but exc luding whether the patient is sitting, standing or resting.

The absorbanc e of a sample is measured at 5 8 0 nm. The spec ific absorbanc e for the substanc e at this wavelength is 1 2 ,5 , the absorbanc e c urve is linear between 5 4 0 and 6 2 0 nm with a slope of 0 .0 0 1 absorbanc e units/nm. The sample was diluted 1 + 9 and the absorbanc e 1 .3 5 . Estimate the c onc entration and its

Sa m pli n g: Literature ( 8 ) postulates an unc ertainty of ± 3 .2 % for S - Gluc ose also based on repeated measurements. For c apillary sampling it is reasonable to add a little, ± 5 % .

unc ertainty if the standard unc ertainties in the volume measurements are 2 % and 1 ,5 % per sample and dilution, respec tively,

Me a su re m e n t: The prec ision of the measurement is related to the

the standard unc ertainty of the spec ific absorbanc e ( u a ) 0 ,5 % ,

instrument. An instrument like Hem c ue c an give an imprec ision

in the reading 1 % , in the slope ( u l ) 5 % and the wavelength of

of ± 3 - 5 % , many of the other instruments the same order and

the filter is given as an interval of 5 8 0 nm ± 5 nm.

magnitude if one uses the same batc h of reagent strips. B etween batc hes the manufac turers may allow as muc h as up to 8 %

There are many fac tors involved in this example and let us argue like this: We must somehow translate the unc ertainty in wave-

imprec ision. Let us assume that we c an manage an interval of ± 5 % unc ertainty that also inc ludes the unc ertainty of c alibration.

length to unc ertainty in spec ific absorbanc e. This c an be done by c alc ulating the spec ific absorbanc e at another wavelength. Let us

Po st a n a lyti c a l If the measurement is made with an instrument

assume 5 6 0 nm. This wavelength has no unc ertainty bec ause we

that really measures B - Gluc ose this sourc e of unc ertainty is not

assume it is without unc ertainty. At exac tly 5 8 0 nm the spec ific

interesting. If, however the instrument measures something

unc ertainty is 0 .1 2 5 . We c an inc lude the unc ertainty in the slope

different then one must divide with a fac tor. Some literature

of the absorbanc e c urve by ( 5 8 0 - 5 6 0 ) x 0 .0 0 1 ± u

suggests 1 .1 7 but it is c loser to 1 .1 1 ( 6 ) . Regardless of the value

l .

The

of the fac tor it is attac hed to an unc ertainty, let us assume 5 % .

Page 20 eJIFCC2001Vol13No1pp016-024

T hree important definitions by ISO :

The smallest differenc e that c an be c alled signific ant is

where D is the differenc e, u c ( y) is the c ombined unc ertainty and k the c overage fac tor. The square root is bec ause two values are

Trueness, agreement between the average value from many observations and the true value

c ompared. A c onvenient abbreviation is D= 3 x u c ( y) .

Precision, agreement between independent results of measurement

References 1 . ISO/DIS 1 5 1 8 9 . Quality management in the medic al laboratory, ISO, Geneva 2 0 0 0 .

Accuracy, agreement between the result of a measurement and a true value of the measurand

2 . ISO/IEC IS 1 7 0 2 5 . General requirements for the c ompetenc e of testing and c alibration laboratories. Geneva 1 9 9 9 . 1 . Guide to the Expression of Unc ertainty in Measurement ( GUM) . ISO Geneva, Sc hweiz 1 9 9 3 ( ISB N 9 2 - 6 7 - 1 0 1 8 8 9)

Note the difference between Trueness and ; the latter refers to one measurement whereas the former to the mean of many measurements. Therefore inaccuracy will include both the bias and the imprecisio n inherent in a specific result. Accuracy

2 . J Kragten. Calc ulating standard deviations and c onfidenc e intervals with a universally applic able spreadsheet tec hnique. Analyst 1 9 9 4 ; 1 1 9 : 2 1 6 1 - 2 1 6 6 . 3 . Eurac hem, Quantifying Unc ertainty in Analytic al Measurement, LGC Information Servic e, Teddington Middlesex 1 9 9 5 ISB N 0 - 948926 - 08 - 2 4 . A. Kallner, J Waldenstrom Does the Unc ertainty of Commonly Performed Gluc ose Measurements Allow Identific ation of Individuals at High Risk for Diabetes? Clin Chem Lab Med 1 9 9 9 ;3 7 :9 0 7 - 9 1 2

Terminology and nomenclature are cornerstones of science and a citation from L Carroll ‘Alice through the looking glass’ might be appropriate to consider:

5 . C Ric os et al. Current databases on biologic al variation: pros, c ons and progress. Sc and J Clin Lab Invest 1 9 9 9 ;5 9 ( 7 ) :4 9 1 500. 6 . X. Fuentes - Arderiu et al. Pre - metrologic al ( Pre - Analytic al)

‘When I use a word’ Humpty Dumpty said in rather a scornful tone, ‘it means just what I choose it to mean ? neither more nor less’

Variation of some B ioc hemic al Quantities. Clin Chem Lab Med 1 9 9 9 ;3 7 :9 8 7 - 9 8 9

‘The question is,’ said Alice, ‘whether you can make words mean different things’.

‘The question is, said Humpty Dumpty, ‘ which is to be master ? that’s all’.

Page 21 eJIFCC2001Vol13No1pp016-024

2VlNHUKHWVEXGJHW

&RS\ULJKW$QGHUV.DOOQHU

Independent variables:

Variable: Value: Stand.unc., constant: Stand.unc., relative:

6WRFN 

Stand.unc.



Variable: 6WRFN 'LOXWLRQ

Value:  



 

'LOXWLRQ  

1.000



0.800

Relativae contributions of uncertatinty

0.900

0.700

 

0.600

0.500

0.400

0.300

0.200

0.100

Nominal  Part:

 

 

0.000 D

E

F

G

H

Comb unc.:  Rel. unc.:  k=  Exp.unc.:  Unc.interval: &RS\FHOO&WRWKHULJKWDVIDUDVWKHUHDUHHQWULHV'RQRWFRS\FHOOVGHOHWHWKHFRQWHQWVDQGUHHQWH '

(

)

*

+

,

Page 22 eJIFCC2001Vol13No1pp016-024

-

.

I

J

K



L



/

([HPSHO

&RS\ULJKW$QGHUV.DOOQHU

Independent variables 7DUD 7RWZHLJKW    

Variable: Value: Stand.unc., constant: Stand.unc., relative:



Stand.unc.: Variable: 7DUD 7RWZHLJKW

Value:  



 

Relativae contributions of uncertatinty

 

1.000 0.900 0.800 0.700 0.600 0.500 0.400 0.300 0.200 0.100 0.000 D

Nominal  Part:

 

E

F

G

H

J

K

L

 

Comb unc.:  Rel. unc.:  k=  Exp.unc.:  Unc.interval: &RS\FHOO&WRWKHULJKWDVIDUDVWKHUHDUHHQWULHV'RQRWFRS\FHOOVGHOHWHWKHFRQWHQWVDQGUHHQWH '

I

(

)

*

+

,

Page 23 eJIFCC2001Vol13No1pp016-024

-

.





/

([HPSHO

&RS\ULJKW$QGHUV.DOOQHU

Independent variables

Variable: Value: Stand.unc., constant: Stand.unc., relative:

/HQJWK  

:LGWK  

+HLJKW  

Stand.unc.:







Variable: /HQJWK :LGWK +HLJKW

Value:   

  

  

Relativae contributions of uncertatinty

  

1.000 0.900 0.800 0.700 0.600 0.500 0.400 0.300 0.200 0.100 0.000 D

Nominal     Part:   

E

F

G

H

I

J

K

L

k=  Exp.unc.:  Unc.intervall   Comb unc.:  Rel. unc.:  &RS\FHOO&WRWKHULJKWDVIDUDVWKHUHDUHHQWULHV'RQRWFRS\FHOOVGHOHWHWKHFRQWHQWVDQGUHHQWH '

(

)

*

+

,

Page 24 eJIFCC2001Vol13No1pp016-024

-

.

/

Suggest Documents