19

Universitas,Volumen 1, Número 2, 2008, 19-26

© 2008 UNAN-León, Editorial Universitaria

Un Experimento sencillo para evaluar la Incertidumbre siguiendo la guia GUM ISO 1995
y utilizando el cálculo simbólico MAPLE 11.0

Gustavo Delgado1* y Bertram Nagel2

1
Facultad de Ciencias y Tecnología, Departamento de Química, Laboratorio de Análisis de Trazas de Metales
Pesados (LATMP), Edificio de Ciencias Básicas, León, Nicaragua. Tel: 505 311 5013, ext: 1132. Fax: 505 311 4012.
Universidad Nacional Autónoma de Nicaragua, León (UNAN-LEÓN).
2
EQL, Consultor PTB (Instituto Alemán de Metrología), Leipzig, Alemania. [email protected]

RESUMEN
El objetivo del presente estudio, es estimar la incertidumbre en la medición del área total de un triángulo aplicando la guía de estimación de la
incertidumbre GUM ISO 1995 y utilizando el programa de cálculo simbólico MAPLE 11. La medición la realizó un grupo de 10 estudiantes de la
maestría en Química y Gestión de la Calidad, utilizando el mismo objeto de medición y el mismo instrumento de medición. Se derivó el modelo
matemático, se identificaron los componentes de las incertidumbres a través del diagrama de causa-efecto, se calcularon los coeficientes de
sensibilidad utilizando el programa de cálculo simbólico Maple 11, se elaboró un balance de las incertidumbres y un diagrama de Pareto donde
se reflejan las variables que más influyen en la incertidumbre del mensurando. Se calcularon las incertidumbres combinada y expandida con
un factor de cobertura igual a t al 95% de nivel de confianza y 22 grados de libertad efectivos. Finalmente, se evaluó el intervalo de confianza,
encontrándose un valor de (50.72 ± 0.39) cm2. La incertidumbre expandida relativa fue menor del 1%, indicando buena precisión en el proceso
de medición.

Palabras claves: Incertidumbre en la medición del área de un triángulo; Cálculo de coeficiente de sensibilidad utilizando Maple;
Incertidumbre expandida relativa (uER); Metrología química; Balance o presupuesto de las incertidumbres.

1. INTRODUCCIÓN del cumplimiento con la especificación del producto

dentro de los límites establecidos [5] en los acuerdos
La armonización de todas las unidades al Sistema comerciales. Estos límites, requieren que cada
Internacional de unidades SI, es de vital importancia laboratorio de ensayo disponga de un procedimiento
en los tratados comerciales. Esta ha sido una de documentado para evaluar la incertidumbre[6-9] en los
las principales barreras técnicas que obstaculiza procesos de medición del parámetro o mensurando
el comercio internacional. La incongruencia de las que indica la calidad del producto.
unidades ha provocado grandes pérdidas para un
país y ganancias para otros. Inclusive, en una misma Dentro de este contexto, la metrología química
nación, si las incongruencias son graves pueden ser desempeña un papel muy importante [10,11] y el
causa de procesos judiciales. Por eso, actualmente personal técnico de un laboratorio de ensayo o de
en los tratados comerciales los países pobres están calibración deberá tener la competencia para evaluar
implementando leyes nacionales para aplicar el uso el mensurando y estimar su incertidumbre asociada[4,5].
del SI basado en el sistema métrico decimal y en sus Con el fin de introducir a los estudiantes de ciencias
unidades básicas, derivadas y suplementarias[1]. químicas y físicas en el tema de la metrología, se
ha reportado el análisis de la incertidumbre en los
Así por ejemplo, en la exportación de productos cálculos del área y el perímetro de un rectángulo
de valor comercial es necesario el cumplimiento utilizando el análisis estadístico clásico y el cálculo
con especificaciones del país comprador basado del error máximo probable[12]. En el presente estudio,
en el sistema SI. De esta manera, el producto va se pretende evaluar el área de un triángulo y estimar
acompañado con un informe de la conformidad[2,3] su incertidumbre siguiendo la guía GUM ISO 1995[6] y
extendido por los organismos competentes (organismos utilizando el programa de cálculo simbólico MAPLE[13].
evaluador de la conformidad, laboratorios de ensayos
acreditados), especificando el valor del parámetro Las mediciones la realizaron un grupo de 10
requerido o mensurando con su incertidumbre a un estudiantes de la maestría en Química y Gestión de la
nivel de confiabilidad establecida [4]. Los compradores Calidad, sobre la misma figura geométrica y utilizando
disponen de una guía donde se describe el método que el mismo instrumento de medición, una regla graduada
tienen que cumplir los laboratorios de ensayos o de en centimetros.
calibración para declarar los resultados y la evaluación
* Autor para correspondencia: [email protected]
20
2. TEORÍA Las derivadas parciales se denominan coeficientes
de sensibilidad[6] (cxi); uxi son las incertidumbres de las
Ley de propagación de la incertidumbre variables xi.

