Chemistry · Libro 3 · Bachelor Year 2

Química universitaria — segundo año

Química universitaria — segundo año · Bachelor Year 2

34Estadística de la medida química

Dos laboratorios miden el plomo de la misma agua del grifo e informan de 9.6 y 11.2 µg/L11.2\,\text{µ}\mathrm{g}/\mathrm{L}. El valor guía para el agua de consumo es 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L}. ¿Es potable el agua? ¿Discrepan siquiera los dos laboratorios, o ambos resultados son el mismo valor visto a través de la dispersión de dos medidas? Un número procedente de un análisis no significa nada sin su incertidumbre, y este capítulo da las herramientas para calcularla: la estadística de las medidas repetidas, la propagación de las incertidumbres a través de un cálculo, el calibrado por mínimos cuadrados y los límites por debajo de los cuales un método no ve nada.

Un laboratorio de análisis de aguas: botellas de muestra, una técnica pipeteando, un analizador elemental al fondo (ilustración).
Un laboratorio de análisis de aguas: botellas de muestra, una técnica pipeteando, un analizador elemental al fondo (ilustración).

Lo que ya sabes

El volumen del primer año: incertidumbre de medida, incertidumbre típica, evaluaciones de tipo A y de tipo B, incertidumbre expandida, incertidumbre relativa, desviación normalizada zz entre dos valores, y la promesa de una demostración de u(xˉ)=s/nu(\bar x) = s/\sqrt n y de una regla general de propagación. El volumen escolar: la recta de calibrado. Capítulo 32: la absorción atómica.

34.1 Medidas repetidas

Si se repite una medida, los resultados se dispersan. Tratamos cada resultado xix_i como una variable aleatoria de media μ\mu (el valor que daría el método en promedio) y desviación típica σ\sigma; cuando se suman muchas causas pequeñas e independientes, los resultados siguen la ley normal (gaussiana), con alrededor del 68 % de ellos dentro de μ±σ\mu \pm \sigma y el 95 % dentro de μ±1.96σ\mu \pm 1.96\sigma.

Definición 34.1 (Dispersión)

Para nn resultados x1,…,xnx_1, \dots, x_n de media xˉ\bar x, la desviación típica experimental es s=∑(xi−xˉ)2/(n−1)s = \sqrt{\sum (x_i - \bar x)^2/(n-1)}; es una estimación de σ\sigma. La desviación típica de la media es la desviación típica de xˉ\bar x sobre muchas series de nn resultados, estimada por s/ns/\sqrt n.

El divisor n−1n - 1 en lugar de nn corrige el hecho de que las desviaciones se toman respecto a xˉ\bar x, calculada a su vez con los mismos datos, lo que las hace un poco demasiado pequeñas: n−1n - 1 es el número de desviaciones independientes, los grados de libertad.

Teorema 34.2 (Varianza de la media)

Si x1,…,xnx_1, \dots, x_n son independientes, cada uno de varianza σ2\sigma^2, entonces var⁡(xˉ)=σ2/n\operatorname{var}(\bar x) = \sigma^2/n: la desviación típica de la media es σ/n\sigma/\sqrt n.

Demostración. Para variables independientes la varianza de una suma es la suma de las varianzas, y var⁡(cX)=c2var⁡(X)\operatorname{var}(cX) = c^2\operatorname{var}(X). Así, var⁡(xˉ)=var⁡(1n∑xi)=1n2∑var⁡(xi)=1n2 nσ2=σ2/n\operatorname{var}(\bar x) = \operatorname{var}\bigl(\tfrac1n\sum x_i\bigr) = \tfrac{1}{n^2}\sum\operatorname{var}(x_i) = \tfrac{1}{n^2}\,n\sigma^2 = \sigma^2/n. Sustituir σ\sigma por su estimación ss da la incertidumbre típica de tipo A u(xˉ)=s/nu(\bar x) = s/\sqrt n admitida en el volumen del primer año. ∎

34.2 Intervalos de confianza

Definición 34.3 (Intervalo de confianza)

Un intervalo de confianza con un nivel de confianza 1−α1 - \alpha (a menudo el 95 %) es un intervalo calculado a partir de los datos mediante una regla que, aplicada a muchas series, contiene el valor verdadero en una fracción 1−α1 - \alpha de ellas. Para la media de nn resultados es xˉ±t s/n\bar x \pm t\,s/\sqrt n, donde tt es el coeficiente de Student para n−1n - 1 grados de libertad y ese nivel.

Proposición 34.4 (Ley de Student)

