Source: https://www.scribd.com/doc/92150291/AQAEXCEL-Validacion-de-Metodos
Timestamp: 2016-12-08 02:46:38+00:00

Document:
AQAEXCEL Validacion de Metodos
BrowseInterestsBiography & MemoirBusiness & LeadershipFiction & LiteraturePolitics & EconomyHealth & WellnessSociety & CultureHappiness & Self-HelpMystery, Thriller & CrimeHistoryYoung AdultBrowse byBooksAudiobooksArticlesSheet MusicBrowse allUploadSign inJoinAplicación de Microsoft Excel a la Química Analítica: validación demétodos analíticos José Marcos Jurado Departamento de Química Analítica 1 de abril de 2008 1. Estadística básica. 1.1. Concepto de población y muestra En estadística, se define población como el conjunto de individuos portadores de información del fenómeno que se estudia. Una muestra sería un conjunto de individuos representativos de dicha población, conjunto que permitiría estudiar el fenómeno sin pérdida significativa de información. Así, si se pretende analizar el contenido en sulfatos de un lago, la población sería el propio lago y la muestra sería una porción de agua representativa del mismo. De este modo, a partir del análisis de la muestra, debe ser posible conocer el contenido en sulfatos de la población. Cuando nos referimos a datos obtenidos experimentalmente en el laboratorio, la población se refiere a la totalidad de medidas posibles, mientras que la muestra sería un conjunto de estas medidas. 1.2. Media aritmética. Se define la media aritmética como la suma de una serie de individuos dividida por el número de individuos considerados. En el caso de una población, es decir, cuando realizásemos una serie infinita de medidas obtendríamos la media de la población, generalmente designada por µ. Cuando se realice una serie de medidas representativas de la población se debe hablar de media muestral. La media muestral se calcula como: 1
(1) siendo x
cada uno de los valores obtenidos y n el número de valores considerados. J. M. Jurado EXCEL para Química Analítica 2 1.3. Desviación estándar. La desviación estándar es un parámetro que da cuenta de la dispersión de los datos obtenidos. En el caso de una población hablaremos de la desviación estándar de la población, σ. Para una muestra se calcula con la siguiente fórmula: 2
(2) Al cuadrado de la desviación estándar se le denomina varianza (s
). A partir de la desviación estándar se puede también calcular el coeficiente de variación o desviación estándar relativa, que es una medida del error relativo y se usa para comparar resultados medidos en distintas escalas. 100
CV RSD
= = ⋅ (3) 1.4. Distribución de medidas repetidas. Cuando realizamos medidas de una población, en ausencia de errores sistemáticos, éstas pueden tomar cualquier valor en torno al valor verdadero o media de la población, dependiendo de la desviación estándar de la misma (del error aleatorio puro). Si se dispusiera de todas las medidas y se representasen las frecuencias con que aparece cada valor, frente a los propios valores, se podría ver la distribución de los datos en torno a la media de la población. En las condiciones mencionadas se obtiene una distribución normal o gaussiana como la de la figura. Figura 1. Distribución normal J. M. Jurado EXCEL para Química Analítica 3 La distribución normal viene descrita por la ecuación { }
exp ( ) / 2
y x µ σ
= − − (4) En una distribución normal de media µ y desviación estándar σ, aproximadamente el 68 % de los valores de la población caen dentro del intervalo µ ± 1σ , el 95 % dentro de µ ± 2σ y el 99.7 % dentro de µ ± 3σ. En realidad solo vamos a realizar una serie finita de medidas (replicados) a partir de las cuales se obtiene la media de la muestra y su desviación estándar. La media de la muestra será una estimación de µ y la desviación estándar de la muestra será una estimación de σ. Por lo tanto, la distribución de esta muestra será también normal y se cumplirán las mismas características anteriormente mencionadas. Gráficamente, donde antes obteníamos una curva ahora obtendremos un histograma. Si dispusiésemos de una serie de medias, cada una obtenida de una serie de medidas individuales, se observaría que estas medias se encuentran más agrupadas en torno a µ que las medidas individuales. La distribución de todas las medias muestrales posibles se denomina distribución muestral de la media o distribución en el muestreo de la media. Su media es la misma que la media de la población original. Su desviación estándar se denomina error estándar de la media (e.e.m.). Existe una relación matemática exacta entre el e.e.m. y la desviación estándar, σ, de la población: e.e.m.= σ/ n . A mayor n, menor será el e.e.m. y menor la dispersión de las medias muestrales en torno a µ. El e.e.m. nos da una medida de la variabilidad de x , además, aunque la distribución de las medidas no sea una distribución normal, la distribución de las medias muestrales tiende a serlo a medida que aumenta n. Esto se conoce como la teoría del Límite Central. Se asume que las medias de muestras muy pequeñas (<5) se distribuyen normalmente. J. M. Jurado EXCEL para Química Analítica 4 1.5. Limites de confianza de una medida. Es poco probable que obtengamos una media de la muestra exactamente igual a la de la población. Generalmente se tiene un intervalo centrado en la media de la muestra que contendrá al valor verdadero con un determinado nivel de confianza. La amplitud de este intervalo depende de dos factores: a) La precisión de los valores individuales, que dependen de la desviación estándar de la población. b) El número de medidas de la muestra, pues el hecho de repetir las medidas implica tener más confianza en el resultado. El límite de confianza de una media viene dado por la expresión: / x zs n ± (5) Donde, en ausencia de errores sistemáticos: / / x zs n x zs n µ − < < + (6) Donde z=1.96 para un 95% de nivel de confianza, 2.58 para el 99% y 2.97 para el 99.7%. 1.6. Mediana Para un conjunto de n datos ordenados de menor a mayor, la mediana es el valor de la observación que ocupa el lugar (n+1)/2 si n es impar o bien la media de los valores que ocupan los lugares n/2 y (n/2+1) si n es par. 1.7. Calibrado. Supongamos que pretendemos realizar la determinación de un determinado analito mediante un método instrumental, de modo que medimos una propiedad física de la muestra que sea directamente proporcional a la concentración de analito. En principio, a partir de la medida de esta propiedad no podremos calcular la concentración de analito. Para poder hacer esto es preciso establecer una relación entre las señales medidas y una serie de patrones de concentración conocida. Esto es lo que se denomina realizar un calibrado. La relación señal-concentración vendrá dada por una función matemática a partir de la cual se puede interpolar el valor de señal obtenido para una muestra y calcular así la concentración de analito en la misma. Generalmente se procura que la relación sea lineal, lo que simplifica mucho los cálculos a realizar. J. M. Jurado EXCEL para Química Analítica 5 Supongamos que disponemos de una serie de patrones de concentración x
2, …x
que, una vez medidos por una determinada técnica presentan unas señales analíticas y
2, …y
, … y
n. Si representamos las señales frente a las concentraciones obtenemos la denominada curva de calibración. En el caso más sencillo, como ya se ha dicho, tendremos una recta de calibrado de ecuación y = bx + a, donde b es la pendiente y a la ordenada en el origen. Para calcular el valor de estos dos parámetros se utiliza el método de los mínimos cuadrados, que consiste en ajustar estos parámetros para minimizar la suma del cuadrado de los residuales: 2
S y a bx
(7) Para poder aplicar este método se debe considerar que no existe error en la variable x, siendo este considerado sólo en la variable y por mediación de la varianza de repetición pura (obtenida a partir de replicados). Además, dicha varianza debe permanecer constante en los distintos niveles (homocedasticidad) y los valores de y deben estar normalmente distribuidos en cada nivel. Esto se puede observar en la gráfica de residuales. J. M. Jurado EXCEL para Química Analítica 6 Loas ecuaciones para calcular la pendiente y ordenada son: ( )( ) { }
(8) a y bx = − (9) Los errores de estos parámetros se calculan a partir de la varianza de regresión, que estima los errores aleatorios en la dirección y. ( )
(10) Donde ˆ
y a bx = + (11) A partir de esta se tiene la desviación estándar de la pendiente y de la ordenada: /
(13) El último parámetro que vamos a considerar es el coeficiente de correlación, r, que da cuenta de la bondad del ajuste de la curva a los datos experimentales. Su valor oscila entre 0 y 1, valores de r próximos a 1 implican un mejor ajuste. J. M. Jurado EXCEL para Química Analítica 7 ( )( ) { }
(14) Si tenemos una muestra de señal y
, se puede calcular la concentración correspondiente a la misma como: 0
= (15) Y el error de esta concentración viene dado por: ( )
(16) 1.8. Algunas fórmulas en EXCEL. Parámetro o función Fórmula Media =PROMEDIO() Desviación estándar =DESVEST() Varianza =VAR() Raíz cuadrada =RAIZ() Número de datos =CONTAR() Mediana =MEDIANA() Pendiente =PENDIENTE() Ordenada en el origen =INTERSECCION.EJE() Coeficiente de correlación =COEF.DE.CORREL() Estimación lineal =ESTIMACION.LINEAL() La última es muy útil cuando se quieren obtener los errores asociados a la pendiente y la ordenada en el origen. Se explicará mas adelante. J. M. Jurado EXCEL para Química Analítica 8 2. Validación de métodos analíticos. Un método analítico es un procedimiento que obtiene unos resultados a partir del análisis químico de muestras. Antes de ser usado, es necesario comprobar que dicho método cumple unos determinados requisitos que dependerán de la aplicación que se le quiere dar. Este proceso de verificación se conoce como validación y debe establecer las características de funcionamiento y las limitaciones del método, así como la identificar las influencias que pueden cambiar dichas características [1]. A continuación se detallan los distintos parámetros relacionados con la calidad de un método analítico que deben ser objeto de estudio en el proceso de validación. 2.1. Selectividad y especificidad. En todo método analítico hay un paso fundamental en el que se produce la medida de una señal relacionada con la concentración (análisis cuantitativo) o la presencia (análisis cualitativo) de un determinado analito. A veces, la presencia de otros compuestos, química o físicamente similares al analito, pueden contribuir a dicha señal, pudiendo causar un error sistemático en la determinación del analito. Estos compuestos pasan a denominarse interferencias. Es por tanto necesario asegurarse de donde procede la señal que se está midiendo. De aquí surge en concepto de selectividad, entendiéndose esta como el parámetro que define de forma cualitativa la extensión en que otras sustancias interfieren en la determinación de un analito usando un método determinado. Al máximo de selectividad se le denomina especificidad, es decir, un método será específico cuando la señal medida solo proviene del analito de interés [2]. Un estudio de selectividad implica investigar la capacidad del método para medir al analito en presencia de interferencias añadidas a la muestra [1], no siendo esto objeto del presente texto. Hasta ahora, solo hemos considerado las interferencias que pueden producir un error sistemático constante, es decir, sustancias que contribuyen a la señal medida con un valor determinado. Pero existen otro tipo de interferencias, relacionadas con efectos inespecíficos, que pueden producir errores proporcionales. Estas son las designadas como efecto matriz y, aunque no existe ninguna aplicación en Excel específica para cuantificar este efecto, hablaremos de el más adelante. 2.2. Exactitud. J. M. Jurado EXCEL para Química Analítica 9 La exactitud da cuenta de la proximidad entre el resultado obtenido por un método y el valor “real”. Como se trata de un parámetro cualitativo (un método será muy o poco exacto) para describir el error asociado a un resultado, hace falta un parámetro que permita medir la exactitud cuantitativamente. Este parámetro es el error que, según la guía para validación de métodos de la EURACHEM [1], se divide en dos componentes, la veracidad y la precisión. La veracidad es un parámetro que expresa la proximidad de la media de una serie de resultados obtenidos con el método al valor real. Generalmente se expresa en términos de error o bias, definido como la diferencia entre el resultado de medida y el valor real [2]. La precisión da cuenta de la proximidad entre sí de esta serie de resultados. Si expresamos la exactitud en función del error E
, para una determinada medida i: ) ( ) ( µ µ − + − = − = x x x x E
(17) En el caso de tener un replicado consistente en n medidas, el error cuadrático se calcularía como: 2 2
) ( ) ( µ − + − =
x n x x E
(18) Si dividimos por n, obtenemos el error cuadrático medio: 2
(19) Si n es suficientemente grande, el primer término tiende a la varianza de la medida (dividiendo entre n-1 en vez de n) y podría expresarse: 2 2 2
( ) E s x µ = + − (20) Con lo cual vemos como el error cuadrático medio correspondiente a la exactitud se equipara a la suma de la varianza de repetición pura (error aleatorio) más el cuadrado de la veracidad (error sistemático). La exactitud es máxima cuando el sesgo es nulo y la dispersión aleatoria es mínima. J. M. Jurado EXCEL para Química Analítica 10 2.3. Veracidad. Como se ha dicho anteriormente, la veracidad da cuenta de lo cercano del resultado medio obtenido con el método y el valor real. Para asegurar la veracidad es necesario compara el valor medio medido con un valor de referencia (que será el valor real o el considerado como tal). Esto se realiza, generalmente comparando los resultados obtenidos por el método con el valor de un material de referencia certificado o comparando con un método ya validado (que suele haber sido sometido a un ejercicio de intercomparación). En caso de no existir materiales de referencia o no tener posibilidad de aplicar el método de referencia, se puede simular la matriz de la muestra añadiendo al analito en una cantidad conocida y efectuar un ensayo de recuperación. Como última instancia, se pueden realizar dichos ensayos de recuperación añadiendo al analito sobre la propia muestra [3]. Como se explicará más adelante, para compara dos valores será necesario conocer la precisión con la que este ha sido obtenido, por ello la exactitud se evalúa a partir de la veracidad en conjunto con la precisión. 2.4. Precisión. Cuando realizamos una medida se deben obtener replicados, a partir de los cuales se puede observar el grado de concordancia que estos resultados tienen entre sí. Como se dijo anteriormente, este grado de concordancia se denomina precisión. La precisión se divide en dos parámetros repetitividad y reproducibilidad. 2.4.1. Repetitividad. Se define la repetitividad como la desviación estándar obtenida al analizar una misma muestra varias veces, en un periodo de tiempo corto, sin cambiar de equipo de medida, reactivos o analista. 2.4.2. Reproducibilidad. La reproducibilidad se define como la desviación estándar obtenida al analizar varias veces la muestra en días distintos, pudiendo variar condiciones tales como el equipo, reactivos o analistas. Se habla de reproducibilidad interlaboratorio cuando las medidas se realizan en laboratorios distintos (ensayos de intercomparación). La reproducibilidad intralaboratorio se realiza en un solo laboratorio, preparando los patrones de calibración cada vez, en distintos días y cambiando en la medida de lo posible de equipo y analista. A esta reproducibilidad intralaboratorio se la conoce como precisión intermedia. J. M. Jurado EXCEL para Química Analítica 11 2.5. Sensibilidad. Según la IUPAC la sensibilidad se define como el cociente entre la señal medida y la concentración de analito. Cuando trabajamos con curvas de calibrado, la sensibilidad es la pendiente de la curva a una concentración dada. En caso de ser una recta, coincide con la pendiente de la misma. 2.6. Límites de detección y cuantificación. Estos parámetros están relacionados con la sensibilidad del método. El límite de detección (LOD) se define como la cantidad de analito presente que se deriva de la menor señal analítica que puede detectarse con certeza razonable (Y
) [4]. Si analizamos un blanco de la muestra varias veces y obtuviéramos la señal media (
Y ) y su desviación estándar, 0
s , las señales mayores a 0
Y Y s = ± caerían fuera de la gausiana del blanco analítico y podríamos asignarlas a la presencia de analito. Si tenemos una recta de calibrado de ecuación y = bx + a, el LOD se calcularía como: LOD
= (21) Como la señal del blanco debe coincidir con la ordenada del blanco analítico (puesto que es la señal que se obtendría a una concentración de analito igual a cero), el LOD se calculará como: 3
= (22) siendo s
la desviación estándar de la ordenada. El límite de cuantificación (LOQ) se define como la cantidad de analito presente derivada de la menor señal analítica que puede ser medida con precisión y exactitud. De manera análoga al LOD, la IUPAC propone como LOQ al valor 0
LOQ Y
Y Y s = ± , que lleva a una concentración, a partir de la recta de calibrado: 10
= (23) 2.7. Linealidad e intervalo dinámico La linealidad a veces se examina a partir del coeficiente de correlación, considerando lineal el calibrado si r > 0.999. Otro método se basa en que la desviación estándar J. M. Jurado EXCEL para Química Analítica 12 relativa de la pendiente no supere el 5%. Esta es la conocida como linealidad online. El porcentaje de linealidad se calcula como [5]: % 100 1
(24) Estos parámetros no son del todo fidedignos, siendo lo correcto realizar un análisis estadístico de los residuales mediante el análisis de la varianza (ANOVA) o la aproximación práctica de Huber [6]. Ambos procedimientos los explicaremos detenidamente cuando estudiemos la aplicación de EXCEL. El intervalo dinámico o rango lineal es el rango de concentraciones donde existe una relación lineal entre estas y las señales analíticas. En este intervalo, los valores deben presentar un nivel aceptable de precisión y exactitud, con lo cual, no pueden ser inferiores al límite de cuantificación. 2.8. Robustez Es la inercia que presenta un método analítico a modificar su señal cuando tienen lugar pequeños cambios en las condiciones ambientales u operativas, que se consideran las variables que gobiernan el experimento de medida. Para el estudio de robustez, estas condiciones se modifican ligeramente, siguiendo un determinado diseño experimental, y se comprueba su influencia en la señal analítica. Con ello, se puede concluir que variables son más significativas a la hora de realizar las medidas y por tanto, deben ser mejor controladas. 2.9. Ámbito de validación de métodos analíticos Como ya hemos dicho, el tipo de validación que debemos realizar va a depender de los requisitos de nuestro método. Parámetro Análisis Cualitativo Análisis Cuantitativo Análisis de trazas Especificidad X X X Precisión X X Veracidad X X Linealidad X X Rango X LOD/LOQ LOD LOQ Robustez X X X J. M. Jurado EXCEL para Química Analítica 13 3. EXCEL como herramienta en el proceso de validación. Las utilidades estadísticas de EXCEL que emplearemos en validación serán: Prueba F de Snedecor-Fisher Prueba t de Student ANOVA de un factor Estadística descriptiva (media, mediana, varianza…) Regresión 3.1. Herramientas para análisis. Aunque todo esto se puede hacer de forma manual, es decir, introduciendo formulas en la hoja de cálculo, EXCEL incluye un paquete de herramientas para realizar este tipo de análisis de manera sencilla. El primer paso será comprobar que está instalado o activado. Para ello hay que abrir el menú Herramientas>Complementos y ver si aparece la opción Herramientas para análisis. Si no aparece debemos activar el paquete en C:\Archivos de programa\Microsoft Office\OFFICE11\Macros\Análisis\ANALYS32. En este momento dispondremos de la herramienta análisis de datos. J. M. Jurado EXCEL para Química Analítica 14 3.2. Prueba F de Snedecor- Fisher. En una prueba estadística empleada para la comparación de varianzas. Sean dos conjuntos de datos: {x
} con varianza
s y tamaño n
} con varianza 2
Consideramos la hipótesis nula que dice que las varianzas no presentan diferencias significativas entre sí. H
σ σ = . Si 2 2
s s ≥ , 2
= ≥ . Se compara con F
(0.05, n
-1). Si F<F
admitimos H
0. Los valores de F para una cola se pueden obtener usando la función =DISTR.F.INV(p;n
-1;n
-1), siendo p la probabilidad (normalmente 0.05, es decir 95% de nivel de confianza) 3.2.1. Aplicaciones - Comprobar que dos métodos son igual de precisos. - Comprobar que las varianzas son iguales para ver que tipo de test de comparación de medias usar posteriormente. - Comprobar que en la reproducibilidad no interviene factores distintos a los encontrados en el estudio de repetitividad. 3.2.2. Resolución con EXCEL Supongamos que hemos realizado una serie de medidas de una muestra con dos métodos analíticos y queremos comprobar si ambos son igual de precisos. Lo primero es tener los datos en EXCEL. J. M. Jurado EXCEL para Química Analítica 15 Del menú Herramientas>Análisis de datos usamos la herramienta “Prueba F para varianza de dos muestras” Pulsamos en la flecha roja del cuadro rango para la variable 1 y seleccionamos el rango A2:A7. En rango para la variable 2 seleccionamos el rango B2:B9. Si marcamos Rótulos, el programa considera que en la primera celda del rango seleccionado se encuentra el nombre del conjunto de datos (deberíamos haber seleccionado A1:A7 y B1:B9). Alfa es el indica la probabilidad o nivel de confianza que usaremos (0.05 es un 95%). En rango de salida se indica la primera celda del rango donde se colocarán los resultados. Es aconsejable que sea en un extremo para que no sobrescriba nada. Otra opción es hacerlo en una hoja nueva. Supongamos que hemos elegido la celda D1 como rango de salida. El resultado será: J. M. Jurado EXCEL para Química Analítica 16 Se comprueba que F < F
, ambos métodos son igual de precisos. Si necesitásemos el valor de F crítico para comparar las varianzas a manos se puede obtener de tablas o bien con la formula =DISTR.F.INV(0.05,5,7), en este caso. 3.3. Prueba t de Student. 3.3.1. Comparación de valores medios de dos conjuntos de datos. Dados dos conjuntos de datos: Conjunto Media Varianza Número de datos 1 1
n Para decidir si las medias muestrales 1
x difieren significativamente se toma la hipótesis nula H
µ µ = . Primero se comprueba que 2 2
s s = mediante la prueba F de Fisher. Así tenemos dos posibilidades: a) Si 2 2
s s = se calcula una varianza promedio ( ) ( )
(25) Se calcula el estadístico t: 1 2
(26) , que tendrá n
– 2 grados de libertad. Este valor de t se compara con el valor tabulado t(0.05, n
- 2). Si t < t
se acepta la hipótesis nula, es decir, las dos medias son estadísticamente iguales. J. M. Jurado EXCEL para Química Analítica 17 b) Si 2 2
s s ≠ se calcula el estadístico t como: 1 2
(27) Ahora los grados de libertad se deben calcular como: ( ) ( )
(28) El valor de t se compara con t(0.05,
ν ). Los valores de t se pueden obtener con la función =DISTR.T.INV(p;g.l.). Estos valores son adecuados para contraste de dos colas, para una cola se toma probabilidad 2*p (siendo p el nivel de confianza deseado). 3.3.1.1. Aplicaciones. Comparación del valor medio obtenido con para un material de referencia aplicando un determinado método y el valor de un certificado. Comparación de los resultados obtenidos con un método en desarrollo y uno de referencia. 3.3.1.2. Resolución con EXCEL. a) Conjunto de datos con varianzas estadísticamente similares. Supongamos que partimos del ejemplo anterior, en el que mediante el test de Fisher se dedujo que las varianzas eran similares. En Herramientas>Análisis de datos se escoge la opción “Prueba t para dos muestras suponiendo varianzas iguales” J. M. Jurado EXCEL para Química Analítica 18 La forma de escoger los rangos, nivel de confianza y rango de salida son similares a la prueba F. Nótese que es posible tomar distintas hipótesis nulas en la casilla Diferencia hipotética entre las medias. El resultado es: Como en valor absoluto el t calculado es menor que el t crítico (contraste de dos colas), la hipótesis nula se toma como cierta y las dos medias muestrales no difieren significativamente. El hecho de usar un contraste de dos colas implica que no nos interesa el sentido en que difieren las medias (es decir si la diferencia es positiva o negativa). Si pretendemos ver solo si un valor aumenta en relación al otro se usa un contraste de una cola. b) conjunto de datos con varianzas distintas. Considérese el conjunto de datos: J. M. Jurado EXCEL para Química Analítica 19 Al aplicar el test F: Las varianzas difieren significativamente. A la hora de realizar el contrate t, se debe elegir “Prueba t para dos muestras suponiendo varianzas desiguales”. Los rangos, probabilidad e hipótesis nula se introducen de igual modo que en el ejemplo anterior. El resultado es: Comparando con el valor de t crítico para dos colas, las medias no difieren significativamente. Nótese que los grados de libertad efectivos son 5 (en realidad 5.2868, pero se ajusta al entero más próximo). 3.3.2. Prueba t para valores emparejados. Supongamos que tenemos un conjunto de muestras de valores 1
x x x . Nótese que se trata de n muestras distintas, por lo que las cantidades 1
x x x ≠ ≠ ≠ ≠ . Se determinan los valores mediante dos métodos y se pretende comprobar que no existen diferencias significativas entre los resultados obtenidos por ambos. Al aplicar cada método se obtienen los valores 11 1 1
x x x con el método 1 y 21 2 2
x x x con el J. M. Jurado EXCEL para Química Analítica 20 método 2. Se calculan las diferencias 1 2 i i i
d x x = − , el valor medio de todas las diferencias y su desviación estándar: n
(30) Si los métodos son equivalentes y levan al mismo resultado las diferencias d
deben ser nulas. Esto es lo mismo que considerar la hipótesis nula H
d = 0. Se calcula el estadístico t: /
t d n s = (31) El número de grados de libertad de t es n-1. Se compara este valor con t(0.05,n-1). 3.3.2.1. Aplicaciones. Se aplica en el caso de que las muestras analizadas por los dos métodos no sean grandes y sea necesario utilizar distintas muestras en el estudio. Es decir, se analizan distintas muestras por ambos métodos y no una misma muestra replicada por ambos métodos. 3.3.2.2. Resolución con EXCEL Consideramos el ejemplo de la determinación de paracetamol en diez pastillas de distintos lotes mediante un método espectrofotométrico y otro de espectroscopía IR [7]. Se utiliza la “Prueba t para dos muestras emparejadas” del menú Herramientas>Análisis de datos. El formulario de entrada de datos es similar a los de los ejemplos anteriores. El resultado es: J. M. Jurado EXCEL para Química Analítica 21 Como el valor experimental (0.88) es menor que el crítico (2.26) los métodos no proporcionan resultados significativamente diferentes. 3.4. ANOVA de un factor. El análisis de la varianza (ANOVA) es un método que se emplea para comparar resultados obtenidos por distintos métodos, laboratorios, analistas, etc, cuando el número de medias obtenidas es superior a dos. También permite separar las contribuciones de uno o más factores ala varianza global del sistema. En el caso de realiza medidas replicadas variando uno de los factores anteriormente mencionados siempre existen dos fuentes de error, la primera, el error aleatorio de medida, la segunda los errores debidos al cambio de método, laboratorio, analista, etc. Mediante el ANOVA se puede controlar el error introducido por esta segunda fuente, con lo que se habla de ANOVA de un factor. Si tuviésemos un caso en que el análisis se realiza por distintos laboratorios utilizando distintos métodos, deberíamos controlar dos factores y se utilizaría un ANOVA de dos factores. Consideremos que realizamos un mismo análisis por i métodos distintos (i varía de 1 a h) y realizamos j replicados (j varia de 1 a n) con cada método. La media obtenida para cada método será: n
(32) La media global obtenida por todos los métodos se calculará como: h
(33) J. M. Jurado EXCEL para Química Analítica 22 Media Método 1 11
x ... 1 j
x ... 1n
x Método 2 21
x ... 2 j
x ... 2n
x . . . . . . . . . . . . . . . . . . . . . . . . Método i 1 i
x ... ij
x ... in
x . . . . . . . . . . . . . . . . . . . . . . . . Método h 1 h
x ... hj
x ... hn
x Media global x Para cada método se tiene una varianza calculada como: ( )
(34) En el ANOVA la hipótesis nula es que todas las muestras se extraen de una población de media µ y varianza 2
σ . Esta varianza se puede estimar de dos formas, estudiando la variación dentro de las muestras y la variación entre muestras. La varianza debida al error puramente aleatorio, 2
s , también conocida como varianza dentro de las muestras (o del grupo), se calcula como promedio de las varianzas individuales 2
s : ( )
(35) La varianza debida al factor (en nuestro caso a los distintos métodos usados), conocida como varianza entre muestras (o entre grupos), es una estimación de la varianza de la población, 2
σ . Si todas las medidas se extraen de una población de varianza 2
σ , entonces sus medias (obtenidas variando el método empleado, es decir, variando un factor controlado) proceden de una población de varianza 2
/ n σ (Véase sección 1.4). En general, la varianza debida al factor se calcula como: J. M. Jurado EXCEL para Química Analítica 23 ( )
(36) Si la hipótesis nula es correcta, las dos estimaciones no deberían diferir significativamente, es decir, que la variación debida al factor no debe se significativamente distinta a la varianza de error puro. Esto es lo mismo que decir que el hecho de variar el factor (el método de medida en nuestro ejemplo) no introduce un error significativo en comparación con el puramente aleatorio. Se comprueba esto mediante un test de Fisher, calculando F como: 2
= (37) Se compara F con el valor tabulado F(0.05, h-1, h(n-1)) para contraste de una cola (puesto que nos interesa saber si la varianza del factor es mayor que la de error puro, no solo si difieren). Si el valor de F calculado es menor que el tabulado, no existen diferencias significativas entre las dos varianzas, es decir, el hecho de variar el factor (cambiar de método) no introduce un error significativo. Si la hipótesis nula es verdadera existe una tercera forma de estimar la varianza que consiste en tratar los datos como si fueran una muestra grande de varianza: ( )
(38) Pero esta fórmula no se usa en el contraste. Hasta ahora hemos supuesto que la hipótesis nula es verdadera. Pero puede que no lo sea, es decir, que F sea mayor que el valor de F crítico. Un resultado significativo en el ANOVA puede deberse a varias cosas: que una media difiera mucho a las demás, que todas las medias difieran entre sí o que existan varios grupos de medias distintas, etc. Para poder ver lo que ocurre se colocan las medias ordenadas de forma creciente y se comparan las diferencias de los valores adyacentes con una cantidad denominada la mínima diferencia significativa. Ésta se calcula como: ( 1)
PE h n
MDS s t
= ⋅ (39) J. M. Jurado EXCEL para Química Analítica 24 Otra posibilidad es comparar las medias dos a dos mediante el test de de la genuina diferencia significativa de Tukey. Si una vez realizados los cálculos de ANOVA se han encontrado diferencias significativas entre la varianza de un factor y la del error puro, debemos averiguar que nivel de ese factor es el responsable de los errores sistemáticos que producen dichas diferencias. Con este fin se realiza el ensayo Tukey-HSD (Honest Significant Difference) [8], por el que se comparan dos a dos las medias obtenidas a cada nivel del factor. Aquel nivel que presente un valor promedio significativamente diferente del resto será el responsable del sesgo detectado y se elimina para recalcular el ANOVA. Supongamos que queremos comparar dos niveles p y q del factor (método en nuestro caso). Para que no haya diferencia significativa entre los valores promedio de estos dos niveles, p
x , se debe demostrar que ambos se corresponden con las medias muestrales de una misma población, de distribución normal, de media µ y varianza σ
. De este modo, si p
x es la media muestral de una población de media µ
x es la media muestral de una población de media µ
, se debe cumplir, como hipótesis nula, que µ
. Se define el intervalo de confianza de Tukey como: , 2
p q p q PE
− ∈ − ± +
(40) Donde, α
ν , r
q es el valor tabulado de la distribución del rango estandarizado [9], que depende del numero de repeticiones, r, y de los grados de libertad utilizados en la desviación estándar que se usa en el calculo de este intervalo, ν , así como del nivel de confianza, α . Teniendo en cuenta que, la varianza utilizada en la expresión anterior es una estimación de la de la población y debe considerar todas las muestras, se sustituye ésta por la varianza del error puro. Además, si en nuestro caso n
q = n, el intervalo de confianza de Tukey queda: J. M. Jurado EXCEL para Química Analítica 25 ,
µ µ − ∈ − ± (41) Como 0 = −
µ µ , el intervalo de confianza de Tukey debe contener el cero, si no es así, se rechaza la hipótesis nula. , ,
p q PE p q PE
x x s x x s
− − < < − + (42) En la práctica, se parte de la hipótesis nula, µ
q, y se calcula un valor de q experimental según: exp
= (43) Se compara este valor experimental de q con el valor tabulado [9]. Si q
se rechazaría la hipótesis nula, con lo que las dos medias presentarían diferencias significativas entre sí. El nivel que presente diferencias significativas con el resto, se considera como causa del efecto significativo en la varianza total. Una tabla de valores de la distribución del rango estandarizado se incluye en el ANEXO 1. 3.4.1. Aplicaciones. El análisis de la varianza se aplica, como se dijo anteriormente, en la comparación de varios métodos, laboratorios, analistas. etc. Es muy útil para detectar resultados anómalos en los ensayos de intercomparación empleados para certificar materiales de referencia o desarrollar métodos de referencia. También se utiliza para calcular la precisión intermedia de un método cuando se analizan muestras de referencia en distintos días. Este ejemplo lo explicaremos más adelante. Otra aplicación tiene que ver con la regresión lineal, puesto que EXCEL realiza un análisis de la varianza de los residuales. Esto lo veremos en el apartado correspondiente a la estimación lineal. Ensayos de recuperación. J. M. Jurado EXCEL para Química Analítica 26 3.4.2. Resolución con EXCEL. Supongamos que realizamos medidas de fluorescencia a muestras de quinina almacenadas de cuatro formas distintas (ejemplo adaptado de [7]): A) disoluciones recién preparadas, B) almacenada en oscuridad, C) con luz tenue y D) luz brillante. Los datos los ordenamos como se observa a continuación: En el menú Herramientas>Análisis de datos seleccionamos “Análisis de varianza de un factor”. El formulario de entradas de datos es el siguiente: Como el rango de entrada seleccionamos A1:D4, señalando la casilla de rótulos para que tome de la primera fila el nombre de los grupos. Hay que cerciorarse que tenemos marcado que los datos se agrupan por columnas. El rango de salida funciona como anteriormente se ha descrito. El resultado final es el siguiente: La columna de promedio de cuadrados indica el valor de varianza de error puro (dentro de los grupos) y la debida al cambio en el factor (entre grupos), que en este caso es el cambio en condiciones de almacenamiento. Como se observa el valor de F calculado (20.66) es mayor que el crítico para un nivel de confianza del 95 % y 3 y 8 grados de J. M. Jurado EXCEL para Química Analítica 27 libertad (4.066). Existen entonces diferencias significativas entre las medias obtenidas tras las distintas formas de almacenamiento de las muestras. Los test de minima diferencia significativa o de Tukey deben realizarse manualmente (EXCEL no implementa los mismos). 3.4.3. La precisión a partir del análisis de la varianza. Hasta ahora se ha aplicado el ANOVA para la comparación de múltiples medias. En este caso se va ha emplear para separar y estimar distintas fuentes de variación. Supongamos que realizamos un estudio de reproducibilidad a un solo nivel de concentración. Tendremos dos fuentes de variación, la debida al error puramente aleatorio, es decir, el error de replicación, y el debido al cambio de condiciones (día, reactivos, analista, laboratorio...) que implica el hecho de estudiar la reproducibilidad. Lo ideal sería añadir otro factor, realizando el estudio en distintos niveles de concentración, pero entonces realizaríamos un ANOVA de dos factores. Es más, al prevalecer el factor de las condiciones de trabajo sobre los niveles de concentración deberíamos realizar un ANOVA de dos factores anidados. Consideremos que realizamos n replicados durante h dias (reproducibilidad entre dias, minimalista). La varianza de repetitividad o de error puro vendrá dada por la ecuación 35, mientras que la debida al factor vendrá dada por la ecuación 36. Puesto que en la ecuación 36 están incluidas las medias de cada grupo, i
x , se estaría considerando el efecto de la repetitividad dentro del efecto del factor. Por ello, se debe calcular una varianza neta “entre grupos” que dependa sólo de la varianza entre las medias de los grupos y la media global. La varianza neta “entre grupos se calcula como: 2 2
2 Factor PE
= (44) La reproducibilidad se calcula finalmente como la suma de la repetitividad (error puro) y la varianza neta debida al factor (neta entre grupos): 2 2 2
R PE B
s s s = + (45) Para comparar la tolerancia correspondiente a esta precisión intermedia se puede usar la conocida como Trompeta de Horwitz [11]. J. M. Jurado EXCEL para Química Analítica 28 Horwitz dedujo una expresión para predecir el valor esperado para la desviación estándar relativa para la precisión intermedia (o interlaboratorio) a partir de la concentración de analito, c (en tanto por 1). (1 0.5log )
= (46) El valor de RSD obtenido se compara con el predicho por Horwitz (RSD
) mediante el parámetro Horrat [12]. R
= (47) Si el valor del parámetro Horrat es igual o menor a 2 se puede decir que el método tiene valores aceptables de precisión intermedia. Además de comparar con los valores de Horwitz, se puede también comparar con los valores establecidos por la AOAC [13] (Asociación Oficial de Químicos Analíticos). Estos valores y los obtenidos por Horwitz se presentan en la siguiente tabla. J. M. Jurado EXCEL para Química Analítica 29 Analyte (%) Analyte fraction Unit Horwitz %RSD AOAC PVM %RSD 100 1 100% 2 1.3 10 10
10% 2.8 1.8 1 10
1% 4 2.7 0.1 10
0.1% 5.7 3.7 0.01 10
100 ppm 8 5.3 0.001 10
10 ppm 11.3 7.3 0.0001 10
1 ppm 16 11 0.00001 10
100 ppb 22.6 15 0.000001 10
10 ppb 32 21 0.0000001 10
1 ppb 45.3 30 Algunos requisitos prácticos relacionados con los estudios interlaboratorio son los siguientes [14]: - Repetitividad: RSD
r = 0.5-0.6 * RSD
- Reproducibilidad: RSD
R = 0.5-2 * RSD
Para validaciones intralaboratorio [15]: - Repetitividad: RSD
r = 0.2-0.3 * RSD
- Reproducibilidad (precisión intermedia): RSD
R = 0.2-1 * RSD
3.4.4. Precisión y veracidad a partir del ANOVA. A partir de un mismo diseño de ANOVA es posible estudiar tanto la precisión (veáse el apartado anterior) como la veracidad. Simplemente hay que considerar que cada medida x
tendrá un sesgo ij ij
x T δ = − (48) Siendo T el valor certificado de un CRM (pues el valor verdadero µ no se conoce). El sesgo se puede entonces calcular como: 1
δ = − = −
(49) La incertidumbre de este sesgo se puede calcular como [16]: J. M. Jurado EXCEL para Química Analítica 30 2 2
= (50) Donde 2
s es la varianza de reproducibilidad (puesto que consideramos un ensayo interlaboratorio), que puede ser sustituida por la varianza de precisión intermedia 2
s en caso de estar realizando la validación intralaboratorio. Se calcula a partir de la ecuación 43. La varianza de repetitividad 2
s no es otra que la varianza de error puro 2
s obtenida en el ANOVA. Se considera que hay ausencia de sesgo para 2
< . 3.4.5. Ensayos de recuperación. Al desarrollar un método totalmente nuevo en una matriz no muy estudiada puede ocurrir que no existan materiales de referencia ni otros métodos validados para comparar el método en desarrollo. Entonces se hace necesario el empleo de ensayos de recuperación [3]. Estos consisten en añadir una cantidad de analito conocida a una matriz que simule a la muestra, analizarla y ver la cantidad de analito que se recupera. Se deben ensayar al menos tres niveles de concentración, uno bajo, otro medio y otro alto. A veces, cuando no es posible preparar una matriz similar a la muestra se realizan los ensayos sobre la propia muestra. Se calcula la recuperación como: obs native
= (51) Donde C
es la concentración medida con nuestro método, C
es la concentración que tiene la muestra (en el caso de ser una matriz preparada debe ser cero) y C
spike es la cantidad que hemos añadido. Si realizamos j replicados en cada uno de los i niveles podemos ver si las recuperaciones son equivalentes en cada uno de los niveles estudiados. 3.5. Regresión Ya se ha explicado la forma de calcular los coeficientes de regresión, cuando la relación entre las variables x e y es lineal. Por todos es conocido que si representamos los datos en Excel y seleccionamos los puntos, pulsando en el botón secundario del ratón, en “Agregar línea de tendencia”, como se indica en la figura, se puede obtener la recta de regresión con su ecuación correspondiente. Los pasos a seguir son los siguientes: J. M. Jurado EXCEL para Química Analítica 31 Se selecciona el tipo de estimación y se pulsa en la pestaña “Opciones”. Marcar “Presentar ecuación en el gráfico” y “Presentar valor de R cuadrado en el gráfico”. Pulsar en aceptar. J. M. Jurado EXCEL para Química Analítica 32 El resultado es el siguiente: Ahora vamos a explicar la herramienta “Regresión” del menú Herramientas>Análisis de datos. Partimos de los datos: El formulario de la herramienta Regresión es el siguiente: J. M. Jurado EXCEL para Química Analítica 33 Loas casillas de rangos de entrada y salida y la de rótulos ya son conocidas por nosotros. La casilla “Constante igual a cero” se utiliza si queremos forzar a la curva a pasar por el cero. Nosotros no la forzaremos. Esta herramienta ofrece siempre el intervalo de confianza de la pendiente y la ordenada en el origen al 95% de nivel de confianza, si se desea otro nivel de confianza se puede activar la casilla “Nivel de confianza, por ejemplo para el 90%. En el cuadro “Residuales” se puede especificar si se quieren tener los residuales, residuales estandarizados, el gráfico de residuales y el de la propia curva ajustada. Los residuos estándares se obtienen como: ,
= (52) Donde Res
es el i-ésimo residual estandarizado, Res
el i-ésimo residual, Res
es la media de los i residuales y Res
su desviación estándar. Finalmente se puede obtener la curva de probabilidad normal. Este último presenta los valores de y ordenados de menor a mayor indicando el percentil al que pertenece. Si los puntos parecen ajustarse a una línea recta, puede decirse que parece indicar que los datos provienen de una distribución normal, pero téngase en cuenta que en algunos casos, aunque los puntos se ajusten a una línea recta puede que los datos no provengan de una distribución normal; por ello se recomienda utilizar métodos objetivos. Seleccionamos B1:B9 para el rango y y A1:A9 para el x. Seleccionamos Rótulos, Residuos, Gráfico de residuales y Curva de regresión ajustada. El resultado será el siguiente: J. M. Jurado EXCEL para Química Analítica 34 En las estadísticas de regresión se tiene el coeficiente de correlación y el error típico (s
). El coeficiente de determinación R
, que coincide con el cuadrado del coeficiente de correlación, se obtiene a partir de la suma de cuadrados de los residuales (SCRes) y la suma de cuadrados total (SCT): 2
ˆ Re ( )
SC s y y = −
SCT y y = −
= − (55) El ajustado se calcula a partir de las varianzas de los residuales y total, también llamados cuadrados medios o promedio, CMRes y CMT. 2
= − (58) Cuanto más próximo a 1 sea el coeficiente de correlación mejor será el ajuste lineal. No es el criterio que tomaremos, trataremos este tema cuando hablemos de linealidad y rango lineal. J. M. Jurado EXCEL para Química Analítica 35 El análisis de la varianza aquí realizado se basa en que la suma de cuadrados totales se descompone en dos componentes, la suma de cuadrados de los residuales y la de regresión: 2 2 2
y y y y y y − = − + −
(59) Estas sumas de cuadrados aparecen en la tercera columna y dividiendo entre los grados de libertad de la segunda se obtienen las correspondientes varianzas (llamadas por EXCEL promedio de cuadrados). Excel solo calcula las varianzas de regresión y la de los residuales (cuarta columna). A partir de ellas calcula un Valor de F (quinta columna) y calcula una probabilidad para este valor (es lo que llama valor crítico de F). 2
(60) En general, si el valor de F es muy elevado, la suma de cuadrados de residuales será mucho menor que la de regresión, con lo que también será mucho menor que la suma de cuadrados totales. Entonces, según la expresión 55, r será próximo a 1. El hecho de que EXCEL trabaje con un valor de probabilidad es similar a comparar con el valor tabulado de F para 1 y n-2 grados de libertad, pero en vez de comparar, calcula directamente la probabilidad de que el F calculado sea mayor que el tabulado. Si la probabilidad es muy pequeña es poco probable que los datos lleven a un valor tan elevado de F, es decir, es poco probable que la suma de cuadrados totales sea mucho mayor que la de residuales y r sea próximo a uno. Esto quiere decir que el ajuste lineal no es fruto de la casualidad, porque es poco probable que los datos se organicen de este modo por casualidad. Como ya hemos dicho, nosotros estamos más acostumbrados a comparar el valor de F experimental con un valor tabulado de F para un nivel de confianza determinado (95 %, por ejemplo) teniendo en cuenta que los grados de libertad para la varianza de regresión es (n-1-n-2=1) y la de la varianza de los residuales n-2. Si F
>F(0.05, 1, n-2) la correlación lineal es significativa y no puede atribuirse a la casualidad. Los parámetros ajustados de pendiente y ordenada se obtienen en: J. M. Jurado EXCEL para Química Analítica 36 Intercepción se refiere a la ordenada y X a la pendiente. En la columna de error típico aparecen las desviaciones estándar de la ordenada y la pendiente. Se calcula un valor de t como el valor del parámetro dividido de su desviación estándar, comparando con t(0.05, n-2) grados de libertad. Igual que en el caso anterior se tiene la probabilidad de que éste t experimental sea mayor que el teórico. Los límites inferiores y superiores se calculan restando y sumando al valor del parámetro el valor obtenido al multiplicar su error por t(0.05,n-2). También se observan los valores de los residuales y los gráficos de residuales y de la curva ajustada. Observando los residuales podemos ver si: a) Si el modelo se ajusta bien a los datos. Los residuos se distribuyen hogéneamente. b) Si los datos son heterocedásticos, es decir, que la varianza no es constante en todo el rango de x. c) Si el modelo es inadecuado. Se observa cuando los residuales de valores próximos en x son todos del mismo signo, cambiando de signo en otro rango de x. J. M. Jurado EXCEL para Química Analítica 37 d) Si existen outliers. 3.5.1. La función =ESTIMACION.LINEAL() Esta función se utiliza de una forma especial, puesto que se trata de una formula matricial. En primer lugar hay que seleccionar un rango de 10 celdas tal y como se indica en la figura. En esta figura se ha escrito en las celdas adyacentes a las del rango el parámetro que la formula devolverá. La pendiente b y su desviación estándar sb, la ordenada a y su desviación estándar sa, el coeficiente de correlación al cuadrado r2 (para obtener r realizar la raíz cuadrada) y la desviación estándar de regresión sy/x son los parámetros fundamentales que se han explicado en la sección 1.7. El valor de F se corresponde con el anteriormente descrito y df son los grados de libertad para los residuales. El parámetro aquí llamado ssres es la suma de los residuales al cuadrado (Ecuación 7). En cuanto a ssreg, es la suma de cuadrados de regresión. El siguiente paso es seleccionar la fórmula ESTIMACION.LINEAL del menú “Insertar /Función” de EXCEL J. M. Jurado EXCEL para Química Analítica 38 En “Conocido_y” se debe seleccionar los valores de y, en “Conocido_x” los de x. En “Constante” se debe poner un 1 o VERDADERO para que calcule la ordenada en el origen. En “Estadistica” se pondrá un 1 o VERDADERO para que el calculo de ssreg y ssres sea de la forma anteriormente explicada. En vez de pulsar aceptar, sino CTRL+MAYÚS+ENTRAR, para que la fórmula de resultados en todas las celdas de la matriz. J. M. Jurado EXCEL para Química Analítica 39 3.5.2. Linealidad: Prueba de la falta de ajuste (lack of fit). Se realiza esta prueba para ver si el modelo se ajusta bien a los datos. Como requisito indispensable haber realizado replicados por cada nivel de concentración. En líneas generales se trata de corregir la varianza de los residuales con la varianza de repetición pura promediada para todos los niveles y comprobar si esta varianza corregida (la varianza de falta de ajuste 2
s ) difiere de la de repetición pura mediante un test F. Para poder promediar las varianzas de repetición en cada nivel de concentración es necesario que estas no difieran entre sí. Para ello se realizaría el test de Cochran para detección de outliers dispersivos. Supongamos que realizamos r replicados en n niveles de concentración. Cuando se tiene una serie de resultados agrupados en replicados, y se sospecha que una serie de replicados puede constituir un outlier dispersivo, es decir, que presenta una varianza sospechosa, por ser la mayor de la serie, se puede usar la prueba de Cochran [17]. De esta forma se decide si dicha varianza puede considerarse o no un outlier. El estadístico C de Cochran se calcula según la siguiente fórmula: 2
(61) El valor obtenido se compara con el tabulado C (h, r, α ). Siendo 2
s la varianza sospechosa, h el número de series de replicados, 2
s las varianzas de cada una de estas series, r el número de replicados y α  el nivel de significación. Si el valor calculado es mayor que el teórico, normalmente para un 95% de nivel de confianza, se considera que la serie de replicados sospechosa constituye un outlier dispersivo. En este caso no tiene sentido la prueba de la falta de ajuste, puesto que los datos no son homocedásticos y no se debería aplicar la regresión lineal ordinaria, sino la ponderada. En el ANEXO 2 se presentan los valores tabulados del parámetro C. Una vez hecho esto, si no hay outliers dispersivos, se calcula una varianza de error puro como la media de las varianzas a cada nivel (suponiendo h niveles). 2
(62) J. M. Jurado EXCEL para Química Analítica 40 Con r-1 grados de libertad (se han sumado h varianzas de r-1 grados de libertad pero se ha dividido entre h). A partir del ANOVA se tiene la varianza de los residuales: ( )
(63) Para obtener la varianza de falta de ajuste se resta la de error puro a la de los residuales. 2 2
n s r s
(64) Tras esto se calcula el factor F de falta de ajuste: 2
= (65) Si (0.05, 1, 1)
LOF crit
F F n r r < − − − la varianza de falta de ajuste no será estadísticamente mayor que la de error puro. Esto quiere decir que no existe un error debido a la falta de ajuste del modelo que no pueda ser atribuido al propiamente aleatorio, por lo tanto, el modelo lineal se ajusta perfectamente a los datos experimentales. 3.5.3. Linealidad e intervalo dinámico: método del factor de respuesta de Huber. El método del factor de respuesta de Huber [18] se aplica cuando tenemos rectas de calibrado corregidas con el blanco para que la ordenada sea prácticamente cero. Así, si la respuesta es prácticamente y=bx, los factores de respuesta R=y/x serán próximos a b. Se representa y/x frente a la concentración de analito (o frente al logaritmo si el intervalo de trabajo es amplio, por ejemplo 1-1000 mg/L). Calculamos la mediana de los factores de respuesta, M, trazando una recta en este valor y dos rectas más con valores 1.05M y 0.95M. Los puntos experimentales que se encuentren entre esas dos líneas corresponden al intervalo lineal. J. M. Jurado EXCEL para Química Analítica 41 3.5.4. Comparación con un método de referencia. Para comparar un método que está siendo validado con uno de referencia usando una recta de regresión se representan los resultados obtenidos por el método a validar frente a los obtenidos por el de referencia, obteniendo la ecuación de la recta ( ) ( )
y b s x a s = ± + ± . Donde y son los resultados obtenidos con el método a validar y x los obtenidos por el de referencia. Si los dos métodos llevan a resultados idénticos, es obvio que la pendiente debe ser igual a la unidad y la ordenada igual a cero. Se calculan los parámetros t de Student: 0
= (66) 1
= (67) Se comparan con el valor tabulado t(0.05,n-2). Si b=1 y a= 0, ambos métodos llevan a resultados estadísticamente iguales. Los métodos son equivalentes. Si b = 1 y a ≠ 0, el método a calibrar da un resultado más alto o más bajo que el de referencia, pero en una cantidad fija. es decir, se producen errores sistemáticos constantes o aditivos. J. M. Jurado EXCEL para Química Analítica 42 Si b ≠ 1 y a = 0, se produce un error sistemático proporcional a la concentración de analito en el método a validar. Si b ≠ 0 y a ≠ 0, existen los dos tipos de errores. 3.5.6. Efecto matriz. El efecto matriz se debe, como su propio nombre indica a la matriz de la muestra. Las muestras, salvo en contadas ocasiones, deben tratarse antes de ser medidas. La medida se realiza con una técnica determinada y la señal obtenida se interpola en los valores obtenidos mediante una recta de calibrado. A veces, el hecho de tratar la muestra no es suficiente, pues puede quedar en una matriz que difiere de los estándares utilizados para realizar el calibrado y es muy posible que esa matriz interfiera en el análisis. Se debe estudiar entonces si existe o no dicho efecto. Si es así, no podremos realizar el el análisis utilizando un calibrado externo, sino que deberemos usar el método de adición patrón. Se puede detectar la presencia de efecto matriz comparando el calibrado externo con el de adición patrón. Se trataría de obtener la pendiente de calibrado externo y la de adición patrón y compararlas. En primer lugar realizamos el cociente entre ambas pendientes y calculamos la desviación estándar asociada al mismo: CE
= (68) 2 2
(69) Donde b
son las pendientes de las rectas de calibrado externo y adición patrón, respectivamente, CE
s y AP
s son sus desviaciones estándar correspondientes. Si ambas pendientes son prácticamente iguales, su cociente, R, no debe ser significativamente distinto de 1. Calculando la t de Student como: 1
= (70) El valor calculado se compara con el factor de cobertura k = 2 para un 95 % de nivel de confianza [10]. Si t
es mayor que 2, las dos pendientes no seran comparables y el análisis deberá realizarse por el método de adición patrón. J. M. Jurado EXCEL para Química Analítica 43 Bibliografía 1. The fitness for purpose of analytical methods. A laboratory guide to method validation and related topics. EURACHEM, LGC, Teddington, UK 1998. http://www.eurachem.org/guides/valid.pdf. 2. R. Compañó, A. Rios. Garantía de la calidad en los laboratorios analíticos. Editorial Síntesis, Madrid, 2002. 3. J. M. Jurado, A. Alcázar, F. Pablos, M. J. Martín, A. G. González. Talanta, 2004, 63, 297-302 4. J. Inczedy, T. Lengyel, A. M. Ure. Compendium of Analytical Nomenclature 3rd edition. Blackwell Science, 1998 5. L. Cuadros, A.M. Garcı a and J.M. Bosque, Anal. Lett. 1996, 29, 1231 6. L. Huber. LC-GC Int. 1998, 11, 96–105 7. J. N. Miller, J. C. Miller. Estadística y quimiometría para química analítica. Prentice Hall, Madrid, 2002. 8. J. A. Rafter, M. L. Abell, J. P. Braselton. Multiple comparison methods for means. SIAM Review, 2002, 44, 259-278. 9. H. L. Harter. Tables of range and studentized range. Annals of Mathematical Statistics, 1960, 31, 1122-1147. 10. A. G. González, M. A. Herrador. A practical guide to analytical method validation, including measurement uncertainty and accuracy profiles. Trends in Analytical Chemistry, 2007, 26, 227-238. 11. W. Horwitz. Evaluation of Analytical Methods used for regulations of foods and drugs. Analytical Chemistry, 1982, 54, 67A-76A. 12. R. Wood. How to validate analytical methods. Trends in Analytical Chemistry, 1999, 18, 624-632. 13. AOAC International. Method Validation Program (OMA/PVM Department), including Appendix D: Guidance for collaborative study procedures to validate characteristics of a method of analysis, 2000 J. M. Jurado EXCEL para Química Analítica 44 14. I. Taverniers, M. De Loose, E. van Bockstaele, Trends in Analytical Chemistry, 2004, 23, 535. 15. E. Pitchard, Quality in the analytical chemistry laboratory. ACOL Series, Wiley, Chichester, West Sussex, UK, 1995. 16. International Organization of Standarization (ISO), ISO/DTS 21748, Guide to te use of repeatability, reproducibility and trueness estimates in measurement uncertainty estimation, ISO, Geneva, Switzerland, 2003 17. W. G. Cochran. The distribution of the largest of a set of estimated variances as a fraction of their total. Annals of Eugenics, 1941, 11, 47-61. 18. L. Huber. Validation of analytical methods: review and strategy. LC-GC International. 1998, 11, 96-105. J. M. Jurado EXCEL para Química Analítica 45 ANEXO 1. Tabla de valores estandarizados para el test de Tukey J. M. Jurado EXCEL para Química Analítica 46 ANEXO 2. Parámétro C para el contraste de Cochran. k: conjunto de datos, n: replicados. 95% n.c. J. M. Jurado EXCEL para Química Analítica 47 J. M. Jurado 1.3. Desviación estándar.
EXCEL para Química Analítica
La desviación estándar es un parámetro que da cuenta de la dispersión de los datos obtenidos. En el caso de una población hablaremos de la desviación estándar de la población, σ. Para una muestra se calcula con la siguiente fórmula:
Al cuadrado de la desviación estándar se le denomina varianza (s2). A partir de la desviación estándar se puede también calcular el coeficiente de variación o desviación estándar relativa, que es una medida del error relativo y se usa para comparar resultados medidos en distintas escalas.
CV = RSD = 1.4. Distribución de medidas repetidas.
s ⋅100 (3) x
Cuando realizamos medidas de una población, en ausencia de errores sistemáticos, éstas pueden tomar cualquier valor en torno al valor verdadero o media de la población, dependiendo de la desviación estándar de la misma (del error aleatorio puro). Si se dispusiera de todas las medidas y se representasen las frecuencias con que aparece cada valor, frente a los propios valores, se podría ver la distribución de los datos en torno a la media de la población. En las condiciones mencionadas se obtiene una distribución normal o gaussiana como la de la figura.
Figura 1. Distribución normal
J. M. Jurado La distribución normal viene descrita por la ecuación
1 exp {−( x − µ ) 2 / 2σ 2 } (4) σ 2π
En una distribución normal de media µ y desviación estándar σ, aproximadamente el 68 % de los valores de la población caen dentro del intervalo µ ± 1σ , el 95 % dentro de µ ± 2σ y el 99.7 % dentro de µ ± 3σ. En realidad solo vamos a realizar una serie finita de medidas (replicados) a partir de las cuales se obtiene la media de la muestra y su desviación estándar. La media de la muestra será una estimación de µ y la desviación estándar de la muestra será una estimación de σ. Por lo tanto, la distribución de esta muestra será también normal y se cumplirán las mismas características anteriormente mencionadas. Gráficamente, donde antes obteníamos una curva ahora obtendremos un histograma.
Si dispusiésemos de una serie de medias, cada una obtenida de una serie de medidas individuales, se observaría que estas medias se encuentran más agrupadas en torno a µ que las medidas individuales. La distribución de todas las medias muestrales posibles se denomina distribución muestral de la media o distribución en el muestreo de la media. Su media es la misma que la media de la población original. Su desviación estándar se denomina error estándar de la media (e.e.m.). Existe una relación matemática exacta entre el e.e.m. y la desviación estándar, σ, de la población: e.e.m.= σ/ n . A mayor n, menor será el e.e.m. y menor la dispersión de las medias muestrales en torno a µ. El e.e.m. nos da una medida de la variabilidad de x , además, aunque la distribución de las medidas no sea una distribución normal, la distribución de las medias muestrales tiende a serlo a medida que aumenta n. Esto se conoce como la teoría del Límite Central. Se asume que las medias de muestras muy pequeñas (<5) se distribuyen normalmente.
1.5. Limites de confianza de una medida.
Es poco probable que obtengamos una media de la muestra exactamente igual a la de la población. Generalmente se tiene un intervalo centrado en la media de la muestra que contendrá al valor verdadero con un determinado nivel de confianza. La amplitud de este intervalo depende de dos factores: a) La precisión de los valores individuales, que dependen de la desviación estándar de la población. b) El número de medidas de la muestra, pues el hecho de repetir las medidas implica tener más confianza en el resultado. El límite de confianza de una media viene dado por la expresión:
x ± zs / n (5)
Donde, en ausencia de errores sistemáticos:
x − zs / n < µ < x + zs / n (6)
Donde z=1.96 para un 95% de nivel de confianza, 2.58 para el 99% y 2.97 para el 99.7%.
1.6. Mediana
Para un conjunto de n datos ordenados de menor a mayor, la mediana es el valor de la observación que ocupa el lugar (n+1)/2 si n es impar o bien la media de los valores que ocupan los lugares n/2 y (n/2+1) si n es par.
1.7. Calibrado.
Supongamos que pretendemos realizar la determinación de un determinado analito mediante un método instrumental, de modo que medimos una propiedad física de la muestra que sea directamente proporcional a la concentración de analito. En principio, a partir de la medida de esta propiedad no podremos calcular la concentración de analito. Para poder hacer esto es preciso establecer una relación entre las señales medidas y una serie de patrones de concentración conocida. Esto es lo que se denomina realizar un calibrado. La relación señal-concentración vendrá dada por una función matemática a partir de la cual se puede interpolar el valor de señal obtenido para una muestra y calcular así la concentración de analito en la misma. Generalmente se procura que la relación sea lineal, lo que simplifica mucho los cálculos a realizar.
que consiste en ajustar estos parámetros para minimizar la suma del cuadrado de los residuales:
S = ∑ ( yi − a − bxi ) 2 (7)
Para poder aplicar este método se debe considerar que no existe error en la variable x. … yn. Esto se puede observar en la gráfica de residuales. y2. x2. dicha varianza debe permanecer constante en los distintos niveles (homocedasticidad) y los valores de y deben estar normalmente distribuidos en cada nivel.
Para calcular el valor de estos dos parámetros se utiliza el método de los mínimos cuadrados. tendremos una recta de calibrado de ecuación y = bx + a. …xi. …yi. Además. siendo este considerado sólo en la variable y por mediación de la varianza de repetición pura (obtenida a partir de replicados). M. En el caso más sencillo. … xn que.J. como ya se ha dicho. Jurado
Supongamos que disponemos de una serie de patrones de concentración x1. una vez medidos por una determinada técnica presentan unas señales analíticas y1. donde b es la pendiente y a la ordenada en el origen.
. Si representamos las señales frente a las concentraciones obtenemos la denominada curva de calibración.
valores de r próximos a 1 implican un mejor ajuste.
. r. que estima los errores aleatorios en la dirección y. M.
ˆ ∑( y − y )
i i n i 2
sy / x =
Donde ˆ yi = a + bxi (11) A partir de esta se tiene la desviación estándar de la pendiente y de la ordenada:
sb = sy / x
sa = s y / x
El último parámetro que vamos a considerar es el coeficiente de correlación. que da cuenta de la bondad del ajuste de la curva a los datos experimentales. Su valor oscila entre 0 y 1. Jurado
Loas ecuaciones para calcular la pendiente y ordenada son:
∑ {( x − x )( y − y )}
a = y − bx (9)
Los errores de estos parámetros se calculan a partir de la varianza de regresión.J.
Algunas fórmulas en EXCEL.CORREL() =ESTIMACION. Se explicará mas adelante.8.
. M.DE. Parámetro o función Media Desviación estándar Varianza Raíz cuadrada Número de datos Mediana Pendiente Ordenada en el origen Coeficiente de correlación Estimación lineal Fórmula =PROMEDIO() =DESVEST() =VAR() =RAIZ() =CONTAR() =MEDIANA() =PENDIENTE() =INTERSECCION.LINEAL()
La última es muy útil cuando se quieren obtener los errores asociados a la pendiente y la ordenada en el origen.J. Jurado
n  n 2 2 ( xi − x )   ∑ ( yi − y )  ∑  i  i 
Si tenemos una muestra de señal y0.EJE() =COEF. se puede calcular la concentración correspondiente a la misma como:
x0 = y0 − a (15) b
Y el error de esta concentración viene dado por:
sx0 = sy / x b 1 1+ + n
( y0 − y )
b 2 ∑ ( xi − x )2
Hasta ahora. relacionadas con efectos inespecíficos. es necesario comprobar que dicho método cumple unos determinados requisitos que dependerán de la aplicación que se le quiere dar.
. Exactitud. Estas son las designadas como efecto matriz y.
Un método analítico es un procedimiento que obtiene unos resultados a partir del análisis químico de muestras. Al máximo de selectividad se le denomina especificidad. A veces. Selectividad y especificidad.J. es decir. Validación de métodos analíticos. Pero existen otro tipo de interferencias.1.2. Estos compuestos pasan a denominarse interferencias. así como la identificar las influencias que pueden cambiar dichas características [1]. 2. Un estudio de selectividad implica investigar la capacidad del método para medir al analito en presencia de interferencias añadidas a la muestra [1]. En todo método analítico hay un paso fundamental en el que se produce la medida de una señal relacionada con la concentración (análisis cuantitativo) o la presencia (análisis cualitativo) de un determinado analito. un método será específico cuando la señal medida solo proviene del analito de interés [2]. solo hemos considerado las interferencias que pueden producir un error sistemático constante. hablaremos de el más adelante. Este proceso de verificación se conoce como validación y debe establecer las características de funcionamiento y las limitaciones del método. M. 2. sustancias que contribuyen a la señal medida con un valor determinado. aunque no existe ninguna aplicación en Excel específica para cuantificar este efecto. De aquí surge en concepto de selectividad. entendiéndose esta como el parámetro que define de forma cualitativa la extensión en que otras sustancias interfieren en la determinación de un analito usando un método determinado. la presencia de otros compuestos. pueden contribuir a dicha señal. Es por tanto necesario asegurarse de donde procede la señal que se está midiendo. que pueden producir errores proporcionales. es decir. no siendo esto objeto del presente texto. A continuación se detallan los distintos parámetros relacionados con la calidad de un método analítico que deben ser objeto de estudio en el proceso de validación. Antes de ser usado. pudiendo causar un error sistemático en la determinación del analito. química o físicamente similares al analito. Jurado
el primer término tiende a la varianza de la medida (dividiendo entre n-1 en vez de n) y podría expresarse:
E 2 = s 2 + ( x − µ )2 (20)
Con lo cual vemos como el error cuadrático medio correspondiente a la exactitud se equipara a la suma de la varianza de repetición pura (error aleatorio) más el cuadrado de la veracidad (error sistemático). obtenemos el error cuadrático medio:
( xi − x ) 2 n + ( x − µ ) 2 (19)
Si n es suficientemente grande. Generalmente se expresa en términos de error o bias. según la guía para validación de métodos de la EURACHEM [1].
La exactitud da cuenta de la proximidad entre el resultado obtenido por un método y el valor “real”. La exactitud es máxima cuando el sesgo es nulo y la dispersión aleatoria es mínima. La precisión da cuenta de la proximidad entre sí de esta serie de resultados. se divide en dos componentes. definido como la diferencia entre el resultado de medida y el valor real [2]. M. hace falta un parámetro que permita medir la exactitud cuantitativamente. la veracidad y la precisión. Este parámetro es el error que. Como se trata de un parámetro cualitativo (un método será muy o poco exacto) para describir el error asociado a un resultado.J. La veracidad es un parámetro que expresa la proximidad de la media de una serie de resultados obtenidos con el método al valor real. el error cuadrático se calcularía como:
∑ Ei2 = ∑ ( xi − x) 2 + n( x − µ ) 2 (18)
Si dividimos por n. Si expresamos la exactitud en función del error Ei. para una determinada medida i: Ei = xi − µ = ( xi − x ) + ( x − µ ) (17) En el caso de tener un replicado consistente en n medidas.
Como última instancia. 2. pudiendo variar condiciones tales como el equipo.3.
Como se ha dicho anteriormente. La precisión se divide en dos parámetros repetitividad y reproducibilidad. se pueden realizar dichos ensayos de recuperación añadiendo al analito sobre la propia muestra [3]. 10
. preparando los patrones de calibración cada vez.J.4. Como se explicará más adelante. Se habla de reproducibilidad interlaboratorio cuando las medidas se realizan en laboratorios distintos (ensayos de intercomparación). este grado de concordancia se denomina precisión. para compara dos valores será necesario conocer la precisión con la que este ha sido obtenido. Reproducibilidad. Como se dijo anteriormente. M.1. Esto se realiza. reactivos o analistas. 2. Para asegurar la veracidad es necesario compara el valor medio medido con un valor de referencia (que será el valor real o el considerado como tal).4. En caso de no existir materiales de referencia o no tener posibilidad de aplicar el método de referencia. en un periodo de tiempo corto. generalmente comparando los resultados obtenidos por el método con el valor de un material de referencia certificado o comparando con un método ya validado (que suele haber sido sometido a un ejercicio de intercomparación). 2. por ello la exactitud se evalúa a partir de la veracidad en conjunto con la precisión.4. Se define la repetitividad como la desviación estándar obtenida al analizar una misma muestra varias veces. la veracidad da cuenta de lo cercano del resultado medio obtenido con el método y el valor real. Precisión. Repetitividad. se puede simular la matriz de la muestra añadiendo al analito en una cantidad conocida y efectuar un ensayo de recuperación. A esta reproducibilidad intralaboratorio se la conoce como precisión intermedia. Cuando realizamos una medida se deben obtener replicados. reactivos o analista. La reproducibilidad se define como la desviación estándar obtenida al analizar varias veces la muestra en días distintos.2. a partir de los cuales se puede observar el grado de concordancia que estos resultados tienen entre sí. Jurado 2. en distintos días y cambiando en la medida de lo posible de equipo y analista. La reproducibilidad intralaboratorio se realiza en un solo laboratorio. sin cambiar de equipo de medida. Veracidad.
Jurado 2. Si analizamos un blanco de la muestra varias veces y obtuviéramos la señal media ( Y0 ) y su desviación estándar.5. Otro método se basa en que la desviación estándar
. Límites de detección y cuantificación.
Según la IUPAC la sensibilidad se define como el cociente entre la señal medida y la concentración de analito. la sensibilidad es la pendiente de la curva a una concentración dada. Si tenemos una recta de calibrado de ecuación y = bx + a. sY0 .6. En caso de ser una recta. Sensibilidad. que lleva a una concentración. De manera análoga al LOD. las señales mayores a YLOD = Y0 ± 3sY0 caerían fuera de la gausiana del blanco analítico y podríamos asignarlas a la presencia de analito. coincide con la pendiente de la misma. el LOD se calculará como:
X LOD = 3s a (22) b
siendo sa la desviación estándar de la ordenada. El límite de detección (LOD) se define como la cantidad de analito presente que se deriva de la menor señal analítica que puede detectarse con certeza razonable (YLOD) [4]. la IUPAC propone como LOQ al valor YLOQ = Y0 ± 10 sY0 . El límite de cuantificación (LOQ) se define como la cantidad de analito presente derivada de la menor señal analítica que puede ser medida con precisión y exactitud.7.999. 2. Estos parámetros están relacionados con la sensibilidad del método.J. considerando lineal el calibrado si r > 0. Cuando trabajamos con curvas de calibrado. Linealidad e intervalo dinámico
10 sa (23) b
La linealidad a veces se examina a partir del coeficiente de correlación. M. el LOD se calcularía como:
X LOD = YLOD − a (21) b
Como la señal del blanco debe coincidir con la ordenada del blanco analítico (puesto que es la señal que se obtendría a una concentración de analito igual a cero). a partir de la recta de calibrado:
X LOQ = 2.
no pueden ser inferiores al límite de cuantificación. El porcentaje de linealidad se calcula como [5]:
 s  % Lin = 100  1 − b  (24) b 