Un proceso de medición puede representarse de Si en el modelo matemático intervienen los parámetros de

acuerdo al “enfoque de proceso” establecido en la norma regresión de un modelo de calibración lineal (intercepto
ISO 9000:2000[14]. Pueden identificarse tres partes: y pendiente) la ecuación (2) se puede expresar de la

variables de entrada x1, x2, ... xn; relación funcional entre siguiente manera:

( )
 variables de entrada; el mensurando y como variable
las
u (2y ) =  c xi u xi + 2(cb 0 )(cb1 )u (b 0 ) u (b1) r(b 0,b1)
2
de salida. El esquema se puede ver en la figura 1
(3) 
Si se considera que no existen variables correlacionadas

Variables de
entrada xi:
Función:
Variable de salida
o Mensurando y:
la incertidumbre uy se simplifica
 a la ecuación (4).
f(x1, x2,...xn)
x1, x2,...xn y=f(xi)

 (c )
2

Figura 1. Proceso de medición de un mensurando. u( y ) = xi u xi ( 4) 
  Para evaluar la incertidumbre expandida, es necesario
Las variables de entrada son aleatorias que tienen una
 el factor de cobertura k. Este
seleccionar  factor puede
distribución de probabilidad conocida y pueden ser tomar valores como los que se muestran en la tabla 1.
independientes o estar correlacionadas entre si[15].
TABLA 1. Factor de cobertura, nivel de confianza y tipos de distribución
de probabilidad
El proceso de medición puede así representarse por un
modelo matemático que relaciona el mensurando y las k Nivel Confianza Tipo de
variables de entradas a través de la función f, tal a como distribución
se muestra en la ecuación (1); se asume ausencia de 2 95.45% Normal
errores sistemáticos[15]. El segundo término representa
1.96 95.% Normal
la incertidumbre expandida del mensurando[6, 9].
Donde Uy=kuy; uy es la incertidumbre combinada; k t0.975,v 95% Student, 2 colas

y = f ( x1 , x2 ,...xn ) ± U y (1) 
Generalmente, cuando se reporta la incertidumbre se
es el factor de cobertura, que depende de los grados selecciona k=2 para 95.45 % de nivel de confianza[7].
de libertad [6-9] y de las diferentes distribuciones
 de Pero si se quiere efectuar una selección rigurosa de k
probabilidad de las variables de entradas implicadas en la es necesario evaluar los grados de libertad efectivos νef
incertidumbre combinada. Para evaluar la incertidumbre [6]
; si éste tiene un valor menor de 30, el factor k sigue
combinada uy se aplica la ley de propagación de la una distribución t de Student. Los grados de libertad
incertidumbre[5,15] a la función f. Esta ley se deduce efectivos se calculan de acuerdo a la ecuación (5).
aplicando la expansión de Taylor [11, 15], tomando uy =∆y 4
u uy
o diferencia entre la función f (xi) y f evaluada en los ν ef = (5) 
promedios de xi. La expresión final que representa esta c i4 u xi4
ley es la siguiente:
 ν
i

Si las variables (xi) son independientes, el segundo cxi son los coeficientes de sensibilidad; uxi son las
2 
incertidumbres 
de la variable xi y νi son los grados
 ∂f  2 N −1 N
 ∂f  ∂f 
u (2y ) =    u xi + 2    u ( xi , xj ) (2) de libertad para la incertidumbre de cada variable. La
 
 ∂xi  i =1 j =i +1 ∂x i  ∂x j  tabla 2 muestra los grados de libertad para cada tipo de
distribución.
término denominado covarianza es despreciable. En los
modelos lineales los coeficientes del modelo (bi) están TABLA 2. Grados de libertad vi en función del tipo de distribución [6,
20]
.
correlacionados y debe tomarse en cuenta en el cálculo
de la incertidumbre combinada. Tipo de distribución vi
Normal 50
Así por ejemplo, los coeficientes intercepto y pendiente
del modelo de una línea recta (b0 y b1) tienen una Rectangular ∞
correlación negativa y debe calcularse el coeficiente Triangular ∞
de correlación (r(b0,b1)) para ser introducido en el término t-Student (univariante) ni-1
covarianza[6-9].
t-Student (bivariante) ni-2
t-Student (p variables) ni-p
 21


 
 

 



  
   

  

  



 
 

 
 

 



  




Figura 2: Diagrama de flujo representando el procedimiento para la estimación de la incertidumbre en un proceso de medición.