Para nn resultados normales independientes, T=(xˉ−μ)/(s/n)T = (\bar x - \mu)/(s/\sqrt n) sigue la ley de Student con n−1n - 1 grados de libertad, cuya densidad es proporcional a (1+t2/ν)−(ν+1)/2(1 + t^2/\nu)^{-(\nu+1)/2} con ν=n−1\nu = n - 1. Es más ancha que la ley normal y tiende a ella cuando nn crece.

Demostración. Admitido a este nivel. ∎

Dividir por ss en lugar de por el desconocido σ\sigma añade la dispersión del propio ss, que es grande cuando nn es pequeño; de ahí las colas más pesadas y los coeficientes mayores que 1.96. La tabla da los coeficientes al 95 %, calculados a partir de la densidad.

Densidades de Student para 2 y 5 grados de libertad frente a la ley normal, y los coeficientes bilaterales al 95 %, calculados integrando la densidad (la última fila es la ley normal).
ν=n−1\nu = n - 1tt (95 %)
112.706
24.303
33.182
42.776
52.571
92.262
202.086
∞\infty1.960
Densidades de Student para 2 y 5 grados de libertad frente a la ley normal, y los coeficientes bilaterales al 95 %, calculados integrando la densidad (la última fila es la ley normal).

Método 34.5 (Un intervalo de confianza para una media)

  1. Calcúlense xˉ\bar x y ss de los nn resultados; compruébese que ningún resultado se aparta groseramente de los demás.
  2. Léase tt para ν=n−1\nu = n - 1 al nivel elegido.
  3. Exprésese xˉ±t s/n\bar x \pm t\,s/\sqrt n, con la incertidumbre redondeada a una o dos cifras significativas y la media a la misma posición decimal.
  4. Para comparar con un valor de referencia xrefx_{\mathrm{ref}}, calcúlese ∣xˉ−xref∣/(s/n)|\bar x - x_{\mathrm{ref}}|/(s/\sqrt n): por encima de tt, la diferencia es significativa a ese nivel; por debajo, los datos no muestran diferencia (lo que no demuestra que no la haya).

34.3 Propagación de la incertidumbre

Teorema 34.6 (Ley de propagación de la incertidumbre)

Si y=f(x1,…,xn)y = f(x_1, \dots, x_n), donde las xix_i son independientes con incertidumbres típicas u(xi)u(x_i) lo bastante pequeñas para que ff sea casi lineal en ellas, entonces

u2(y)=∑i=1n(∂f∂xi)2u2(xi).u^2(y) = \sum_{i=1}^{n}\left(\frac{\partial f}{\partial x_i}\right)^2 u^2(x_i).

Demostración. A primer orden (Taylor), y−f(μ1,…,μn)≈∑ci(xi−μi)y - f(\mu_1, \dots, \mu_n) \approx \sum c_i(x_i - \mu_i) con ci=∂f/∂xic_i = \partial f/\partial x_i en las medias. La varianza de una combinación lineal de variables independientes es ∑ci2var⁡(xi)\sum c_i^2\operatorname{var}(x_i) (como en el Teorema 34.2); con var⁡(xi)=u2(xi)\operatorname{var}(x_i) = u^2(x_i) se obtiene el resultado. ∎

Corolario 34.7 (Sumas y productos)

Para una suma o una diferencia, las incertidumbres típicas se suman cuadráticamente: u2(x1±x2)=u2(x1)+u2(x2)u^2(x_1 \pm x_2) = u^2(x_1) + u^2(x_2). Para un producto o un cociente y=x1ax2b⋯y = x_1^{a}x_2^{b}\cdots, las incertidumbres relativas se suman cuadráticamente, cada una ponderada por su exponente: (u(y)/y)2=a2(u(x1)/x1)2+b2(u(x2)/x2)2+⋯\bigl(u(y)/y\bigr)^2 = a^2\bigl(u(x_1)/x_1\bigr)^2 + b^2\bigl(u(x_2)/x_2\bigr)^2 + \cdots

Demostración. Para una suma las derivadas parciales valen ±1\pm1. Para y=x1ax2by = x_1^ax_2^b, ∂y/∂x1=a y/x1\partial y/\partial x_1 = a\,y/x_1 y ∂y/∂x2=b y/x2\partial y/\partial x_2 = b\,y/x_2; el teorema da u2(y)=a2y2u2(x1)/x12+b2y2u2(x2)/x22u^2(y) = a^2y^2u^2(x_1)/x_1^2 + b^2y^2u^2(x_2)/x_2^2; divídase por y2y^2. ∎