Estos parámetros no son del todo fidedignos. se puede concluir que variables son más significativas a la hora de realizar las medidas y por tanto. que se consideran las variables que gobiernan el experimento de medida. y se comprueba su influencia en la señal analítica.
2. siguiendo un determinado diseño experimental.
2. deben ser mejor controladas. Jurado
relativa de la pendiente no supere el 5%. Ámbito de validación de métodos analíticos
Como ya hemos dicho.8. Para el estudio de robustez. Ambos procedimientos los explicaremos detenidamente cuando estudiemos la aplicación de EXCEL. M. con lo cual.9. Robustez
Es la inercia que presenta un método analítico a modificar su señal cuando tienen lugar pequeños cambios en las condiciones ambientales u operativas. los valores deben presentar un nivel aceptable de precisión y exactitud. Esta es la conocida como linealidad online. El intervalo dinámico o rango lineal es el rango de concentraciones donde existe una relación lineal entre estas y las señales analíticas.J. En este intervalo. Parámetro Especificidad Precisión Veracidad Linealidad Rango LOD/LOQ Robustez LOD X X Análisis Cualitativo X Análisis Cuantitativo X X X X Análisis de trazas X X X X X LOQ X
. el tipo de validación que debemos realizar va a depender de los requisitos de nuestro método. estas condiciones se modifican ligeramente. Con ello. siendo lo correcto realizar un análisis estadístico de los residuales mediante el análisis de la varianza (ANOVA) o la aproximación práctica de Huber [6].
EXCEL como herramienta en el proceso de validación.
Si no aparece debemos activar el paquete en C:\Archivos de programa\Microsoft Office\OFFICE11\Macros\Análisis\ANALYS32.
Aunque todo esto se puede hacer de forma manual. es decir. EXCEL incluye un paquete de herramientas para realizar este tipo de análisis de manera sencilla. mediana. Para ello hay que abrir el menú Herramientas>Complementos y ver si aparece la opción Herramientas para análisis. El primer paso será comprobar que está instalado o activado.
Las utilidades estadísticas de EXCEL que emplearemos en validación serán: Prueba F de Snedecor-Fisher Prueba t de Student ANOVA de un factor Estadística descriptiva (media. Jurado
3. Herramientas para análisis. varianza…) Regresión
3. En este momento dispondremos de la herramienta análisis de datos.J. introduciendo formulas en la hoja de cálculo. M.1.
Comprobar que en la reproducibilidad no interviene factores distintos a los encontrados en el estudio de repetitividad. F = 2 s2 ≥ 1 .2. Aplicaciones .INV(p.05. Resolución con EXCEL Supongamos que hemos realizado una serie de medidas de una muestra con dos métodos analíticos y queremos comprobar si ambos son igual de precisos.1. 3.05. es decir 95% de nivel de confianza) 3.
.Fisher.F. Lo primero es tener los datos en EXCEL. .n2-1. . Sean dos conjuntos de datos: {x1} con varianza s12 y tamaño n1
2 {x2} con varianza s2 y tamaño n2
Consideramos la hipótesis nula que dice que las varianzas no presentan diferencias
2 2 significativas entre sí.2. n1-1). Jurado
3. Si F<Ftab admitimos H0.n1-1). s12
n2-1.J.2. M.2.Comprobar que las varianzas son iguales para ver que tipo de test de comparación de medias usar posteriormente. Se compara con Ftab(0.Comprobar que dos métodos son igual de precisos. Prueba F de Snedecor. Los valores de F para una cola se pueden obtener usando la función =DISTR. siendo p la probabilidad (normalmente 0. Si s2 ≥ s12 . H0: σ 12 = σ 2 .
En una prueba estadística empleada para la comparación de varianzas.
M. El resultado será:
. Otra opción es hacerlo en una hoja nueva. Si marcamos Rótulos. Alfa es el indica la probabilidad o nivel de confianza que usaremos (0.J. En rango de salida se indica la primera celda del rango donde se colocarán los resultados.05 es un 95%). el programa considera que en la primera celda del rango seleccionado se encuentra el nombre del conjunto de datos (deberíamos haber seleccionado A1:A7 y B1:B9). Jurado
Del menú Herramientas>Análisis de datos usamos la herramienta “Prueba F para varianza de dos muestras”
Pulsamos en la flecha roja del cuadro rango para la variable 1 y seleccionamos el rango A2:A7. Es aconsejable que sea en un extremo para que no sobrescriba nada. En rango para la variable 2 seleccionamos el rango B2:B9. Supongamos que hemos elegido la celda D1 como rango de salida.
05. en este caso. es decir.05.1.
2 Primero se comprueba que s12 = s2 mediante la prueba F de Fisher.J. n1 + n2 .3. Así tenemos dos
2 a) Si s12 = s2 se calcula una varianza promedio
s Se calcula el estadístico t:
( n1 − 1) s12 + ( n2 − 1) s22 =
x1 − x2 1 1 s + n1 n2
.INV(0.3. Comparación de valores medios de dos conjuntos de datos. Si necesitásemos el valor de F crítico para comparar las varianzas a manos se puede obtener de tablas o bien con la formula =DISTR. ambos métodos son igual de precisos. Dados dos conjuntos de datos: Conjunto 1 2 Media Varianza
Para decidir si las medias muestrales x1 y x2 difieren significativamente se toma la hipótesis nula H0: µ1 = µ2 . Si t < tcrit se acepta la hipótesis nula.
3. Este valor de t se compara con el valor tabulado t(0.2). Prueba t de Student. M. Jurado
Se comprueba que F < Fcrit. las dos medias son estadísticamente iguales.5.7).F. que tendrá n1 + n2 – 2 grados de libertad.
3. 3. en el que mediante el test de Fisher se dedujo que las varianzas eran similares. M.
En Herramientas>Análisis de datos se escoge la opción “Prueba t para dos muestras suponiendo varianzas iguales”
.3.1.T.05.2.1. Resolución con EXCEL.g. a) Conjunto de datos con varianzas estadísticamente similares. Comparación de los resultados obtenidos con un método en desarrollo y uno de referencia.ν eff ).J.). Comparación del valor medio obtenido con para un material de referencia aplicando un determinado método y el valor de un certificado.1. Jurado
2 b) Si s12 ≠ s2 se calcula el estadístico t como:
Ahora los grados de libertad se deben calcular como:
2  s12 s2  +    n1 n2  2
  s14 s4 + 2 2  2   n1 ( n1 − 1) n2 ( n2 − 1) 