Conociendo el valor k se puede estimar la incertidumbre El proceso de calibración, se debe llevar cabo bajo
expandida y el intervalo de confianza: y ± Uy. condiciones ambientales específicas (temperatura y
humedad relativa, medidas también con instrumentos
Procedimiento para la evaluación de la calibrados trazables). Si se asume que, los valores
incertidumbre de las variables de entrada han sido obtenidos con
instrumentos de medición calibrados con patrones
Existe un procedimiento para estimar la incertidumbre[6,9], trazables al SI, el mensurando y de la ecuación (1) es
que se presenta en el diagrama de flujo de la figura 2. también considerado trazable al Sistema Internacional
En primer lugar, se define el mensurando, se deriva
[16]
.
el modelo matemático, se identifican las fuentes de
incertidumbres (las de tipo A, obtenidas a partir de la Si las condiciones ambientales (temperatura) durante
desviación estándar de una serie de observaciones y el proceso de medición del mensurando, son diferentes
las de tipo B, obtenidas a partir de las distribuciones de las de calibrado del instrumento de medición, es
normal, rectangular o triangular[6-9], se cuantifican necesario realizar las correcciones necesarias. Por
las incertidumbres, se calculan los coeficientes de ejemplo, si la regla graduada en centímetros es el
sensibilidad, se hace un balance o presupuesto de instrumento de medición, la corrección a realizar debe
las incertidumbres, se combinan las incertidumbres a tomar en cuenta el coeficiente de dilatación lineal, cuyo
través de la ley de la propagación de la incertidumbre, valor dependerá de la naturaleza de material del cual
se revisa si no hay incongruencias para reevaluar los está hecho este instrumento, considerando la diferencia
componentes de la incertidumbre, se selecciona el factor entre la temperatura de medición y la de calibrado.
de cobertura, se evalúa la incertidumbre expandida y se
reporta el intervalo de confianza con el número correcto 3. PARTE EXPERIMENTAL
de cifras significativas.
Materiales
Trazabilidad
• Computadora equipada con el programa de cálculo
Todo instrumento de medición debe ser calibrado con simbólico MAPLE 11[13].
patrones certificados trazables al Sistema Internacional
de Unidades[4,16]. Para medir la longitud de un segmento • Instrumento de medición: regla de madera
es necesario que el instrumento de medición, por ejemplo impermeabilizada de 0.30 m de longitud, graduada en
una regla graduada en centímetros, esté calibrado con centímetros.
patrones certificados trazables a patrones nacionales y
éstos a su vez trazables al SI. • Coeficiente de dilatación lineal del orden de 10-6/
C .
o [17]
22
• Objeto de medición: como el representado en la figura Definición del mensurando
geométrica del anexo.
El mensurando, es el área total de la figura geométrica
• Termómetro digital calibrado TESTO 925, con del anexo, constituida por los dos triángulos rectángulos,
resolución 0.1ºC. y se designa simbólicamente como A.

4. Procedimiento Modelo matemático

Resolución del instrumento de medición El modelo para el área total del objeto de medición, es
la suma de las áreas de los dos triángulos. La expresión
La resolución de la regla graduada en centímetros fue matemática se presenta en la ecuación (6) admitiendo
evaluada por simple observación visual. El valor mínimo ausencia de errores sistemáticos:
leído en la escala fue 0.05 cm.
c
A= (b + d ) ( 6) 
Proceso de medición 2
Identificación de los componentes de la
El objeto de medición y el instrumento de medición incertidumbre
fueron proporcionados a un grupo de 10 estudiantes.
Cada estudiante midió los segmentos de recta de Los componentes de las incertidumbres para cada
la figura geométrica, que representa dos triángulos variable del modelo se muestran en el diagrama
utilizando la misma regla y ubicando con precisión, el causa-efecto de la figura 3. Cada variable tiene dos
cero y el valor del extremo superior del segmento. Las componentes: la incertidumbre debido a la resolución
mediciones fueron realizadas tratando de minimizar el de la regla y la de la serie de observaciones de los 10
error de paralaje y los resultados se registraron según la estudiantes.
matriz de datos de la tabla 3. La variación de temperatura
durante el proceso de medición fue de 2ºC.

TABLA 3. Diseño de los resultados de la medición en cm de los

segmentos a, b, .. de la figura geométrica en anexo.

No. a b ......
1
2
.
.
.
10

Figura 3. Diagrama causa-efecto para componentes de la incertidumbre
4. RESULTADOS Y DISCUSIÓN en la medición del área de la figura geométrica del anexo.
 
Los resultados que a continuación se presentan, Cuantificación de los componentes de la
siguen el procedimiento descrito en el flujograma incertidumbre
de la figura 1. Dado que, el coeficiente de dilatación
lineal del instrumento de medición es muy pequeño, a) Incertidumbre en la resolución uR
del orden de 10-6 por cada ºC, el efecto de la variación
de temperatura, no influyó significativamente en las La incertidumbre en la resolución (de tipo B) sigue una
mediciones de los segmentos del triángulo. Como la distribución rectangular[15], la cual se expresa así:
regla está impermeabilizada, se asume que la humedad
relativa no afecta la medición. Se admite que, la regla R
uR = (7 )  
está calibrada y que el error sistemático, es mucho 2 3
menor que el error aleatorio.
Sustituyendo R=0.05 cm se tiene un valor de uR=
0.0144
 23
b) Incertidumbre en la serie de observaciones uo El subíndice xi representa al segmento b, c o d. R es la
resolución de la regla graduada. sxi es la incertidumbre
La incertidumbre de una serie de observaciones (de estándar del segmento leído. n es el número de
tipo A) está dada por la desviación estándar de la media observaciones. Sustituyendo los valores de R=0.05 cm
o clásicamente conocida como error estándar[18]. La y los valores de sxi y n de la tabla 4 en las ecuaciones
ecuación (8) representa la fórmula para la incertidumbre (7), (8) y (10), se tienen los valores de incertidumbres
de las observaciones. en las mediciones b, c y d que se muestran en la tabla
s 5.
uo = s x = (8)  
n TABLA 5. Incertidumbres y grados de libertad para las variables b, c
y d.
Donde sx es la desviación estándar de la media, s es
xi uR(v) uo(v) uxi(vi)


la desviación estándar de las observaciones y n es el