Método 34.8 (Propagar)

  1. Escríbase el resultado como fórmula de las magnitudes medidas; enumérese cada dato de entrada con su incertidumbre típica (tipo A o B).
  2. Para productos y cocientes, trabájese con incertidumbres relativas; para sumas, con absolutas; en otro caso, tómense derivadas parciales.
  3. Combínese cuadráticamente; búsquese el término dominante: mejorar los demás es un esfuerzo inútil.
  4. Multiplíquese por el factor de cobertura (2, o el tt de Student cuando el término dominante procede de pocas repeticiones) para obtener una incertidumbre expandida.

34.4 Calibrado

Definición 34.9 (Curva de calibrado)

Una curva de calibrado relaciona la señal yy de un instrumento con la concentración xx, a partir de patrones de concentración conocida. La recta de mínimos cuadrados y=a+bxy = a + bx es la recta que minimiza S=∑(yi−a−bxi)2S = \sum (y_i - a - bx_i)^2; cada ei=yi−a−bxie_i = y_i - a - bx_i es un residuo.

Teorema 34.10 (Recta de mínimos cuadrados)

Con xˉ\bar x, yˉ\bar y las medias, Sxx=∑(xi−xˉ)2S_{xx} = \sum(x_i - \bar x)^2 y Sxy=∑(xi−xˉ)(yi−yˉ)S_{xy} = \sum(x_i - \bar x)(y_i - \bar y), la recta de mínimos cuadrados tiene

b=SxySxx,a=yˉ−bxˉ.b = \frac{S_{xy}}{S_{xx}}, \qquad a = \bar y - b\bar x.

Demostración. SS es una función cuadrática de aa y bb. Sus derivadas parciales se anulan en el mínimo: ∂S/∂a=−2∑(yi−a−bxi)=0\partial S/\partial a = -2\sum(y_i - a - bx_i) = 0 da a=yˉ−bxˉa = \bar y - b\bar x; sustituyendo en ∂S/∂b=−2∑xi(yi−a−bxi)=0\partial S/\partial b = -2\sum x_i(y_i - a - bx_i) = 0 se obtiene ∑xi((yi−yˉ)−b(xi−xˉ))=0\sum x_i\bigl((y_i - \bar y) - b(x_i - \bar x)\bigr) = 0, es decir, Sxy=bSxxS_{xy} = bS_{xx} (pues ∑xˉ(yi−yˉ)=∑xˉ(xi−xˉ)=0\sum \bar x(y_i - \bar y) = \sum \bar x(x_i - \bar x) = 0). Las derivadas segundas forman la matriz

(2n2∑xi2∑xi2∑xi2),\begin{pmatrix} 2n & 2\sum x_i\\ 2\sum x_i & 2\sum x_i^2\end{pmatrix},

cuya diagonal es positiva y cuyo determinante 4nSxx4nS_{xx} es positivo cuando las xix_i no son todas iguales: el punto estacionario es un mínimo. ∎

Proposición 34.11 (Leer una concentración)

Para nn patrones con desviación típica residual sy/x=∑ei2/(n−2)s_{y/x} = \sqrt{\sum e_i^2/(n-2)}, una muestra leída mm veces con señal media yˉ0\bar y_0 tiene x0=(yˉ0−a)/bx_0 = (\bar y_0 - a)/b y desviación típica

sx0=sy/xb1m+1n+(yˉ0−yˉ)2b2Sxx,s_{x_0} = \frac{s_{y/x}}{b}\sqrt{\frac1m + \frac1n + \frac{(\bar y_0 - \bar y)^2}{b^2S_{xx}}},

que se usa con el tt de Student para n−2n - 2 grados de libertad.

Demostración. Admitido a este nivel. ∎

Los tres términos bajo la raíz son la dispersión de las lecturas de la muestra, la incertidumbre de la altura de la recta y la de su pendiente; el último se anula en el centro del intervalo de calibrado, que es donde mejor se leen las muestras.

Calibrado del plomo por absorción atómica en horno de grafito a 283.3\, nm (datos de ejercicio, usados en el problema). Izquierda: los cinco patrones y la recta de mínimos cuadrados. Derecha: los residuos, pequeños y sin tendencia: una recta es adecuada. Calibrado del plomo por absorción atómica en horno de grafito a 283.3\, nm (datos de ejercicio, usados en el problema). Izquierda: los cinco patrones y la recta de mínimos cuadrados. Derecha: los residuos, pequeños y sin tendencia: una recta es adecuada.
Calibrado del plomo por absorción atómica en horno de grafito a 283.3 nm283.3\,\mathrm{nm} (datos de ejercicio, usados en el problema). Izquierda: los cinco patrones y la recta de mínimos cuadrados. Derecha: los residuos, pequeños y sin tendencia: una recta es adecuada.