El valor de t se compara con t(0. Supongamos que partimos del ejemplo anterior.3.INV(p. para una cola se toma probabilidad 2*p (siendo p el nivel de confianza deseado). Aplicaciones. Los valores de t se pueden obtener con la función =DISTR.l. Estos valores son adecuados para contraste de dos colas.
M. Nótese que es posible tomar distintas hipótesis nulas en la casilla Diferencia hipotética entre las medias. Jurado
La forma de escoger los rangos.J. b) conjunto de datos con varianzas distintas. Considérese el conjunto de datos:
. El resultado es:
Como en valor absoluto el t calculado es menor que el t crítico (contraste de dos colas). El hecho de usar un contraste de dos colas implica que no nos interesa el sentido en que difieren las medias (es decir si la diferencia es positiva o negativa). la hipótesis nula se toma como cierta y las dos medias muestrales no difieren significativamente. Si pretendemos ver solo si un valor aumenta en relación al otro se usa un contraste de una cola. nivel de confianza y rango de salida son similares a la prueba F.
3.... ≠ xi ≠ .. se debe elegir “Prueba t para dos muestras suponiendo varianzas desiguales”..x2i . ≠ xn . Los rangos....... Prueba t para valores emparejados..... pero se ajusta al entero más próximo).x1i . A la hora de realizar el contrate t.2. por lo que las cantidades x1 ≠ .x2 n con el
Comparando con el valor de t crítico para dos colas. Se determinan los valores mediante dos métodos y se pretende comprobar que no existen diferencias significativas entre los resultados obtenidos por ambos.J. Nótese que se trata de n muestras distintas.. las medias no difieren significativamente...x1n con el método 1 y x21 .2868. Jurado
Al aplicar el test F:
Las varianzas difieren significativamente.xn . Al aplicar cada método se obtienen los valores x11 . M. probabilidad e hipótesis nula se introducen de igual modo que en el ejemplo anterior....3. Nótese que los grados de libertad efectivos son 5 (en realidad 5..xi . Supongamos que tenemos un conjunto de muestras de valores x1 .
Se aplica en el caso de que las muestras analizadas por los dos métodos no sean grandes y sea necesario utilizar distintas muestras en el estudio. M.2. 3. el valor medio de todas las diferencias y su desviación estándar:
∑(d
Si los métodos son equivalentes y levan al mismo resultado las diferencias di deben ser nulas.3. Se compara este valor con t(0.
Se utiliza la “Prueba t para dos muestras emparejadas” del menú Herramientas>Análisis de datos.2.2. 3.3. Es decir.J. Se calcula el estadístico t:
t = d n / sd (31)
El número de grados de libertad de t es n-1. Esto es lo mismo que considerar la hipótesis nula H0: µd = 0. Aplicaciones.1.n-1). Se calculan las diferencias d i = x1i − x2i . Jurado
método 2. El formulario de entrada de datos es similar a los de los ejemplos anteriores. Resolución con EXCEL Consideramos el ejemplo de la determinación de paracetamol en diez pastillas de distintos lotes mediante un método espectrofotométrico y otro de espectroscopía IR [7]. se analizan distintas muestras por ambos métodos y no una misma muestra replicada por ambos métodos. El resultado es:
26) los métodos no proporcionan resultados significativamente diferentes.
El análisis de la varianza (ANOVA) es un método que se emplea para comparar resultados obtenidos por distintos métodos.88) es menor que el crítico (2. laboratorios. laboratorio. la primera. En el caso de realiza medidas replicadas variando uno de los factores anteriormente mencionados siempre existen dos fuentes de error. Si tuviésemos un caso en que el análisis se realiza por distintos laboratorios utilizando distintos métodos. Jurado
Como el valor experimental (0. con lo que se habla de ANOVA de un factor. analistas.J. También permite separar las contribuciones de uno o más factores ala varianza global del sistema.
La media obtenida para cada método será:
La media global obtenida por todos los métodos se calculará como:
.4. deberíamos controlar dos factores y se utilizaría un ANOVA de dos factores. cuando el número de medias obtenidas es superior a dos. analista. etc. ANOVA de un factor. etc.
3. Consideremos que realizamos un mismo análisis por i métodos distintos (i varía de 1 a
h) y realizamos j replicados (j varia de 1 a n) con cada método. el error aleatorio de medida. Mediante el ANOVA se puede controlar el error introducido por esta segunda fuente. M. la segunda los errores debidos al cambio de método.
. Si todas las medidas se extraen de una población de varianza σ 2 . xhn
. .. .. Método h
x11 x21 x12 x22
. Método i . . también conocida como varianza
dentro de las muestras (o del grupo). .J. sPE . xh
Para cada método se tiene una varianza calculada como:
∑ ( xij − xi )
En el ANOVA la hipótesis nula es que todas las muestras se extraen de una población de media µ y varianza σ 2 . .. conocida como varianza entre muestras (o entre grupos). . es decir. . . . . . . .4). .
x1n x2 n
. xi . se calcula como promedio de las varianzas individuales si2 :
2 sPE =
ij − xi )
h(n − 1)
La varianza debida al factor (en nuestro caso a los distintos métodos usados). xhj
. . entonces sus medias (obtenidas variando el método empleado.. . . . . ..
2 La varianza debida al error puramente aleatorio. xh 2
EXCEL para Química Analítica Media Método 1 Método 2 ..
x1 j x2 j
.... . . xij . . σ 2 . M. . es una estimación de la varianza de la población. . ... la varianza debida al factor se calcula como:
. . En general. . estudiando la variación dentro de las muestras y la variación entre muestras. .. .. . . . .. variando un factor controlado) proceden de una población de varianza σ 2 / n (Véase sección 1. .. xi 2 . .. xi1 . xh1
. . . xin . . Esta varianza se puede estimar de dos formas. . .
que F sea mayor que el valor de F crítico. Se comprueba esto mediante un test de Fisher.05. Un resultado significativo en el ANOVA puede deberse a varias cosas: que una media difiera mucho a las demás. que todas las medias difieran entre sí o que existan varios grupos de medias distintas. Ésta se calcula como:
MDS = sPE
2 ⋅ th ( n −1) (39) n
.J. es decir. que la variación debida al factor no debe se significativamente distinta a la varianza de error puro. Pero puede que no lo sea. Hasta ahora hemos supuesto que la hipótesis nula es verdadera. M. las dos estimaciones no deberían diferir significativamente. Si el valor de F calculado es menor que el tabulado. es decir. es decir. Jurado
n ∑ ( xi − x )
2 sFactor =
Si la hipótesis nula es correcta. no solo si difieren). calculando F como:
2 sFactor (37) 2 sPE
Se compara F con el valor tabulado F(0. Si la hipótesis nula es verdadera existe una tercera forma de estimar la varianza que consiste en tratar los datos como si fueran una muestra grande de varianza:
2 sTotal = i j
Pero esta fórmula no se usa en el contraste. Para poder ver lo que ocurre se colocan las medias ordenadas de forma creciente y se comparan las diferencias de los valores adyacentes con una cantidad denominada la mínima diferencia significativa. Esto es lo mismo que decir que el hecho de variar el factor (el método de medida en nuestro ejemplo) no introduce un error significativo en comparación con el puramente aleatorio. h-1. no existen diferencias significativas entre las dos varianzas. etc. h(n-1)) para contraste de una cola (puesto que nos interesa saber si la varianza del factor es mayor que la de error puro. el hecho de variar el factor (cambiar de método) no introduce un error significativo.
Si una vez realizados los cálculos de ANOVA se han encontrado diferencias significativas entre la varianza de un factor y la del error puro. que µ p = µ q.ν
 1 1 2 sPE  +  n   p nq 
Donde. si x p es la media muestral de una población de media µ p y xq es la
media muestral de una población de media µ q. De este modo. Jurado
Otra posibilidad es comparar las medias dos a dos mediante el test de de la genuina diferencia significativa de Tukey. r. Teniendo en cuenta que. debemos averiguar que nivel de ese factor es el responsable de los errores sistemáticos que producen dichas diferencias. Para que no haya diferencia significativa entre los valores promedio de estos dos niveles.ν es el valor tabulado de la distribución del rango estandarizado [9]. Se define el intervalo de confianza de Tukey como:
µ p − µq ∈ ( x p − xq ) ±
α qr . α . Supongamos que queremos comparar dos niveles p y q del factor (método en nuestro caso).J. se debe cumplir. el intervalo de confianza de Tukey queda:
. Además. por el que se comparan dos a dos las medias obtenidas a cada nivel del factor. de distribución normal. de media µ y varianza
σ 2. se sustituye ésta por la varianza del error puro. x p y xq . ν . qα. si en nuestro caso np = nq = n. Con este fin se realiza el ensayo Tukey-HSD (Honest Significant Difference) [8]. como hipótesis nula. que r depende del numero de repeticiones. M. se debe demostrar que ambos se corresponden con las medias muestrales de una misma población. y de los grados de libertad utilizados en la desviación estándar que se usa en el calculo de este intervalo. Aquel nivel que presente un valor promedio significativamente diferente del resto será el responsable del sesgo detectado y se elimina para recalcular el ANOVA. la varianza utilizada en la expresión anterior es una estimación de la de la población y debe considerar todas las muestras. así como del nivel de confianza.
El análisis de la varianza se aplica. También se utiliza para calcular la precisión intermedia de un método cuando se analizan muestras de referencia en distintos días. El nivel que presente diferencias significativas con el resto.ν
EXCEL para Química Analítica 1 (41) n
Como µ p − µ q = 0 . etc. analistas. Una tabla de valores de la distribución del rango estandarizado se incluye en el ANEXO 1. el intervalo de confianza de Tukey debe contener el cero.ν
( x p − xq ) −
q 1 1 < 0 < ( x p − xq ) + r .
. Este ejemplo lo explicaremos más adelante. Aplicaciones. y se calcula un valor de q experimental según:
qexp = x p − xq sPE
Se compara este valor experimental de q con el valor tabulado [9].J. 3.ν sPE n n 2
En la práctica. puesto que EXCEL realiza un análisis de la varianza de los residuales. Es muy útil para detectar resultados anómalos en los ensayos de intercomparación empleados para certificar materiales de referencia o desarrollar métodos de referencia. Si qexp>qtab se rechazaría la hipótesis nula. en la comparación de varios métodos. se rechaza la hipótesis nula. se considera como causa del efecto significativo en la varianza total. laboratorios. µ p = µ q. con lo que las dos medias presentarían diferencias significativas entre sí. Jurado
α qr . se parte de la hipótesis nula. Esto lo veremos en el apartado correspondiente a la estimación lineal. como se dijo anteriormente. Ensayos de recuperación.
α qr . Otra aplicación tiene que ver con la regresión lineal. si no es así. M.1.4.
J.4. señalando la casilla de rótulos para que tome de la primera fila el nombre de los grupos. El resultado final es el siguiente:
La columna de promedio de cuadrados indica el valor de varianza de error puro (dentro de los grupos) y la debida al cambio en el factor (entre grupos). El formulario de entradas de datos es el siguiente:
Como el rango de entrada seleccionamos A1:D4. C) con luz tenue y D) luz brillante. M.66) es mayor que el crítico para un nivel de confianza del 95 % y 3 y 8 grados de
. Como se observa el valor de F calculado (20. Hay que cerciorarse que tenemos marcado que los datos se agrupan por columnas. B) almacenada en oscuridad. que en este caso es el cambio en condiciones de almacenamiento.2. Resolución con EXCEL.
Supongamos que realizamos medidas de fluorescencia a muestras de quinina almacenadas de cuatro formas distintas (ejemplo adaptado de [7]): A) disoluciones recién preparadas. Jurado 3. Los datos los ordenamos como se observa a continuación:
En el menú Herramientas>Análisis de datos seleccionamos “Análisis de varianza de un factor”. El rango de salida funciona como anteriormente se ha descrito.
Supongamos que realizamos un estudio de reproducibilidad a un solo nivel de concentración. minimalista). La varianza neta “entre grupos se calcula como:
2 sB = 2 2 sFactor − sPE (44) n
La reproducibilidad se calcula finalmente como la suma de la repetitividad (error puro) y la varianza neta debida al factor (neta entre grupos):
2 2 2 sR = sPE + sB (45)
Para comparar la tolerancia correspondiente a esta precisión intermedia se puede usar la conocida como Trompeta de Horwitz [11]. realizando el estudio en distintos niveles de concentración. La varianza de repetitividad o de error puro vendrá dada por la ecuación 35. al prevalecer el factor de las condiciones de trabajo sobre los niveles de concentración deberíamos realizar un ANOVA de dos factores anidados. En este caso se va ha emplear para separar y estimar distintas fuentes de variación. Consideremos que realizamos n replicados durante h dias (reproducibilidad entre dias. Puesto que en la ecuación 36 están incluidas las medias de cada grupo. Lo ideal sería añadir otro factor. y el debido al cambio de condiciones (día.) que implica el hecho de estudiar la reproducibilidad. es decir. Jurado
libertad (4..3. Existen entonces diferencias significativas entre las medias obtenidas tras las distintas formas de almacenamiento de las muestras. La precisión a partir del análisis de la varianza.4. se debe calcular una varianza neta “entre grupos” que dependa sólo de la varianza entre las medias de los grupos y la media global. Por ello.
. pero entonces realizaríamos un ANOVA de dos factores. Hasta ahora se ha aplicado el ANOVA para la comparación de múltiples medias. la debida al error puramente aleatorio. se estaría considerando el efecto de la repetitividad dentro del efecto del factor. 3.. mientras que la debida al factor vendrá dada por la ecuación 36. el error de replicación.J. Los test de minima diferencia significativa o de Tukey deben realizarse manualmente (EXCEL no implementa los mismos). laboratorio. analista. reactivos. xi . M. Es más. Tendremos dos fuentes de variación.066).
. se puede también comparar con los valores establecidos por la AOAC [13] (Asociación Oficial de Químicos Analíticos). Además de comparar con los valores de Horwitz.5log c ) (46)
El valor de RSD obtenido se compara con el predicho por Horwitz (RSDH) mediante el parámetro Horrat [12].J. Jurado
Horwitz dedujo una expresión para predecir el valor esperado para la desviación estándar relativa para la precisión intermedia (o interlaboratorio) a partir de la concentración de analito. c (en tanto por 1).
RSDH = 2(1−0.
Horrat = RSDR (47) RSDH
Si el valor del parámetro Horrat es igual o menor a 2 se puede decir que el método tiene valores aceptables de precisión intermedia. Estos valores y los obtenidos por Horwitz se presentan en la siguiente tabla.
1% 100 ppm 10 ppm 1 ppm 100 ppb 10 ppb 1 ppb Horwitz %RSD AOAC PVM %RSD 2 2.6 * RSDH Reproducibilidad: RSDR = 0.J.0000001 1 10
Unit 100% 10% 1% 0.3 11 15 21 30
10−2 10−3 10−4 10−5 10 10
Algunos requisitos prácticos relacionados con los estudios interlaboratorio son los siguientes [14]: Repetitividad: RSDr = 0. Precisión y veracidad a partir del ANOVA.2-1 * RSDH
3.7 3.0001 0.3 7.8 4 5.2-0.001 0.4. El sesgo se puede entonces calcular como:
1 h n ∑∑ xij − T = x − T (49) hn i j
La incertidumbre de este sesgo se puede calcular como [16]:
. A partir de un mismo diseño de ANOVA es posible estudiar tanto la precisión (veáse el apartado anterior) como la veracidad.7 8 11.3 1.000001 0.5-0.8 2.00001 0.7 5. Simplemente hay que considerar que cada medida xij tendrá un sesgo
δ ij = xij − T (48)
Siendo T el valor certificado de un CRM (pues el valor verdadero µ no se conoce). M.1 0.4.01 0.3 16 22.3 * RSDH Reproducibilidad (precisión intermedia): RSDR = 0.6 32 45.3 1.5-2 * RSDH
Para validaciones intralaboratorio [15]: Repetitividad: RSDr = 0. Jurado
Analyte (%) Analyte fraction 100 10 1 0.
u (δ )
Al desarrollar un método totalmente nuevo en una matriz no muy estudiada puede ocurrir que no existan materiales de referencia ni otros métodos validados para comparar el método en desarrollo. A veces. cuando no es posible preparar una matriz similar a la muestra se realizan los ensayos sobre la propia muestra. en “Agregar línea de tendencia”. Se calcula a partir de la ecuación
2 43. Regresión Ya se ha explicado la forma de calcular los coeficientes de regresión. M. Estos consisten en añadir una cantidad de analito conocida a una matriz que simule a la muestra. como se indica en la figura. se puede obtener la recta de regresión con su ecuación correspondiente. Ensayos de recuperación. Se deben ensayar al menos tres niveles de concentración.5.4. analizarla y ver la cantidad de analito que se recupera. Se calcula la recuperación como: R= Cobs − Cnative (51) Cspike
Donde Cobs es la concentración medida con nuestro método. otro medio y otro alto. pulsando en el botón secundario del ratón. uno bajo. Jurado
 n −1  2 sR − sr2    n  (50) u 2 (δ ) = h
2 Donde sR es la varianza de reproducibilidad (puesto que consideramos un ensayo 2 interlaboratorio).J. Si realizamos j replicados en cada uno de los i niveles podemos ver si las recuperaciones son equivalentes en cada uno de los niveles estudiados. Entonces se hace necesario el empleo de ensayos de recuperación [3].5. Por todos es conocido que si representamos los datos en Excel y seleccionamos los puntos. 3. que puede ser sustituida por la varianza de precisión intermedia sIP en
caso de estar realizando la validación intralaboratorio. Se considera que hay ausencia de sesgo para 3. Los pasos a seguir son los siguientes:
. cuando la relación entre las variables x e y es lineal.
δ < 2. La varianza de repetitividad sr2 no es otra que la varianza de error puro sPE obtenida
en el ANOVA. Cnative es la concentración que tiene la muestra (en el caso de ser una matriz preparada debe ser cero) y Cspike es la cantidad que hemos añadido.
Marcar “Presentar ecuación en el gráfico” y “Presentar valor de R cuadrado en el gráfico”. M.
.J. Jurado
Se selecciona el tipo de estimación y se pulsa en la pestaña “Opciones”. Pulsar en aceptar.
Partimos de los datos:
El formulario de la herramienta Regresión es el siguiente:
. M. Jurado
Ahora vamos a explicar la herramienta “Regresión” del menú Herramientas>Análisis de datos.J.
Este último presenta los valores de y ordenados de menor a mayor indicando el percentil al que pertenece.J. Jurado
Loas casillas de rangos de entrada y salida y la de rótulos ya son conocidas por nosotros. Residuos. Resmean es la media de los i residuales y Ressd su desviación estándar. puede decirse que parece indicar que los datos provienen de una distribución normal. Finalmente se puede obtener la curva de probabilidad normal. Gráfico de residuales y Curva de regresión ajustada. por ello se recomienda utilizar métodos objetivos. Seleccionamos Rótulos. Nosotros no la forzaremos. pero téngase en cuenta que en algunos casos. M. residuales estandarizados.i = Re si − Re smean (52) Re ssd
Donde Resst. En el cuadro “Residuales” se puede especificar si se quieren tener los residuales. si se desea otro nivel de confianza se puede activar la casilla “Nivel de confianza.i es el i-ésimo residual estandarizado. Los residuos estándares se obtienen como:
Re sst . Esta herramienta ofrece siempre el intervalo de confianza de la pendiente y la ordenada en el origen al 95% de nivel de confianza. La casilla “Constante igual a cero” se utiliza si queremos forzar a la curva a pasar por el cero. aunque los puntos se ajusten a una línea recta puede que los datos no provengan de una distribución normal. por ejemplo para el 90%. Seleccionamos B1:B9 para el rango y y A1:A9 para el x. El resultado será el siguiente:
. Resi el i-ésimo residual. el gráfico de residuales y el de la propia curva ajustada. Si los puntos parecen ajustarse a una línea recta.
se obtiene a partir de la suma de cuadrados de los residuales (SCRes) y la suma de cuadrados total (SCT):
ˆ SC Re s = ∑ ( yi − yi ) 2 (53)
SCT = ∑ ( yi − y )2 (54)
SC Re s (55) SCT
El ajustado se calcula a partir de las varianzas de los residuales y total. también llamados cuadrados medios o promedio. El coeficiente de determinación R2. que coincide con el cuadrado del coeficiente de correlación. Jurado
En las estadísticas de regresión se tiene el coeficiente de correlación y el error típico (sy/x).J. CMRes y CMT.
. M. No es el criterio que tomaremos. trataremos este tema cuando hablemos de linealidad y rango lineal.
CM Re s =
CMT = n −1
R 2′ = 1 −
CM Re s (58) CMT
Cuanto más próximo a 1 sea el coeficiente de correlación mejor será el ajuste lineal.
porque es poco probable que los datos se organicen de este modo por casualidad. Si Fexp>F(0. si el valor de F es muy elevado. pero en vez de comparar. r será próximo a 1. El hecho de que EXCEL trabaje con un valor de probabilidad es similar a comparar con el valor tabulado de F para 1 y n-2 grados de libertad. Excel solo calcula las varianzas de regresión y la de los residuales (cuarta columna). ˆ ∑( y − y)
(n − 1) − (n − 2) ˆ ∑(y − y )
n−2 En general. calcula directamente la probabilidad de que el F calculado sea mayor que el tabulado. es decir.J. Como ya hemos dicho. es poco probable que la suma de cuadrados totales sea mucho mayor que la de residuales y r sea próximo a uno. A partir de ellas calcula un Valor de F (quinta columna) y calcula una probabilidad para este valor (es lo que llama valor crítico de F). la suma de cuadrados de residuales será mucho menor que la de regresión. la suma de cuadrados de los residuales y la de regresión:
ˆ ˆ ∑ ( y − y) = ∑ (y − y) + ∑( y − y )
Estas sumas de cuadrados aparecen en la tercera columna y dividiendo entre los grados de libertad de la segunda se obtienen las correspondientes varianzas (llamadas por EXCEL promedio de cuadrados). Jurado
El análisis de la varianza aquí realizado se basa en que la suma de cuadrados totales se descompone en dos componentes. Si la probabilidad es muy pequeña es poco probable que los datos lleven a un valor tan elevado de F. n-2) la correlación lineal es significativa y no puede atribuirse a la casualidad. M. Entonces. con lo que también será mucho menor que la suma de cuadrados totales. 1.05. según la expresión 55. por ejemplo) teniendo en cuenta que los grados de libertad para la varianza de regresión es (n-1-n-2=1) y la de la varianza de los residuales n-2. Los parámetros ajustados de pendiente y ordenada se obtienen en:
. Esto quiere decir que el ajuste lineal no es fruto de la casualidad. nosotros estamos más acostumbrados a comparar el valor de F experimental con un valor tabulado de F para un nivel de confianza determinado (95 %.
05. En la columna de error típico aparecen las desviaciones estándar de la ordenada y la pendiente. M.
Observando los residuales podemos ver si: a) Si el modelo se ajusta bien a los datos. 36
. También se observan los valores de los residuales y los gráficos de residuales y de la curva ajustada. comparando con t(0. c) Si el modelo es inadecuado. Igual que en el caso anterior se tiene la probabilidad de que éste t experimental sea mayor que el teórico. n-2) grados de libertad. Se observa cuando los residuales de valores próximos en x son todos del mismo signo. cambiando de signo en otro rango de x.n-2).05. b) Si los datos son heterocedásticos. Los residuos se distribuyen hogéneamente. que la varianza no es constante en todo el rango de x. Los límites inferiores y superiores se calculan restando y sumando al valor del parámetro el valor obtenido al multiplicar su error por t(0. Se calcula un valor de t como el valor del parámetro dividido de su desviación estándar. Jurado
Intercepción se refiere a la ordenada y X a la pendiente.J. es decir.
1. Jurado d) Si existen outliers.LINEAL del menú “Insertar /Función” de EXCEL
3. El valor de F se corresponde con el anteriormente descrito y df son los grados de libertad para los residuales. En cuanto a ssreg. la ordenada a y su desviación estándar sa. puesto que se trata de una formula matricial.7. El parámetro aquí llamado ssres es la suma de los residuales al cuadrado (Ecuación 7). es la suma de cuadrados de regresión.
La pendiente b y su desviación estándar sb.J. La función =ESTIMACION. En primer lugar hay que seleccionar un rango de 10 celdas tal y como se indica en la figura. El siguiente paso es seleccionar la fórmula ESTIMACION. En esta figura se ha escrito en las celdas adyacentes a las del rango el parámetro que la formula devolverá.LINEAL() Esta función se utiliza de una forma especial.5. el coeficiente de correlación al cuadrado r2 (para obtener r realizar la raíz cuadrada) y la desviación estándar de regresión sy/x son los parámetros fundamentales que se han explicado en la sección 1.
En “Estadistica” se pondrá un 1 o VERDADERO para que el calculo de ssreg y ssres sea de la forma anteriormente explicada. sino CTRL+MAYÚS+ENTRAR.
. en “Conocido_x” los de x. Jurado
En “Conocido_y” se debe seleccionar los valores de y. M.J. En “Constante” se debe poner un 1 o VERDADERO para que calcule la ordenada en el origen. para que la fórmula de resultados en todas las celdas de la matriz. En vez de pulsar aceptar.
si2 las varianzas de cada una de estas series. Siendo smax la varianza
sospechosa.J. que presenta una varianza sospechosa. es decir. normalmente para un 95% de nivel de confianza. El estadístico C de Cochran se calcula según la siguiente fórmula:
2 smax
2 El valor obtenido se compara con el tabulado C (h. puesto que los datos no son homocedásticos y no se debería aplicar la regresión lineal ordinaria. se calcula una varianza de error puro como la media de las varianzas a cada nivel (suponiendo h niveles). si no hay outliers dispersivos.
Para poder promediar las varianzas de repetición en cada nivel de concentración es necesario que estas no difieran entre sí.2. se considera que la serie de replicados sospechosa constituye un outlier dispersivo. Linealidad: Prueba de la falta de ajuste (lack of fit). r el número de replicados y α el nivel de significación. De esta forma se decide si dicha varianza puede considerarse o no un outlier. por ser la mayor de la serie. Jurado
3. se puede usar la prueba de Cochran [17]. M. Se realiza esta prueba para ver si el modelo se ajusta bien a los datos. y se sospecha que una serie de replicados puede constituir un outlier dispersivo. En este caso no tiene sentido la prueba de la falta de ajuste. Cuando se tiene una serie de resultados agrupados en replicados. En el ANEXO 2 se presentan los valores tabulados del parámetro C.5. α ). Como requisito indispensable haber realizado replicados por cada nivel de concentración.
. h el número de series de replicados. Una vez hecho esto. sino la ponderada. Supongamos que realizamos r replicados en n niveles de concentración. r. En líneas generales se trata de corregir la varianza de los residuales con la varianza de repetición pura promediada para todos los niveles y comprobar si esta varianza corregida (la
2 varianza de falta de ajuste sLOF ) difiere de la de repetición pura mediante un test F. Si el valor calculado es mayor que el teórico. Para ello se realizaría el test de Cochran para detección de outliers dispersivos.
Calculamos la mediana de los factores de respuesta. Se representa y/x frente a la concentración de analito (o frente al logaritmo si el intervalo de trabajo es amplio.
. 3.3.05. A partir del ANOVA se tiene la varianza de los residuales: ˆ ∑( y − y )
i i i n 2
2 sRES =
Para obtener la varianza de falta de ajuste se resta la de error puro a la de los residuales. por lo tanto.95M.05M y 0. trazando una recta en este valor y dos rectas más con valores 1. Jurado
Con r-1 grados de libertad (se han sumado h varianzas de r-1 grados de libertad pero se ha dividido entre h). Así. por ejemplo 1-1000 mg/L).J. Linealidad e intervalo dinámico: método del factor de respuesta de Huber. el modelo lineal se ajusta perfectamente a los datos experimentales. El método del factor de respuesta de Huber [18] se aplica cuando tenemos rectas de calibrado corregidas con el blanco para que la ordenada sea prácticamente cero. si la respuesta es prácticamente y=bx. Esto quiere decir que no existe un error debido a la falta de ajuste del modelo que no pueda ser atribuido al propiamente aleatorio. r − 1) la varianza de falta de ajuste no será estadísticamente mayor que la de error puro. M. M. Los puntos experimentales que se encuentren entre esas dos líneas corresponden al intervalo lineal. n − r − 1.5.
2 sLOF = 2 2 (n − 2) sRES − (r − 1) sPE (64) (n − 2) − (r − 1)
Tras esto se calcula el factor F de falta de ajuste:
FLOF =
2 sLOF (65) 2 sPE
Si FLOF < Fcrit (0. los factores de respuesta R=y/x serán próximos a b.
Si b = 1 y a ≠ 0. ambos métodos llevan a resultados estadísticamente iguales.05. Los métodos son equivalentes. Comparación con un método de referencia. pero en una cantidad fija. el método a calibrar da un resultado más alto o más bajo que el de referencia. Jurado
3. se producen errores sistemáticos constantes o aditivos. Se calculan los parámetros t de Student:
ta = a−0 (66) sa
b −1 (67) sb
Se comparan con el valor tabulado t(0. Si los dos métodos llevan a resultados idénticos. Donde y son los resultados obtenidos con el método a validar y x los obtenidos por el de referencia. Para comparar un método que está siendo validado con uno de referencia usando una recta de regresión se representan los resultados obtenidos por el método a validar frente a los obtenidos por el de referencia.n-2). es decir. es
obvio que la pendiente debe ser igual a la unidad y la ordenada igual a cero. 41
.J. M. obteniendo la ecuación de la recta
y = (b ± sb ) x + (a ± sa ) . Si b=1 y a= 0.4.5.
Si b ≠ 0 y a ≠ 0. El efecto matriz se debe.J. existen los dos tipos de errores. Las muestras. Se trataría de obtener la pendiente de calibrado externo y la de adición patrón y compararlas. no debe ser significativamente distinto de 1. Si es así. 42
. no podremos realizar el el análisis utilizando un calibrado externo. Se puede detectar la presencia de efecto matriz comparando el calibrado externo con el de adición patrón. como su propio nombre indica a la matriz de la muestra. respectivamente. M. su cociente. salvo en contadas ocasiones. sbCE y sbAP son sus desviaciones estándar correspondientes. La medida se realiza con una técnica determinada y la señal obtenida se interpola en los valores obtenidos mediante una recta de calibrado. Si ambas pendientes son prácticamente iguales. Se debe estudiar entonces si existe o no dicho efecto. deben tratarse antes de ser medidas. R. Si tcal es mayor que 2.5. el hecho de tratar la muestra no es suficiente. se produce un error sistemático proporcional a la concentración de analito en el método a validar.6. En primer lugar realizamos el cociente entre ambas pendientes y calculamos la desviación estándar asociada al mismo:
R= bCE (68) bAP
2 2  sbCE sb S R = R  2 + 2AP  bCE bAP 
Donde bCE y bAP son las pendientes de las rectas de calibrado externo y adición patrón. las dos pendientes no seran comparables y el análisis deberá realizarse por el método de adición patrón. 3. Efecto matriz. sino que deberemos usar el método de adición patrón. Jurado
Si b ≠ 1 y a = 0. A veces. pues puede quedar en una matriz que difiere de los estándares utilizados para realizar el calibrado y es muy posible que esa matriz interfiera en el análisis. Calculando la t de Student como:
R −1 (70) SR
El valor calculado se compara con el factor de cobertura k = 2 para un 95 % de nivel de confianza [10].
Lengyel. 1231 6. 624-632. Trends in Analytical Chemistry. L. González. Miller. P. Harter. 29. Method Validation Program (OMA/PVM Department). 8. 1960. R. Ure. Tables of range and studentized range. 44. Multiple comparison methods for means. Miller.M. 12. 227-238. Pablos. J. 63. M. 11. M. UK 1998. 26. C.
13. AOAC International. 2002. Teddington.pdf. EURACHEM. Alcázar. González. Herrador. A. J. G. J. Garantía de la calidad en los laboratorios analíticos. including measurement uncertainty and accuracy profiles. Huber. http://www. 2. A. 67A-76A. 259-278. Jurado
1. Estadística y quimiometría para química analítica. M. Bosque. 3. 2004.
10. J. The fitness for purpose of analytical methods. A. Talanta. W. A. A. Madrid. 96–105 7. T. Braselton. A practical guide to analytical method validation. A laboratory guide to method validation and related topics. Rios.eurachem. 2007. Annals of Mathematical
Statistics. How to validate analytical methods. M. J. Anal. 1998 5. including Appendix D: Guidance for collaborative study procedures to validate characteristics of a method of analysis.
11. Compañó. 2002. Prentice Hall. 31. Abell. 2002. LGC. Lett.J. L. 1998. A. 2000
. G. J. Madrid. Trends in Analytical
Chemistry. M. H. R. Garcıa and J. J. L. M. Rafter. Wood. Martín. 1996. Cuadros. Evaluation of Analytical Methods used for regulations of foods and drugs. 1982. 1122-1147. Horwitz. L. 54. Inczedy. Blackwell Science. A. N. 297-302 4.
9. Analytical Chemistry. Jurado. 18. F. A. Compendium of Analytical Nomenclature 3rd edition.M. LC-GC Int.org/guides/valid.
SIAM Review. Editorial Síntesis.
47-61. G. 18. Pitchard. E. M. I. Guide to te use of repeatability. E. Wiley. Validation of analytical methods: review and strategy.
. Switzerland. ISO/DTS 21748.
2004. 1995. 535. M.
15. W. Annals of Eugenics. 2003 17. Jurado
14. ACOL Series. 23. Trends in Analytical Chemistry. The distribution of the largest of a set of estimated variances as a fraction of their total. Huber. 96-105. UK. Chichester. Cochran. Quality in the analytical chemistry laboratory.J. L. 1998. Taverniers. West Sussex. Geneva. International Organization of Standarization (ISO). van Bockstaele. 11. 16. De Loose. 11. 1941. ISO. reproducibility and trueness estimates in measurement uncertainty estimation. LC-GC
ANEXO 1. Tabla de valores estandarizados para el test de Tukey
k: conjunto de datos. M. Parámétro C para el contraste de Cochran. n: replicados. 95% n.c. Jurado
ANEXO 2.J.
More From This UserGestión de ResiduosCont LumínicaCont Acusticarelleno sanitario mava.pdfTLYA HIDRÓGENOTlya Solar FotovoltaicaEnergía SolarEnergía Oceánica Energía HidráulicaAprov Energía Geotérmica - Energía HidráulicaControl de emisiones Almacenamiento RSU 2016Segregación de RS Parques eólicosPrincipios Energía EólicaAire - ContaminantesAerogeneradoresMetodos naturales TLYA Biomasa-AprovechamientoTLYA Biomasa-GeneralidadesTLYA Producción Limpia Caracterización de Residuos SólidosRESIDUOS SÓLIDOS, CLASIFICACIÓN Y SITUACIÓN ACTUALTecnologías Limpias y Alternativas 1Residuo Sólido
AQAEXCEL Validacion de Metodos by Mercedes Villa33 viewsEmbedDownloadRead on Scribd mobile: iPhone, iPad and Android.Copyright: Attribution Non-Commercial (BY-NC)List price: $0.00Download as PDF, TXT or read online from ScribdFlag for inappropriate contentMore informationShow less
RelatedRiesgos en el uso de Quimicosby Pelado Luis AndresAFOROby Mercedes VillaCurtido_cuerosby eliseodelapenaTeoria de Medicion de Caudalby api-3832818Huella Ecologicaby Miguel Gonzalez SaborioCuero, Calzado y Marroquineriaby edithVaporby Jorge Alberto CruzMediciones Electricas - Juan Sabato - En Españolby alvarado02Instrumentos y Mediciones Eléctricasby lautaroezequielrinaldiHuella de aguaby L Carolina UscangaHuellaby René Alberto Martínez MagdalenoHUELLA DE CARBONOby Alexandra Ake SVEl trenzador de tientosby jacuña_67x6pinneapleby gordonsanweycueroby Patricio Toledo EstradaCaso Exitoso de Producción Más Límpia en la ciudad de Medellínby gonabeTalabartería y Cesteríaby Jonathan MedranoHuella de Carbonoby Javier FajardoManual Práctico de Curtido Natural de Cueros y Producción de Artesaníasby Permacultura501Manual Gestion de La Calidad Ambientalby Arquitecto JAMES DAVID ROJAS QUISPELa Materia Primaby Paola Meravigliacueroby kote_yesCurso Trabajo del cueroby fenix_rsTrabajo del cueroby seancarArtesanias en Cueroby maria547Artesania en Cueroby Nelson David Beltranquimica analiticaby Nicolas NuevaCreager, Caucasiana y Lechoby Mercedes VillaSimilar to EXCEL en química analíticaRiesgos en el uso de QuimicosAFOROCurtido_cuerosTeoria de Medicion de CaudalHuella EcologicaCuero, Calzado y MarroquineriaVaporMediciones Electricas - Juan Sabato - En EspañolInstrumentos y Mediciones EléctricasHuella de aguaHuellaHUELLA DE CARBONOEl trenzador de tientos7x6pinneaplecueroCaso Exitoso de Producción Más Límpia en la ciudad de MedellínTalabartería y CesteríaHuella de CarbonoManual Práctico de Curtido Natural de Cueros y Producción de ArtesaníasManual Gestion de La Calidad AmbientalLa Materia PrimacueroCurso Trabajo del cueroTrabajo del cueroArtesanias en CueroArtesania en Cueroquimica analiticaCreager, Caucasiana y Lechocurso de cueroel arte gaucho del cuero crudo hilario faudoneEXCEL en química analítica

References: Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución