45 Feinberg Valid
45 Feinberg Valid
45 Feinberg Valid
Max Feinberg1
Rsum : partir des donnes collectes pour construire le profil dexactitude, il est possible de
calculer divers critres complmentaires souvent exiges pour les dossiers de validation, comme la
limite de quantification ou lincertitude. Sappuyant sur les concepts dcrits dans dautres articles
consacrs la validation des mthodes : principe et mise en uvre du profil dexactitude cet article
prsente diffrents complments et, en particulier, le mode de calcul de la rptabilit et de la fidlit.
Notations et symboles
Une srie de mesurages i, (i variant de 1 I) contient des niveaux de concentration k (k variant
de 1 K) pour lesquels des rptitions ont t faites j (j variant de 1 J) les indices seront
repris dans lordre i, j, k. Les variables alatoires sont notes en majuscules et les valeurs
quelles prennent en minuscules.
Symbole Description
xijk Valeur de rfrence assigne un talon pour la srie i (1 i I), la rptition j
(1 j J) et le niveau k (1 k K)
ou
Valeur de rfrence assigne un chantillon de validation pour la srie i, la
rptition j et le niveau k.
yijk Mesure de la rponse instrumentale ou exprimentale observe pour un talon ou
un tmoin de validation pour la srie i, la rptition j et le niveau k.
zijk Valeur retrouve dans lchantillon de validation de la srie i, pour la rptition j
et le niveau k, obtenue soit par prdiction inverse laide dun modle
dtalonnage, soit par une mesure directe.
bijk Biais exprimant lerreur de justesse pour un tmoin de validation entre sa valeur
retrouve et sa valeur de rfrence bijk = zijk xijk
1 UR1204 Met@risk - Mthodologie danalyse de risque alimentaire INRA F-75231 Paris cedex 05
01 44 08 16 52 max.feinberg@paris.inra.fr
zi =1 j =1
ij
( z
i =1 j =1
Z )2 = ( zij Zi )2 + J ( Zi Z )2
ij
i =1 j =1 i =1
(2)
SCEt SCEr SCEB
Cette quation est traditionnellement crite sous une forme abrge, faisant appel trois
sommes de carrs dcarts (SCE).
SCEt = SCEB + SCEr quation gnrale de lanalyse de variance
O chacune des sommes est dfinie comme suit pour faciliter linterprtation :
SCEt Somme totale des carts la moyenne gnrale du niveau
SCEB Somme des carts inter-sries
SCEr Somme des carts intra-srie
Comme le propose la norme ISO 5725, il nest pas ncessaire de dvelopper les trois sommes
de carrs pour effectuer les calculs ; le calcul de SCEB se fait par diffrence. Cette mthode
peut poser des problmes si le rsultat est ngatif : forcer alors la valeur de SCEB 0.
SCEB = SCEt SCEr si SCEB > 0
SCEB = 0 si SCEB 0
n 2
i
N* = N i =1
I
Avec N et N = ni
i =1
Diffrents exemples de calcul sont disponibles dans la littrature (Feinberg M., 2009 -
Commission SFSTP 2003, 2006). Lexemple prsent porte sur une mthode de dosage du
nicotinamide (vitamine B3) dans le lait par HPLC avec une dtection fluorimtrique. Cest
une mthode indirecte qui ncessite un talonnage pour raliser la quantification. Les 10
tapes de la procdure se droulent ainsi.
4. Plan dtalonnage (tableau 2). Le plan dtalonnage comporte I = 3 sries (jours) pour K
= 2 niveaux de concentration et J = 2 rptitions par jour et par niveau, soit un total de 12
essais. Les deux niveaux sont cods Bas et Haut.
6. Calcul des concentrations prdites inverses et calcul de la justesse (tableau 4). partir des
coefficients du tableau 3 et des mesures des surfaces de pic des chantillons de validation
rassembles au tableau 1, calculer les concentrations prdites inverses de chaque rptition
en prenant soin dutiliser les modles obtenus pour chaque srie.
Tableau 4. Concentrations retrouves par prdiction inverse partir des donnes du plan de
validation (exprimes en mg/l)
Srie Niveau Concentration X Rponse Y Retrouve Z Biais absolu Biais relatif
Jour 1 A 0.40 22.6 0.40 -0.004 -1.06%
Jour 1 A 0.40 22.1 0.39 -0.011 -2.82%
Jour 1 A 0.40 22.4 0.39 -0.007 -1.76%
Jour 2 A 0.40 23.3 0.40 0.004 0.89%
Jour 2 A 0.40 24.1 0.42 0.015 3.75%
Jour 2 A 0.40 23.9 0.41 0.012 3.04%
Jour 3 A 0.40 23.8 0.43 0.026 6.47%
Jour 3 A 0.40 23.6 0.42 0.023 5.75%
Jour 3 A 0.40 23.5 0.42 0.022 5.39%
Jour 1 B 2.00 135.0 1.98 -0.021 -1.04%
Jour 1 B 2.00 135.1 1.98 -0.019 -0.97%
Jour 1 B 2.00 129.9 1.91 -0.093 -4.63%
Jour 2 B 2.00 137.6 2.04 0.037 1.85%
Jour 2 B 2.00 135.2 2.00 0.003 0.14%
Jour 2 B 2.00 138.8 2.05 0.054 2.71%
Jour 3 B 2.00 136.3 2.04 0.043 2.13%
Jour 3 B 2.00 135.1 2.03 0.025 1.27%
Jour 3 B 2.00 134.4 2.02 0.015 0.77%
Jour 1 C 4.00 275.2 3.95 -0.046 -1.14%
Jour 1 C 4.00 276.9 3.98 -0.022 -0.55%
Jour 1 C 4.00 261.3 3.76 -0.242 -6.04%
Jour 2 C 4.00 268.1 3.90 -0.098 -2.45%
Jour 2 C 4.00 269.7 3.92 -0.075 -1.88%
Jour 2 C 4.00 276.9 4.03 0.028 0.70%
Jour 3 C 4.00 271.6 3.99 -0.013 -0.32%
Jour 3 C 4.00 273.3 4.01 0.011 0.29%
Jour 3 C 4.00 275.0 4.04 0.036 0.90%
7. Calculer les biais individuels. On observe que les biais individuels varient entre -6.04% et
6.47% et permettent didentifier les mesures ventuellement atypiques. Aucune donne
atypique na t identifie.
8. Calcul des donnes de fidlit et des intervalles de tolrance. Calculer les critres de
fidlit, comme dcrit dans le 1et les limites des intervalles de tolrance selon la mthode
dcrite dans Feinberg M. (2010b). Rassembler les valeurs trouves dans un tableau
rcapitulatif organis comme le tableau 5.
Tableau 5 : tableau rcapitulatif des critres de validation et des lments du profil dexactitude
Probabilit tolrance () 80%
Limite d'acceptabilit () 10%
Niveaux Niveau A Niveau B Niveau C
Concentration moyenne thorique 0.4 2.0 4.0
Concentration moyenne retrouve 0.409 2.005 3.953
Fidlit
cart-type de rptabilit (sr) 0.004 0.030 0.081
cart-type inter-sries (sB) 0.015 0.039 0.033
cart-type de fidlit (sFI) 0.016 0.049 0.087
Coefficient de variation de fidlit 3.90% 2.45% 2.21%
Justesse
Biais (%) 2.18% 0.25% -1.17%
Intervalle de tolrance
Nombre de degrs libert 2.197 3.374 6.826
Facteur de couverture (kIT) 1.837 1.599 1.419
cart-type de tolrance (sIT) 0.018 0.055 0.093
Limite intervalle tolrance basse 0.375 1.917 3.821
Limite intervalle tolrance haute 0.442 2.093 4.086
Limite d'acceptabilit basse 0.3600 1.8000 3.6000
Limite d'acceptabilit haute 0.4400 2.2000 4.4000
Profil d'exactitude
Recouvrement % 102.2% 100.2% 98.8%
Limite intervalle tolrance basse en % 93.8% 95.9% 95.5%
Limite intervalle tolrance haute en % 110.6% 104.6% 102.2%
Limite d'acceptabilit basse en % 90.0% 90.0% 90.0%
Limite d'acceptabilit haute en % 110.0% 110.0% 110.0%
110%
105%
Recouvrement (%)
100%
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
95%
90%
85% Concentration
Figure 1 : profil dexactitude pour le dosage de la vitamine B3 dans le lait par HPLC
A B C D E F G
1 Niveau A
2 Concentration moyenne thorique 0.4000
3 Sries Rptition 1 Rptition 2 Rptition 3 ni SCEi
4 Jour (Srie) 1 0.396 0.389 0.393 3 0.00003
5 Jour (Srie) 2 0.404 0.415 0.412 3 0.00007
6 Jour (Srie) 3 0.426 0.423 0.422 3 0.00001
7
8 Nombre de sries (I) 3.0000 =NB(B4:B6)
9 Nombre de mesures (IJ) 9.0000 =NB(B4:D6)
10 Nombre de rptitions (J) 3.0000 =SI(F4*B8<>B9;"Calcul impossible";F4)
11 SCE rsiduelle 0.00011 =SOMME(G4:G8)
12 SCE totale 0.00156 =SOMME.CARRES.ECARTS(B4:D6)
13 SCE inter-sries 0.00145 =B12-B11
14 Calcul intermdiaire de sB 0.00024 =((B13/(B8-1))-B15)/B10
15 Variance de rptabilit (sr) 0.00002 =B11/(B9-B8)
16 Variance inter-sries (sB) 0.00024 =SI(B14<0;0;B14)
17 Variance de fidlit (sFI) 0.00025 =SOMME(B15:B16)
18 Fidlit
19 Concentration moyenne retrouve 0.40873 =MOYENNE(B4:D6)
20 Ecart-type de rptabilit (sr) 0.00419 =RACINE(B15)
21 Ecart-type inter-sries (sB) 0.01536 =RACINE(B16)
22 Ecart-type de fidlit (sFI) 0.01592 =RACINE(B17)
23 Justesse
24 Biais (%) 2.18% =B19/B2-1
25 Facteur de couverture de l'intervalle de tolrance
26 Rapport des variances (R) 13.40469 =B16/B15
27 Coefficient B 0.34951 =(B26+1)/(B10*B26+1)
28 Nombre de degrs libert 2.19709 =(B26+1)^2/((B26+1/B10)^2/(B8-1)+(1-1/B10)/B9)
29 Probabilit tolrance (bta) 80%
30 t Student bas 1.88562 =LOI.STUDENT.INVERSE(1-B29;ARRONDI.INF(B28;0))
31 t Student haut 1.63774 =LOI.STUDENT.INVERSE(1-B29;ARRONDI.SUP(B28;0))
32 Facteur de couverture (kIT) 1.83676 =B30-(B30-B31)*(B28-ARRONDI.INF(B28;0))
33 Ecart-type de tolrance (sIT) 0.01828 =B22*RACINE(1+1/(B9*B27))
34 Intervalle de tolrance
35 Limite intervalle tolrance basse 0.37516 =B19-B33*B32
36 Limite intervalle tolrance haute 0.44230 =B19+B33*B32
Afin de faciliter la mise en uvre, les principales formules, telles que saisies, sont indiques
dans la colonne C et se rfrent au contenu de la colonne B.
En ce qui concerne ltalonnage et le calcul des modles, il suffit dutiliser la fonction de
rgression DROITEREG pour trouver les paramtres. Lorganisation de la feuille doit tenir
compte du fait quon calcule un modle par srie.
Principe
Sur la figure 1 on note quil existe une concentration partir de laquelle au moins une des
bornes de lintervalle de tolrance coupe une des limites dacceptabilit. Au-dessus de ce
seuil, lanalyste peut garantir que la mthode produira une proportion de rsultats au moins
gale . Selon les rgles de dcision proposes, la mthode est donc valide puisquelle
permet de garantir les objectifs fixs. En revanche, en dessous de ce seuil, il nest plus
Les formules utilises sont indiques dans les colonnes D et E. Finalement, la LQ vaut 0,434 mg/l.
Bien sr, la limite de quantification dpend des valeurs choisies des limites dacceptabilit et
de la proportion . Si ces deux critres nont pas pu tre fixs lavance par manque de
consignes ou de rglementation, mais si une LQ limite est attendue, par exemple pour un
contrle rglementaire, il est possible de les calculer en retour en se fixant une LQ prtablie.
Ce mode de calcul inverse se dduit directement des formules fournies.
Niveau B Niveau A
Z% Niveau B
Z
Niveau A
2.0 110%
109%
1.5
108%
1.0
107%
0.5
106%
LQ X LQ X
0.0 105%
0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5
4.5
4.0
Concentration retrouve Z (mg/L)
3.5
3.0
2.5
Z = 0.52X + 0.02
2.0
1.5
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
Valeur de rfrence X (mg/L)
Figure 3 : reprsentation des concentrations retrouves en fonction des valeurs de rfrence permettant
de mettre en vidence un effet de matrice et de calculer une droite moyenne de recouvrement
Une mthode plus exacte consiste reporter sur un graphique les concentrations retrouves en
fonction des valeurs de rfrence. La figure 3 montre comment mettre en vidence un effet de
matrice. On peut alors calculer, sous la forme dune droite de rgression au moindres carrs,
lquation dune droite dont la pente reprsente un taux moyen de recouvrement et lordonne
lorigine un biais systmatique. Elle permet ainsi de dfinir deux facteurs de correction
possibles :
linverse de la pente permet de corriger un taux de recouvrement non satisfaisant. Dans
lexemple on prendra 1/0,52 2.0 ;
lordonne lorigine permet de corriger un biais systmatique, par exemple gnr par
lintgration dun pic mal rsolu en chromatographie. Dans lexemple ce type de biais
nest pas important et ne sera pas pris en compte.
Le facteur de correction choisi sera alors inclus dans le mode opratoire de la mthode, par
exemple, dans la partie relative lexpression du rsultat. Afin de vrifier la pertinence du
facteur de correction choisi, il suffit de recalculer le profil dexactitude avec les donnes
corriges. Cependant, pour assurer lindpendance des mesures, il est prfrable dutiliser
deux chantillons de validation et de raliser lensemble des plans dexprience en parallle.
La figure 4 illustre leffet du facteur de correction 2,0. Lchantillon A dont les mesures nont
pas t corriges produit un profil dexactitude qui mne conclure que la mthode nest pas
valide. En revanche, les mesures faites sur un deuxime chantillon B puis corriges par 2,0,
conduisent la conclusion que la mthode est valide si les rsultats sont corrigs.
140%
chantillon B corrig
130%
120%
110%
100%
Recouvrement %
80%
60%
50%
40%
30%
Concentration thorique (mg/L)
Cette approche est volontairement empirique et ne fait pas appel des tests de signification de
la pente ou de lordonne lorigine.
2 1
sFI 1 sr2
u ( D) = I
J
Bibliographie
AFNOR Norme NF ISO 5725-2:1994, Exactitude (justesse et fidlit) des rsultats et
mthodes de mesure Partie 2 : Mthode de base pour la dtermination de la rptabilit
et de la reproductibilit d'une mthode de mesure normalise. Afnor, Paris
AFNOR Norme NF X 06-041, Application de la statistique Fidlit des mthodes d'essai
Dtermination de la rptabilit et de la reproductibilit d'une mthode d'essai
normalise par essais inter-laboratoires. Afnor, Paris
AFNOR Norme FD X07-021:1999 Normes fondamentales - Mtrologie et applications de la
statistique - Aide la dmarche pour l'estimation et l'utilisation de l'incertitude des mesures
et des rsultats d'essais, Afnor, Paris
AFNOR Norme FD V03-116:2001 Analyse des produits agricoles et alimentaires - Guide
d'application des donnes mtrologiques, Afnor, Paris
FD ISO/TS 21748 :2004. Lignes directrices relatives l'utilisation d'estimations de la
rptabilit, de la reproductibilit et de la justesse dans l'valuation de l'incertitude de
mesure, ISO, Genve
Feinberg M., (2009) Labo-Stat : Guide de validation des mthodes danalyse Lavoisier
Tec&Doc, Paris
Feinberg M. (2010a) Principes et vocabulaire in Validation des mthodes danalyse
quantitative par le profil dexactitude. Le Cahier des Techniques de lInra, numro
spcial : 13-26