Método 34.12 (Calibrar)

  1. Prepárense cinco o más patrones que abarquen el intervalo esperado, en la misma matriz que las muestras, y un blanco.
  2. Ajústese la recta de mínimos cuadrados; represéntense los residuos. Una tendencia curva significa que la respuesta no es lineal: estréchese el intervalo o ajústese una curva. Un residuo grande aislado señala un patrón defectuoso.
  3. Léanse las muestras, preferiblemente cerca del centro del intervalo, por replicado; calcúlense x0x_0 y sx0s_{x_0}.
  4. Exprésese x0±t sx0x_0 \pm t\,s_{x_0} tras corregir cada dilución, añadiendo la incertidumbre propagada de la dilución si no es despreciable.

Definición 34.13 (Adiciones estándar)

En el método de las adiciones estándar, se añaden cantidades conocidas de analito a porciones iguales de la propia muestra; se representa la señal frente a la concentración añadida, y la concentración del analito en la muestra es la distancia del origen al punto en que la recta extrapolada a señal nula corta el eje, a/ba/b.

Método 34.14 (Adiciones estándar)

  1. Úsese cuando la matriz de la muestra cambia la respuesta (sales, ácidos, materia orgánica), de modo que patrones en agua pura darían una pendiente errónea.
  2. Añádanse a porciones iguales cantidades conocidas crecientes, llévense todas al mismo volumen y mídanse.
  3. Ajústese la recta; la concentración en las disoluciones medidas es a/ba/b; corríjase la dilución.
  4. La respuesta debe ser lineal desde cero: compruébense los residuos.
Adiciones estándar (datos de ejercicio): absorbancia de una muestra a la que se añaden 0, 2, 4 y 6\, mg/ L de analito. La recta corta la señal nula en -3.11\, mg/ L: la disolución medida contenía 3.11\, mg/ L.
Adiciones estándar (datos de ejercicio): absorbancia de una muestra a la que se añaden 0, 2, 4 y 6 mg/L6\,\mathrm{mg}/\mathrm{L} de analito. La recta corta la señal nula en −3.11 mg/L-3.11\,\mathrm{mg}/\mathrm{L}: la disolución medida contenía 3.11 mg/L3.11\,\mathrm{mg}/\mathrm{L}.

34.5 Detección, cuantificación y validación

Definición 34.15 (Límites)

El límite de detección de un método es la menor concentración cuya señal puede distinguirse de la de un blanco con un riesgo pequeño y declarado de falsa detección; el límite de cuantificación es la menor concentración que puede medirse con una incertidumbre relativa aceptable.

Proposición 34.16 (Estimaciones habituales)

Con sblanks_{\mathrm{blank}} la desviación típica de señales repetidas del blanco y bb la pendiente de la recta de calibrado, xLOD≈3sblank/bx_{\mathrm{LOD}} \approx 3s_{\mathrm{blank}}/b y xLOQ≈10sblank/bx_{\mathrm{LOQ}} \approx 10s_{\mathrm{blank}}/b.

Argumento. Si las señales del blanco son normales con desviación típica sblanks_{\mathrm{blank}}, un blanco supera su media en más de 3sblank3s_{\mathrm{blank}} con una probabilidad de alrededor del 0.13 % (cola unilateral de la ley normal más allá de 3 desviaciones típicas): considerar “detectada” una señal por encima de ese umbral rara vez confunde un blanco con una muestra. Una señal 10sblank10s_{\mathrm{blank}} por encima del blanco tiene una desviación típica relativa de alrededor del 10 %, el umbral habitual para un resultado cuantitativo. Dividir por bb convierte las señales en concentraciones. ∎

Definición 34.17 (Validación)

La veracidad es la proximidad entre la media de un gran número de resultados y el valor verdadero; su diferencia es el sesgo. La precisión es la proximidad entre resultados independientes, expresada como una desviación típica: repetibilidad en las mismas condiciones (mismo operador, instrumento, laboratorio, poco tiempo), reproducibilidad en condiciones distintas (laboratorios diferentes). La validación de un método es el conjunto de experimentos que establecen, para un método, su veracidad, su precisión, su linealidad e intervalo, y sus límites de detección y de cuantificación.

La veracidad se comprueba frente a un material de referencia certificado o por adición de analito; la precisión, con réplicas en un laboratorio y con comparaciones entre laboratorios. Un método puede ser preciso pero sesgado, con cada resultado próximo a los demás y todos ellos erróneos: la precisión no prueba la veracidad.

Historia — Student, 1908

William Sealy Gosset, un químico que trabajaba para una cervecera, tenía que juzgar materias primas a partir de muy pocas muestras, demasiado pocas para la ley normal. En 1908 publicó la distribución de la media dividida por la desviación típica muestral, con el seudónimo “Student” porque su empresa no permitía al personal publicar con su nombre. La ley de Student se usa cada vez que un analista da un intervalo a partir de tres o cuatro réplicas. (Fotografía, 1908: dominio público; Wikimedia Commons.)