número de observaciones igual a 10. Los resultados b 0.0144 (∞) 0.00748331 (9) 0.021794 (9)
de las 10 observaciones del grupo de estudiantes se c 0.0144 (∞) 0.00763763 (9) 0.021741 (9)
presentan en la tabla 4.
d 0.0144 (∞) 0.00816497 (9) 0.021985 (9)
TABLA 4: Valores de los segmentos b, c, d en cm, media, desviaciones
estándares y grados de libertad vi. Lecturas de 10 estudiantes. Los grados de libertad vi se calculan combinando
los grados de libertad de R y o según la expresión
Estudiante b c d representada en la ecuación 5. De acuerdo a la tabla
1 8.3 7.85 4.55 2, los grados de libertad para uR con una distribución
triangula[15] tienen un valor infinito y los grados
2 8.3 7.9 4.6
de libertad para uo obtenidos para una serie de
3 8.3 7.9 4.6 observaciones tienen una distribución t a un nivel de
4 8.3 7.9 4.6 confianza de 95% y n-1 grados de libertad.
5 8.25 7.9 4.6
Balance de las incertidumbres
6 8.25 7.85 4.6
7 8.3 7.85 4.55 El balance de las incertidumbres establece las
8 8.29 7.9 4.55 incertidumbres de las variables del modelo uxi, los
coeficientes de sensibilidad cxi y el índice de contribución
9 8.3 7.9 4.55
de cada variable xi en la incertidumbre de la función del
10 8.25 7.9 4.6 modelo o el mensurando buscado, que se representa
Media 8.284 7.885 4.58 como porcentaje (Ind%). Los coeficientes de sensibilidad
s 0.02366432 0.02415229 0.02581989 son la derivada parcial de la función A con respecto a
cada variable xi.
uo=s/√n 0.00748331 0.00763763 0.00816497
v 9 9 9 Cálculo de los coeficientes de sensibilidad cxi:
c) Incertidumbre de las variables b, c y d (uxi)
∂A
c xi = (11) 
La incertidumbre en las mediciones de los segmentos b, ∂ xi
c y d, se puede ver en la ecuación (9). Las derivadas parciales se pueden calcular fácilmente,
dado que la función A de la ecuación (6) es simple.
u xi = 2u R2 + u o2 (9)  
Cuando las funciones son complejas el cálculo simbólico
Maple 11 simplifica el tiempo de cálculo. Además la
Generalizando se tiene la siguiente ecuación (10): precisión es ilimitada y no utiliza redondeos. Aplicando
1 36 s xi Maple para este ejemplo se tienen los siguientes
u xi = 6R 2 + (10)  resultados formales: 
6 n

El coeficiente 2 del primer término en la raíz de la Se escribe la función A:
ecuación (9) representa la doble lectura en los extremos
 c
del segmento xi, el cero leído en el extremo inferior y el > A := * (b + d ); 
valor de xi (en cm) leído en el otro extremo. 2
c
(b + d ) 
 2 

 
24
Se deriva la función A con respecto a los segmentos La tabla 6 muestra el balance de las incertidumbres,
c, b y d, utilizando las expresiones
 preestablecidas de donde se pueden ver las fuentes de incertidumbres (b,
derivadas parciales: c y d), el valor de cada variable xi, las incertidumbres
∂ y los coeficientes de sensibilidad correspondientes, los
> c1 := A; 1 1 valores absolutos y los cuadrados del producto de los
∂c b+ d coeficientes por las incertidumbres y por último, el índice
2 2
de contribución de la incertidumbre de cada variable en

la incertidumbre del mensurando.
 
∂ La contribución de cada variable a la incertidumbre
> c 2 := A; 1
∂b c del área total de la figura geométrica se evalúa con el

2  índice de contribución expresado en porcentaje según
la ecuación siguiente:
 
∂ (c xi u xi ) 2
> c3 := A; 1 Ind % = (100) (12) (12) 
∂d c  xi xi
( c u ) 2
2 
Si trazamos el diagrama de Pareto [19] se tiene la
  figura 4, donde se puede ver que la variable que