34.6 Ejercicios

Ejercicio 34.1 ★

Cinco valoraciones dan 12.45, 12.50, 12.48, 12.52 y 12.46 mL12.46\,\mathrm{mL}. Calcúlense la media, la desviación típica experimental y la desviación típica de la media.

Solución

Solución de Ejercicio 34.1.

Vˉ=62.41/5=12.482 mL\bar V = 62.41/5 = 12.482\,\mathrm{mL}; ∑(Vi−Vˉ)2=0.00328\sum(V_i - \bar V)^2 = 0.00328, s=0.00328/4≈0.029 mLs = \sqrt{0.00328/4} \approx 0.029\,\mathrm{mL}; s/5≈0.013 mLs/\sqrt5 \approx 0.013\,\mathrm{mL}.

Ejercicio 34.2 ★

Dese el intervalo de confianza al 95 % de la media del ejercicio 1.

Solución

Solución de Ejercicio 34.2.

t=2.776t = 2.776 para 4 grados de libertad: 12.482±2.776×0.012812.482 \pm 2.776 \times 0.0128, es decir, (12.48±0.04)(12.48 \pm 0.04) mL.

Ejercicio 34.3 ★

Se prepara una disolución disolviendo m=0.5000 gm = 0.5000\,\mathrm{g} (u=0.0002 gu = 0.0002\,\mathrm{g}) de un sólido en un matraz de V=250.0 mLV = 250.0\,\mathrm{mL} (u=0.15 mLu = 0.15\,\mathrm{mL}). Calcúlese la incertidumbre típica relativa de c=m/(MV)c = m/(MV), siendo la masa molar suficientemente exacta.

Solución

Solución de Ejercicio 34.3.

Las incertidumbres relativas de la masa y del volumen son del 0.04 % y del 0.06 %; combinadas cuadráticamente,

u(c)c=(0.00020.5000)2+(0.15250.0)2≈0.07 %.\frac{u(c)}{c} = \sqrt{\Bigl(\frac{0.0002}{0.5000}\Bigr)^2 + \Bigl(\frac{0.15}{250.0}\Bigr)^2} \approx 0.07~\%.

Ejercicio 34.4 ★

Ajústese la recta de mínimos cuadrados que pasa por (0,0.02)(0, 0.02), (1,0.21)(1, 0.21), (2,0.39)(2, 0.39), (3,0.62)(3, 0.62).

Solución

Solución de Ejercicio 34.4.

xˉ=1.5\bar x = 1.5, yˉ=0.31\bar y = 0.31, Sxx=5S_{xx} = 5, Sxy=0.99S_{xy} = 0.99: b=0.198b = 0.198, a=0.31−0.198×1.5=0.013a = 0.31 - 0.198 \times 1.5 = 0.013.

Ejercicio 34.5 ★★

Una disolución de referencia certificada contiene 10.00 mg/L10.00\,\mathrm{mg}/\mathrm{L} de un elemento. Cinco análisis dan una media de 10.12 mg/L10.12\,\mathrm{mg}/\mathrm{L} con s=0.08 mg/Ls = 0.08\,\mathrm{mg}/\mathrm{L}. ¿Está sesgado el método al nivel del 95 %?

Solución

Solución de Ejercicio 34.5.

∣10.12−10.00∣/(0.08/5)=0.12/0.036≈3.4>2.776|10.12 - 10.00|/(0.08/\sqrt5) = 0.12/0.036 \approx 3.4 > 2.776: la diferencia es significativa, y el método está sesgado (en unos +0.12 mg/L+0.12\,\mathrm{mg}/\mathrm{L}).

Ejercicio 34.6 ★★

Una concentración de iones hidronio se conoce con una incertidumbre típica relativa del 2 %. ¿Cuál es la incertidumbre típica del pH?

Solución

Solución de Ejercicio 34.6.

pH=−log⁡10c\mathrm{pH} = -\log_{10}c, d pH/dc=−1/(cln⁡10)\mathrm d\,\mathrm{pH}/\mathrm dc = -1/(c\ln10): u(pH)=(u(c)/c)/ln⁡10=0.02/2.303≈0.009u(\mathrm{pH}) = (u(c)/c)/\ln10 = 0.02/2.303 \approx 0.009.

Ejercicio 34.7 ★★

Diez blancos tienen s=0.0012s = 0.0012 unidades de absorbancia; la pendiente de calibrado es 0.0098 L/µg0.0098\,\mathrm{L}/\text{µ}\mathrm{g}. Estímense los límites de detección y de cuantificación.

Solución

Solución de Ejercicio 34.7.

xLOD=3×0.0012/0.0098≈0.37 µg/Lx_{\mathrm{LOD}} = 3 \times 0.0012/0.0098 \approx 0.37\,\text{µ}\mathrm{g}/\mathrm{L}; xLOQ=10×0.0012/0.0098≈1.2 µg/Lx_{\mathrm{LOQ}} = 10 \times 0.0012/0.0098 \approx 1.2\,\text{µ}\mathrm{g}/\mathrm{L}.

Ejercicio 34.8 ★★

En el experimento de adiciones estándar de la figura, las cuatro porciones se habían preparado diluyendo 10.0 mL10.0\,\mathrm{mL} de muestra hasta 50.0 mL50.0\,\mathrm{mL}. Dese la concentración en la muestra.

Solución

Solución de Ejercicio 34.8.

La recta da 3.11 mg/L3.11\,\mathrm{mg}/\mathrm{L} en las disoluciones medidas; la muestra se había diluido cinco veces: 3.11×50.0/10.0≈15.6 mg/L3.11 \times 50.0/10.0 \approx 15.6\,\mathrm{mg}/\mathrm{L}.

Ejercicio 34.9 ★★

En un estudio interlaboratorio, cada laboratorio obtiene s=0.05 mg/Ls = 0.05\,\mathrm{mg}/\mathrm{L} en sus propias réplicas, pero los resultados de los distintos laboratorios se dispersan con s=0.15 mg/Ls = 0.15\,\mathrm{mg}/\mathrm{L}. Nómbrense las dos magnitudes y explíquese la diferencia.

Solución

Solución de Ejercicio 34.9.

0.05 mg/L0.05\,\mathrm{mg}/\mathrm{L} es la repetibilidad, y 0.15 mg/L0.15\,\mathrm{mg}/\mathrm{L}, la reproducibilidad. Cada laboratorio tiene su propio pequeño sesgo (calibrado, patrones, instrumento); estos sesgos difieren de un laboratorio a otro y se suman a la dispersión, de modo que la reproducibilidad es siempre la mayor.

Ejercicio 34.10 ★★★

Demuéstrese que la recta de mínimos cuadrados pasa por (xˉ,yˉ)(\bar x, \bar y) y que sus residuos suman cero.

Solución

Solución de Ejercicio 34.10.

∂S/∂a=0\partial S/\partial a = 0 se escribe ∑(yi−a−bxi)=0\sum(y_i - a - bx_i) = 0: los residuos suman cero. Dividiendo por nn: yˉ=a+bxˉ\bar y = a + b\bar x, de modo que (xˉ,yˉ)(\bar x, \bar y) está sobre la recta.

Ejercicio 34.11 ★★★

Una disolución madre se diluye cien veces, bien en una sola etapa (1.000 mL1.000\,\mathrm{mL} en 100.0 mL100.0\,\mathrm{mL}), bien en dos etapas de diez veces (10.00 mL10.00\,\mathrm{mL} en 100.0 mL100.0\,\mathrm{mL}, dos veces). Incertidumbres típicas: pipeta de 1 mL 0.006 mL0.006\,\mathrm{mL}, pipeta de 10 mL 0.02 mL0.02\,\mathrm{mL}, matraz de 100 mL 0.08 mL0.08\,\mathrm{mL} (datos de ejercicio). ¿Qué vía es más precisa?

Solución

Solución de Ejercicio 34.11.

En una etapa:

(0.006/1.000)2+(0.08/100.0)2≈0.61 %.\sqrt{(0.006/1.000)^2 + (0.08/100.0)^2} \approx 0.61~\%.

Una etapa de diez veces:

(0.02/10.00)2+(0.08/100.0)2≈0.22 %,\sqrt{(0.02/10.00)^2 + (0.08/100.0)^2} \approx 0.22~\%,

y dos de ellas combinadas cuadráticamente dan 0.222≈0.30 %0.22\sqrt2 \approx 0.30~\%. La vía en dos etapas es el doble de precisa: la pipeta pequeña domina la vía en una etapa.

Ejercicio 34.12 ★★★

Los dos laboratorios de la introducción informan de 9.6 y 11.2 µg/L11.2\,\text{µ}\mathrm{g}/\mathrm{L} con incertidumbres típicas 0.4 y 0.5 µg/L0.5\,\text{µ}\mathrm{g}/\mathrm{L} (datos de ejercicio). ¿Son compatibles sus resultados? ¿Qué dice cada uno sobre el valor guía?

Solución

Solución de Ejercicio 34.12.