donde c1, c2 y c3 indican los coeficientes de sensibilidad más contribuye a la incertidumbre es la medición en
para c, b y d. Asignando los valores promedios c, b y d a el segmento c. No obstante, las tres mediciones de
partir de la tabla 4, se tienen los resultados numéricos los segmentos contribuyen significativamente a la
de cxi. incertidumbre en el mensurando A, admintiendo como
límite el 90% de la contribución acumulativa.
100
Estimación de la incertidumbre combinada uA
80 La incertidumbre de A (uA) es la combinación de
todas las incertidumbres que influyen en la medición
60
Ind%

del mensurando A. Se asume que las variables de la

función son independientes y se aplica la ecuación (4).
40 La incertidumbre del mensurando A es la raíz cuadrada
de la suma de los productos (cxiuxi)2 en la columna 6 de
20 la tabla 6, cuyo valor es:

0 uA = 0.1858
c d b
Figura 4. Diagrama de Pareto indicando la contribucion de las
Selección del factor de cobertura k
incertidumbres en la medición del área total del triangulo.
Si se observa la tabla 5, la distribución t es la
predominante
 en la estimación de uxi. Por lo tanto, si se
aplica la ecuación (5) y tomando los datos de la tabla 6,
se obtienen los grados de libertad efectivos:

TABLA 6. Balance de las incertidumbres en la medición del área de un triángulo como

objeto de medición.
Fuente de
incertidumbre Valor de mxi cxi /cxi mxi/ (cxi mxi)2 Ind%
x xi
i
b 8.284 0.021794 3.9425 0.08571 0.0073468 21.29
c 7.885 0.021741 6.432 0.14018 0.0196651 56.94
d 4.580 0.021985 3.9425 0.086675 0.0075126 21.77
Suma= 0.0345104 100
 25
(0.1858) 4
función área con respecto a cada una de las variables o
ν ef = segmentos que constituyen el triángulo. Para evaluar la
(3.9425) (0.021794) (6.432) 4 (0.021741) 4 (3.9425) 4 (0.021985) 4 
4 4
+ + incertidumbre combinada en el proceso de medición, se
9 9 9 aplicó la ecuación simplificada de la ley de propagación
de la incertidumbre, considerando que las variables son
νef = 22 independientes.

El valor de k es igual a t al 95% de nivel de confianza a El intervalo de confianza del área del triángulo, se
dos colas y 22 grados de libertad[6,18]. evalúo la incertidumbre expandida utilizando el factor
de cobertura k como el valor de t al 95 % de nivel
k=t0.975,22 = 2.074 de confianza y 22 grados de libertad efectivos. El
valor encontrado del intervalo de confianza fue de
Evaluación de la incertidumbre expandida (50.72 ± 0.39) cm2.

La incertidumbre expandida, se calcula multiplicando el 6. Reconocimiento

factor k por la incertidumbre combinada UA. El cálculo
se muestra a continuación: Los autores agradecen a Gustavo Montiel, Director
del Laboratorio Nacional de Metrología de Nicaragua
UA= k uA=(2.074)(0.1858)=0.38535 cm2 (LANAMET) y a Miguel Carrero Gálvez, Consultor de
Asesores de Calidad y Medioambiente (CACM-España),
Estimación del intervalo de confianza del área del por sus valiosos comentarios y sugerencias al presente
triángulo trabajo.

El intervalo de confianza del área del triángulo con una 7. Bibliografía

probabilidad del 95% y tomando en cuenta el número de
cifras significativas correcto, se calcula de la siguiente 1. LEY N° 225 (1996). “Ley sobre Metrología”, La
manera: Gaceta No. 135, República de Nicaragua.
7.885
y = A ±U A = (8.284 + 4.58) ± 0.38535  2. ISO/IEC 17000 (ES) (2004). “Evaluación de la
2
conformidad — Vocabulario y Principios Generales”,
y= (50.72 ± 0.39) cm2 Ginebra.

La precisión de la medición puede estimarse evaluando 3. ISO/IEC 17011(ES) (2004). “Evaluación de la