z=1.6/0.42+0.52≈2.5>2z = 1.6/\sqrt{0.4^2 + 0.5^2} \approx 2.5 > 2: no son compatibles; uno de los laboratorios (al menos) tiene un sesgo. Con incertidumbres expandidas (k=2k = 2), el primero da 9.6±0.89.6 \pm 0.8, que incluye 10; el segundo, 11.2±1.011.2 \pm 1.0, enteramente por encima de 10. Antes de cualquier veredicto sobre el agua hay que resolver la discrepancia, por ejemplo con un material de referencia.

34.7 Problema: plomo en el agua del grifo

Problema 34.1

Problema de fin de semana — un calibrado de absorción atómica por mínimos cuadrados, la concentración de una muestra con su intervalo de confianza, la dilución y los límites, y un veredicto frente al valor guía

Datos de ejercicio. El plomo se mide por absorción atómica en horno de grafito a 283.3 nm283.3\,\mathrm{nm}. Patrones (µg/L\text{µ}\mathrm{g}/\mathrm{L}; absorbancia): 0.0; 0.003, 5.0; 0.051, 10.0; 0.101, 15.0; 0.148, 20.0; 0.199. El agua del grifo (50.0 mL50.0\,\mathrm{mL}) se acidula con 0.50 mL0.50\,\mathrm{mL} de ácido nítrico; tres lecturas de esta disolución dan 0.104, 0.098 y 0.101. Diez blancos: 0.002, 0.004, 0.003, 0.001, 0.003, 0.005, 0.002, 0.003, 0.004, 0.003. Incertidumbres típicas de los volúmenes: 0.03 mL0.03\,\mathrm{mL} para 50.0 mL, 0.005 mL0.005\,\mathrm{mL} para 0.50 mL. Valor guía: 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L}.

Parte I — Calibrado.

  1. ¿Por qué se preparan los patrones en el mismo ácido nítrico que las muestras?
  2. Calcúlense la pendiente y la ordenada en el origen de la recta de mínimos cuadrados.
  3. Calcúlense los cinco residuos.
  4. ¿Muestran los residuos alguna tendencia? Conclúyase sobre la linealidad.
  5. Calcúlese sy/xs_{y/x}.
  6. ¿Por qué la ordenada en el origen no es nula?
  7. ¿Qué mide la pendiente?

Parte II — La muestra.

  1. Calcúlense la media y la desviación típica de las tres lecturas.
  2. Calcúlese la concentración de plomo de la disolución medida.
  3. Calcúlese sx0s_{x_0}.
  4. ¿Cuántos grados de libertad, y qué coeficiente de Student?
  5. Dese el intervalo al 95 % de la concentración en la disolución medida.
  6. ¿Por qué leer la muestra tres veces y no una?

Parte III — Dilución y límites.

  1. Calcúlense el factor de dilución y la concentración en el agua del grifo.
  2. Propáguese la incertidumbre de los dos volúmenes al factor de dilución.
  3. ¿Es significativa esa incertidumbre frente a la del calibrado?
  4. Calcúlense la desviación típica de los blancos y el límite de detección.
  5. Calcúlese el límite de cuantificación.
  6. ¿Está la muestra por encima del límite de cuantificación?

Parte IV — Veredicto.

  1. ¿Contiene el intervalo el valor guía?
  2. ¿Qué habría sugerido una sola lectura de 0.104?
  3. ¿Qué significa aquí “confianza del 95 %”?
  4. ¿Cómo podría el laboratorio decidir con más firmeza?
  5. ¿Cómo comprobaría que el método no está sesgado?
  6. El valor guía es provisional “sobre la base del rendimiento del tratamiento y de la factibilidad analítica”. ¿Qué significa la segunda razón, a la luz de este problema?
  7. Indíquese la concentración de plomo del agua del grifo con su intervalo al 95 %, y compárese con 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L}.
Solución

Solución de Problema 34.1.