la incertidumbre expandida relativa[11] (UER) expresada conformidad — Requisitos generales para los
en porcentaje, la cual resulta de dividir la incertidumbre organismos de acreditación que realizan la acreditación
expandida entre el valor del mensurando o área total de organismos de evaluación de la conformidad”, Norma
multiplicada por 100. El valor calculado de esta manera Internacional, Ginebra
es de 0.77%, lo que indica una buena precisión en
la medición si tomamos como referencia 1% como
4. ISO/IEC 17025(ES) (2005). “Requisitos Generales
máximo.
para la Competencia de los Laboratorios de Ensayo y
de Calibración”, Norma Internacional, Ginebra, Segunda
5. CONCLUSIONES
Edición.
Siguiendo la guía para la estimación de la incertidumbre
establecida en GUM ISO 1995, se aplicó el procedimiento 5. ILAC-G8 (1966). “Guidelines on Assessment and
para evaluar la incertidumbre en la medición del área Reporting of compliance with specification”.
total de un triángulo, realizando las mediciones un grupo
de 10 estudiantes de maestría en química y gestión de 6. BIPM, IEC, IFCC, IUPAC, OIML (1995). Guide for to
la calidad de la Facultad de Ciencias y Tecnología de la the Expression of Uncertainty in Measurement (GUM),
UNAN-León, Nicaragua. ISO, Ginebra.

En el cálculo de los coeficientes de sensibilidad se 7. NIST (1993). “Guidelines for evaluation and
utilizó el programa de cálculo simbólico MAPLE 11, expressing the uncertainty of NIST measurement
como una aplicación de las derivadas parciales de la results”, Gaithersburg, USA.
26
8. SCHMID, W.A. y LAZOS, R.J. (2000), “Guía para 8. Anexo
estimar la incertidumbre de la medición”, CENAM,
México. Figura geométrica como objeto de medición. La escala
de la figura original utilizada durante el experimento,
9. EURACHERM/CITAC GUIDE, (2000), Quantifying puede variar al pasar de MS Word a la configuración
Uncertainty in Analytical Chemistry, 2da edición, UK. del editor de texto utilizado en la revista científica.

10. DELGADO M, VANEGAS M y DELGADO G. (2007),

“Metrología Química I: Calibración de un pH-metro y
control de calidad”, UNIVERSITAS UNAN-León, 1, 14-
20.
a e
11. DELGADO M, VANEGAS M y DELGADO G. (2007)
“Metrología Química II: Estimación de la incertidumbre
en la medición de pH en aguas”, UNIVERSITAS UNAN- c
Léon, 2, 19-26.
Área 1
Área 2
12. YATES P. C. (2001), “A simple for illustrating
uncertainty analysis”, Journal of Chemical Education,
78, 770-771.

13. a) CHAR B. W. (1992), Maple V, Springer-Verlag, b d

New York, pag. 41; b) Maple 11, software.

14. ISO 9000:2000 (ES), (2000), Sistemas de Gestión

de la Calidad. Fundamentos y Vocabulario.

15. KIMOTHI, S.K., (2002), The Uncertainty of

Measurements, ASQ Quality Press, Wisconsin , pag.
205.

16. EURACHEM/CITAC GUIDE, (2003), Treceability

in chemical measurement, Guía de Referencia, 1ra.
Edición, UK.

17. DEAN, J. A. (1979), Lange´s Hanbook of Chemistry,

McGraw Hill, pag. 10-119.

18. NEULLY, M. (1998), Modélisation et estimation des

erreurs de mesure, Technique & Documentation, 2ª
edición, pag. 296.

19. WODSWORTH, H. M., STEPHENS, K. y GODFREY,

A. B. (2005), Métodos de Control de Calidad, México,
Compañía Editorial Continental, 1ra. Edición, p. 362.

20. CHRZANOW (2002). “For evaluation of uncertainty

in mesasurements”, Joint IMEKO TC-1 & XXXIV
Conference 2002, 11- 12 septiembre, Warsaw.