1. El ácido modifica la atomización y la respuesta; los patrones y las muestras deben tener la misma matriz para que la pendiente sea aplicable. 2. xˉ=10\bar x = 10, yˉ=0.1004\bar y = 0.1004, Sxx=250S_{xx} = 250, Sxy=2.445S_{xy} = 2.445: b=0.00978b = 0.00978 por µg/L\text{µ}\mathrm{g}/\mathrm{L}, a=0.1004−0.0978=0.0026a = 0.1004 - 0.0978 = 0.0026. 3. +0.0004+0.0004, −0.0005-0.0005, +0.0006+0.0006, −0.0013-0.0013, +0.0008+0.0008. 4. Los signos alternan, sin curvatura: la recta es adecuada en el intervalo 0–20 µg/L20\,\text{µ}\mathrm{g}/\mathrm{L}. 5. sy/x=3.1×10−6/3≈0.00102s_{y/x} = \sqrt{3.1 \times 10^{-6}/3} \approx 0.00102. 6. El blanco (ácido, agua, horno) da por sí mismo una pequeña absorbancia. 7. La sensibilidad: absorbancia por µg/L\text{µ}\mathrm{g}/\mathrm{L} de plomo. 8. yˉ0=0.1010\bar y_0 = 0.1010; s=0.0030s = 0.0030. 9. x0=(0.1010−0.0026)/0.00978≈10.06 µg/Lx_0 = (0.1010 - 0.0026)/0.00978 \approx 10.06\,\text{µ}\mathrm{g}/\mathrm{L}. 10. sx0=(0.00102/0.00978)1/3+1/5+(0.1010−0.1004)2/(0.009782×250)≈0.104×0.730≈0.076 µg/Ls_{x_0} = (0.00102/0.00978)\sqrt{1/3 + 1/5 + (0.1010 - 0.1004)^2/(0.00978^2 \times 250)} \approx 0.104 \times 0.730 \approx 0.076\,\text{µ}\mathrm{g}/\mathrm{L}. 11. n−2=3n - 2 = 3; t=3.182t = 3.182. 12. 10.06±3.182×0.076=(10.06±0.24)10.06 \pm 3.182 \times 0.076 = (10.06 \pm 0.24) µg/L\text{µ}\mathrm{g}/\mathrm{L}. 13. El término 1/m1/m baja de 1 a 1/31/3: sx0s_{x_0} disminuye en alrededor de un tercio. 14. f=50.50/50.0=1.010f = 50.50/50.0 = 1.010; 10.06×1.010≈10.16 µg/L10.06 \times 1.010 \approx 10.16\,\text{µ}\mathrm{g}/\mathrm{L}. 15. f=1+V2/V1f = 1 + V_2/V_1: u2(f)=(u(V2)/V1)2+(V2u(V1)/V12)2=(0.005/50.0)2+(0.50×0.03/2500)2u^2(f) = (u(V_2)/V_1)^2 + (V_2u(V_1)/V_1^2)^2 = (0.005/50.0)^2 + (0.50 \times 0.03/2500)^2, u(f)≈1.0×10−4u(f) \approx 1.0 \times 10^{-4}, un 0.01 % relativo. 16. No: 0.01 % frente a 0.076/10.06≈0.75 %0.076/10.06 \approx 0.75~\%. 17. sblank≈0.00115s_{\mathrm{blank}} \approx 0.00115; xLOD=3×0.00115/0.00978≈0.35 µg/Lx_{\mathrm{LOD}} = 3 \times 0.00115/0.00978 \approx 0.35\,\text{µ}\mathrm{g}/\mathrm{L}. 18. 10×0.00115/0.00978≈1.2 µg/L10 \times 0.00115/0.00978 \approx 1.2\,\text{µ}\mathrm{g}/\mathrm{L}. 19. Sí, con mucho. 20. Sí: de 9.92 a 10.40 µg/L10.40\,\text{µ}\mathrm{g}/\mathrm{L} en el agua del grifo. 21. (0.104−0.0026)/0.00978×1.010≈10.5 µg/L(0.104 - 0.0026)/0.00978 \times 1.010 \approx 10.5\,\text{µ}\mathrm{g}/\mathrm{L} y, leída sola, un “por encima del valor guía” categórico que los datos no respaldan. 22. La regla usada para construir el intervalo captura el valor verdadero en el 95 % de las series a las que se aplica; no es una probabilidad sobre este valor concreto. 23. Más lecturas replicadas y más patrones cerca de 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L} (los términos 1/m1/m y 1/n1/n), o un método más sensible; el intervalo solo se estrecha lentamente, como 1/nuˊmero1/\sqrt{\text{número}}. 24. Analizando un agua de referencia certificada, o añadiendo a la muestra una cantidad conocida de plomo y comprobando la recuperación. 25. Cerca de 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L} la incertidumbre de los métodos de rutina es de unos pocos por ciento, como aquí: un valor guía más bajo no podría ser comprobado de forma fiable por los laboratorios ordinarios. 26. (10.16±0.24)\boldsymbol{(10.16 \pm 0.24)} µg/L\text{µ}\mathrm{g}/\mathrm{L} al 95 %: el intervalo contiene 10 µg/L10\,\text{µ}\mathrm{g}/\mathrm{L}, de modo que la medida no demuestra ni que el agua supere el valor guía ni que lo cumpla con margen.

Términos definidos en este capítulo

Ver los 852 términos del glosario