1_TomaMuestra
1_TomaMuestra
1_TomaMuestra
1. Introducción
Parece que existe consenso en clasificar una muestra dentro de uno de los dos
grandes grupos siguientes: lote de muestra a granel y lote de muestra manufacturada.
Luego vendría la cadena de transformación del lote de muestra hasta la porción que se
analiza (figura 1). Esta cadena de transformación depende de la naturaleza de la muestra
problema inicial, de los analitos y del proceso analítico posterior. Parte de esta cadena de
transformación tendrá lugar en la zona de toma de muestra (fuera del laboratorio) y en el
laboratorio de análisis.
1
http://www.rand.org/pubs/monograph_reports/MR1418/
2
https://chem.libretexts.org/Bookshelves/Analytical_Chemistry/Book%3A_Analytical_Chemistry_2.0_(Ha
r vey)/07%3A_Collecting_and_Preparing_Samples/7.2%3A_Designing_a_Sampling_Plan
00 01 02 98 99
00
01
02
98
99
Punto caliente
4.1 Caso 1
Supongamos que la variabilidad más importante de las etapas totales de un
análisis se encuentra en la etapa de toma de muestra (lo habitual). Supongamos, además,
que nuestra población objetivo está constituida por dos tipos de partículas A y B. Las
partículas de tipo A contienen una concentración determinada en el analito de interés,
Canalito,A , y en las partículas B el analito está ausente. Las dimensiones de los dos tipos de
partículas suelen ser las mismas ya que, normalmente el lote proviene de un
triturado/cribado previo que afecta de las misma forma a todas ellas.
Para este caso, dependiendo de cuántas partículas de tipo A recogiéramos al hacer
el muestreo así sería la concentración media que diríamos que tiene la población. Si,
hipotéticamente cogemos una muestra de 1 partícula y ésta es de tipo A, estimaríamos la
concentración promedio del analito en la población como Canalito,A; si, por el contrario,
hipotéticamente, cogiéramos en el muestreo una sola partícula y esta fuera de tipo B,
estimaríamos que la concentración del analito en la población sería nula. En general,
n A C analito, A + n B 0
C analito = = p A C analito, A
n
donde n es el número de partículas tomadas, nA y nB es el número de partículas tomadas
de A y de B, respectivamente, pA es la probabilidad de tomar una partícula de tipo A y
Canalito es la concentración media del analito estimada. Por lo tanto, la desviación estándar
relativa al cuadrado (Sr2) en Canalito será la suma de la desviación estándar relativa al
cuadrado asociada a pA y la desviación estándar relativa al cuadrado asociada a la
concentración intrínseca del analito en las partículas de tipo A.
S r2,Canalito = S r2, p A + S r2,Canalito, A S r2, p A
En la mayoría de los casos la variabilidad en la concentración existente del analito
en las partículas de tipo A (y su determinación analítica) suele ser mucho menor que la
asociada al conocimiento de la proporción de partículas de A existentes en la muestra,
con lo que la variabilidad final en la concentración media estimada para el analito sólo
dependerá de la precisión con la que podamos conocer la proporción de partículas de tipo
A que existen en la población, la cual corresponde a la varianza (S2) de una distribución
poblacional binomial:
S p2A = n p A (1 − p A )
Por lo tanto,
n p A (1 − p A )
S r2,Canalito S r2, p A =
(n p A ) 2
Despejando n
(1 − p A ) 1
n=
p A S r2, p A
Por lo tanto, para alcanzar un grado de precisión prefijado (Sr2) será necesario
tomar n partículas. Lógicamente, para unas condiciones fijas de densidad y volumen de
cada partícula podrá establecerse una correlación directa entre el número de partículas y
la masa de muestra para tomar.
Ejemplo 0: Sea un suelo donde las partículas conteniendo al analito representan
tan solo 1x10-7% de la población. ¿Cuántas partículas habría que tomar para alcanzar
una desviación estándar relativa en el muestreo del ±1%? Considérese que la fuente de
variabilidad principal está en la toma de muestra.
10−7
(1 − )
100 1
n= = 1013 partículas
10 −7
( 0.01) 2
100
Supongamos que ese número de partículas corresponde a una masa de 80 g.
Como esa cantidad de masa es muy grande y no es práctico manejar muestras de esa
envergadura para nuestros propósitos ¿Cómo podríamos hacer para trabajar con
muestras de 0.60 g, sin perder precisión en la toma de muestra (mantener la misma
desviación estándar relativa)?
Para no perder precisión necesitamos pulverizar más las partículas de forma que
podamos tomar ese mismo número de partículas (1013) tomando menos cantidad de masa
(0.60g). Si suponemos esféricas las partículas de radio r y del mismo tamaño las de tipo
A que las del resto de tipos (por esos procesos de triturado/cribado),
4
masa = volumen n r 3
3
masa n r 3
1/ 3
0.60g 1
rparticulas finales = rparticulas finiciales = rparticulas finiciales
80g 5.1
Es decir, haciendo las partículas 5 veces más pequeñas, podríamos tomar tan sólo
muestras de 0.60g y ya nos llevaríamos 1013 partículas, con lo que la desviación estándar
relativa en la toma de muestra se mantendría en el ±1%.
4.2 Caso 2
Como se ha visto en el ejemplo, el caso 1 (distribución binomial) también es
aplicable a poblaciones que tengan más de dos tipos partículas, siempre y cuando el
analito esté presente en un solo tipo de ellas, quedando el resto de los tipos de partículas,
como “partícula B” de la distribución binomial.
Ahora bien, una población real puede contener más de dos partículas estando el
analito presente (a distintos niveles de concentración) en cada uno de esos tipos de
partículas y la expresión de la distribución binomial ya no sería aplicable. En muchos de
estos casos, sin embargo, si suponemos que la muestra heterogénea de laboratorio está
muy bien mezclada, de manera que aunque sea heterogénea (partículas de diferente
tamaño y tipo) las partículas puedan considerarse que están aleatoriamente distribuídas en
la muestra y aleatoriamente tomadas durante la recogida de sub-muestras: es decir que
una muestra combinada formada a partir de dos muestras simples (cada una de peso m)
tenga la misma variabilidad que una única muestra simple de peso 2m. En estos casos,
que sería como tratar con muestras heterogéneas que (a la escala a la que se lleva a cabo
el muestreo) se comportan como homogéneas3 se cumple la relación4
1
m = Ks 2
S r ,Canalito 1002
donde m es la masa de la muestra tomada y Sr2 es la desviación estándar relativa en la
Canalito al cuadrado. Esta expresión nos indica que se establece una relación inversa entre
la masa tomada de una muestra de este tipo “uniforme” y la varianza relativa en la toma
de muestra. A la constante de proporcionalidad, Ks, se le denomina constante de
Ingamells, la cual hay que calcular experimentalmente para cada unas condiciones
concretas de una muestra de este tipo. Ks correspondería, por tanto, a aquella masa de
muestra que produce una desviación estándar relativa de la toma de muestra del ±1%.
Eso sí, igual que en el caso anterior, se considera que la fuente de variabilidad principal
está en la toma de muestra.
Ejemplo 1: En la determinación del contenido residual de ceniza inorgánica en
cereales para el desayuno se obtuvieron los siguientes resultados preliminares:
Masa cereal (g) 0.9956 0.9981 1.0036 0.9994 1.0067
% peso ceniza 1.34 1.29 1.32 1.26 1.28
Estimar el valor de Ks y determinar qué tamaño de muestra de cereal es
necesario para alcanzar una desviación estándar relativa en el muestreo del ±2.0%.
Considérese que la fuente de variabilidad principal está en la toma de muestra.
Para esos análisis que han sido realizados a un mismo nivel de masa de muestra
podemos suponer que su varianza es la misma para las cinco muestras analizadas y
calcular la desviación estándar absoluta y relativa:
1.34 + 1.29 + 1.32 + 1.26 + 1.28
C analito = = 1.298%( p / p)
5
S Canalito = 0.0319%( p / p)
S r ,Canalito = 0.0246
Ese valor de desviación estándar relativa se da para un nivel de masa de muestra
de cereal de 1.000g de cereal (valor medio). Por lo tanto, la constante Ks estimada para
este tipo de muestras será:
Ks = m ( S r2,Canalito 100) 2 = 1.000g (0.0246) 2 1002 =
= 6.06g
Luego la masa de muestra de cereal necesaria, como mínimo, para alcanzar una
precisión del ±1% (o inferior) es de 6 gramos o superior. Si, en nuestro estudio, nos
conformamos con una precisión en el muestreo del ±2.0%, el tamaño que tendríamos que
coger de muestra de cereal sería, como mínimo, de:
3
Es decir, en las porciones de masa tomadas que no sean muy pequeñas para que no se manifieste dicha
heterogeneidad.
4
C.O.Ingamells, Talanta 1973 (20) 547-568.
1 1
m = Ks = 6.06g = 1.51g de cereal
S 2
r ,Canalito 100 2
(0.020) 2 1002
1/ 3
1.00g 1
rparticulas finales = rparticulas finiciales = rparticulas finiciales
3.462g 3.5
queda que el número de muestras para tomar a fin de alcanzar un error indeterminado
concreto sería:
t 2 S sampl
2
n sampl =
,Canalito
C2 analito
n rep =
,Canalito
C2 analito
0.40 0.070
r ,C (estrategia 2/5) = 2.262
+ = 1.00%
analito
2 2x5
Por lo tanto, la estrategia primera resulta más adecuada ya que se minimiza la
parte del método que mayor variabilidad posee, que, en este caso es la toma de muestra.
El nrep lo vamos variando de forma que vamos obteniendo distintos valores de r
(tanteo) hasta que el error calculado sea igual o inferior al deseado (la diferencia se hace
cero o negativa). El valor de la t-Student es para el 0.05 de significación y nsampl*nrep-1
grados de libertad.
Vemos que con una sola muestra tomada no es posible, aun haciendo multitud de
réplicas de la misma llegar a igualar o mejorar el error deseado. A partir de la toma de
dos muestras, haciendo las réplicas necesarias conseguimos igualar o mejorar el error
deseado. Las estrategias (nsampl/nrep) 2/12; 3/4; 4/3 y 5/2 cumplen el requisito de coste.
Eso sí, la 5/2 es la estrategia donde se sitúa, en este caso, el mínimo del coste.
Con los datos de la primera estrategia podemos estimar una varianza global de la
determinación que incluye tanto la varianza del método como la de la toma de muestra:
S glob,Canalito = 0.0750 ppm S glob
2
,Canalito = 0.005625(ppm)
2
,Canalito = S sampl ,Canalito + S meth ,Canalito S sampl ,Canalito = S glob,Canalito − S meth ,Canalito
2 2 2 2 2 2
S glob
0.00460 0.001030
C (estrategia 2/3) = 2.571 + = 0.127 ppm
analito
2 2 x3
O, en forma porcentual:
0.093 ppm
r ,C (%)(estrategia 5/1) = x100 = 8.86%
analito
1.05 ppm
0.127 ppm
r ,C (estrategia 2/3) = = 12.1%
analito
1.05 ppm
que también puede calcularse, obviamente, a partir de la expresión inicial del error
indeterminado a partir de los valores de las desviaciones estándar relativas, a ese nivel
promedio de concentración de 1.05 ppm:
5 5 5x1
minicial ( S r2, sampl ,Canalito ) inicial (100) 2 = m final ( S r2, sampl ,Canalito ) final (100) 2
(2)
(3)
Figura 6. Sondas para la toma de muestras particuladas sólidas estáticas. (1) y (2) sondas
acanaladas; (3) sonda con orificios y con tubo interno.
En las sondas acanaladas (figura 6.1 y figura 6.2) la sonda se introduce, se rota y
se retira el contenido con el canal colocado hacia la parte superior. La sonda con tubo
interior permite cerrar los orificios del tubo externo hasta que se haya introducido la
sonda. Al girar el tubo interior estos orificios se abren y entran las partículas. Luego se
gira de nuevo para cerrar los orificios y así aislar la muestra tomada, antes de sacar de
nuevo la sonda. Las sondas suelen ser de unos 40-80 cm de largo y de 15-50 mm de
diámetro interno. Además han de ser de un material duro e inerte que no contamine el
interior de la muestra porque se fragmente al ser introducida la sonda.
Para que funcione bien este tipo de toma de muestra es necesario que el material
sea fino para que pueda rellenar sin problemas la sonda.
Figura 7. Sonda con dispositivo perforador para toma de muestras sólidas compactas.
Para muestras que no sea posible taladrar, por ejemplo un bloque de metal
precioso, si el material es suficientemente homogéneo basta con tomar una porción de un
extremo o de la superficie; cuando son, por ejemplo obras de arte y tampoco es posible
taladrar, sino lo que se necesita es una mínima invasión del objeto, se emplean
microsondas.
7.1 Criterios sanitarios de la calidad del agua de consumo humano (BOE 21-02-
2003)
Ejemplo:
Consideremos una línea de producción que genera lotes de 1000 unidades sobre el
que deseamos realizar un muestreo simple de 80 unidades del lote con un número de
aceptación de 4, es decir, que se aceptará el lote si de las 80 muestras tomadas hay como
mucho 4 defectuosas. Si se sabe que la línea de producción en cuestión ha presentado en
lotes de similares características un 6% de unidades defectuosas, nos piden estimar la
probabilidad que tenemos de que un lote “tipo” de estas características sea aceptado en la
inspección.
N= 1000
p(promedio)=pm=0.06
(n,c)= (80,4)
Se cumplen las condiciones para utilizar la Distribución de Poisson. Utilizando,
por ejemplo, el excel, obtenemos la probabilidad acumulada HASTA el valor de “x”=4 y
para un valor promedio de la distribución de n*p=80*0.06=4.8.
PAceptación=POISSON(4;4.8;VERDADERO)=0.476
Es decir, se estima que la probabilidad de aceptación de un lote “tipo” con ese
plan de muestreo es del 47.6%. Si yo quiero que la probabilidad de aceptación del lote
aumente, pero con la misma tasa de aceptación (4/80), ¿qué pasaría si incremento el
tamaño de la muestra tomada o lo disminuyo?
(n,c)=(60,3) PAceptación=POISSON(3;3.6;VERDADERO)=0.515
(n,c)=(100,5) PAceptación=POISSON(5;6;VERDADERO)=0.446
Vemos que tampoco cambia mucho ya que le pedimos a la muestra tomada la
misma tasa de defectuosos. Sólo cuando bajamos mucho el tamaño de la muestra se
pueden encontrar diferencias
(n,c)=(20,1) PAceptación=POISSON(1;1.2;VERDADERO)=0.663
Eso sí, se incrementa la probabilidad de aceptación del lote para muestras de 20
elementos, respecto a una toma de muestra de 80 elementos (para una misma tasa de
defectuosos (4/80=1/20)), pero eso no quiere decir que el cliente quede más satisfecho o
que resulte más barato el coste del plan de muestreo, como más adelante veremos.
2
MUESTREO
0
3
1
4 0
PBinom(0)=DISTR.BINOM(0;10;0.02;FALSO)=0.8171
PBinom(1)=DISTR.BINOM(1;10;0.02;FALSO)=0.1667
Por lo tanto:
PAceptación=0.3679+0.3679+0.1839+0.0613*0.8171+0.0613*0.1667+
+0.0153*0.8171=0.9925
Se concluye que la probabilidad de aceptación del lote “tipo” es de un 99.25%
con ese plan de muestreo doble así diseñado.
Nota que podríamos haber calculado la probabilidad de aceptación del muestreo
simple anterior de forma similar, sumando probabilidades puntuales en lugar de la
acumulada:
0
2
MUESTREO
3
PPoisson(0)=POISSON(0;4.8;FALSO)=0.0082
PPoisson(1)=POISSON(1;4.8;FALSO)=0.0395
PPoisson(2)=POISSON(2;4.8;FALSO)=0.0948
PPoisson(3)=POISSON(3;4.8;FALSO)=0.1517
PPoisson(4)=POISSON(4;4.8;FALSO)=0.1820
PAceptación=0.0082+0.0395+0.0948+0.1517+0.1820=0.4762
MUESTREO 0
2
1 0
3 0 0
PBinom(0)=DISTR.BINOM(0;7;0.04;FALSO)=0.7514
PBinom(1)=DISTR.BINOM(1;7;0.04;FALSO)=0.2192
PBinom(2)=DISTR.BINOM(2;7;0.04;FALSO)=0.027396
PBinom(3)=DISTR.BINOM(3;7;0.04;FALSO)=0.001925
PPoisson(0)=POISSON(0;50*0.04;FALSO)=0.1353
PPoisson(1)=POISSON(1;50*0.04;FALSO)=0.27067
P’Poisson(0)=POISSON(0;80*0.04;FALSO)=0.04076
PAceptación=0.7514+0.2192+0.027396*0.1353+0.027396*0.27067*0.04076+
+0.001925*0.1353*0.04076=0.9746
Es decir, con ese plan de muestreo diseñado para ese sistema de producción (o lo
que sea) funcionando, existe una probabilidad del 97.46% de que un lote sea aceptado.
Ejemplo:
Consideremos primero una línea de producción y un plan de muestreo simple
como el que se indica:
N= 5000
p(%)(defectuosos promedio)=pm=4%
(n,c)= (50,4)
Plan de muestreo simple con rectificación.
Estímese el coste del plan de muestreo diseñado.
0
2
MUESTREO
3
Supongamos que tenemos “m” lotes. De esos m lotes, el número de ellos de tipo
“aceptado por haberse encontrado 0 unidades defectuosas en las n inspeccionadas” será
de m*PPoisson(0;50*0.04;FALSE). El coste de inspección para un lote de este tipo será de
n*I que se explica como sigue: la multiplicación del primer factor por el segundo
corresponde al coste de inspeccionar n elementos del lote considerado. Además, se
incurrirá en el coste (promedio) asociado debido a aquellos elementos defectuosos que
contenga el lote aceptado y que van a llegar al consumidor, que será de será (N-n)*pm*A
que se explica como sigue: la multiplicación del primer factor por el segundo
corresponde al número de elementos defectuosos (promedio) que todavía tendrá el lote
considerado (nota que se emplea N-n en lugar de n porque se supone que aquellos
elementos que se lleguen a encontrar defectuosos en la inspección de los n elementos ya
se habrán sustituido por conformes, lógicamente).
Operando de forma análoga los lotes (de esos m) que sean de tipo “aceptado por
haberse encontrado 1 unidad defectuosa en las n inspeccionadas”, y del resto de tipos de
lotes aceptados, obtenemos los costes asociados a aquellos lotes de los m que han sido
ACEPTADOS:
Coste de los lotes aceptados=
= m*(PPoisson(0)+...+ PPoisson(4)) * (n*I+(N-n)*pm*A)
Ahora nos queda estimar el coste de aquellos lotes RECHAZADOS en el plan de
muestreo planteado. De los “m” lotes iniciales, el número de ellos que será de tipo
“rechazado por haberse encontrado más de 4 unidades defectuosas en las n
inspeccionadas” será de m*(1-(PPoisson(0)+...+ PPoisson(4))). Y, para un lote de este tipo, se
incurrirá, primero, en el coste de inspección n*I y, además, SI SE HACE
RECTIFICACIÓN (lo habitual) se incurrirá en el coste adicional de inspeccionar el resto
de elementos del lote (N-n) a fin de sustituir aquellas unidades defectuosas por unidades
conformes. Es decir: (N-n)*I.
Coste de los lotes rechazados al inspeccionarlos para rectificarlos=
= m*(1-(PPoisson(0)+...+ PPoisson(4)))* (n*I+(N-n)*I)
Por lo tanto, la estimación del COSTE GLOBAL (por lote) del plan de muestreo
con rectificación planteado será:
Coste (por lote) global=
[m*(PPoisson(0)+...+ PPoisson(4)) * (n*I+(N-n)*pm*A) +
+ m*(1-(PPoisson(0)+...+ PPoisson(4)))* (n*I+(N-n)*I] * (1/m) =
Ejemplo:
Consideremos ahora una línea de producción y un plan de muestreo múltiple
como el que se indica:
N= 3000
p(%)(defectuosos promedio)=pm=4%
(n1,c1)= (7,1)
(n2,c2)= (50,2)
(n3,c3)= (80,3)
Plan de muestreo múltiple con rectificación.
Estímese el coste del plan de muestreo planteado.
1er muestreo 2do muestreo 3er muestreo
MUESTREO 0
2
1 0
3 0 0
Sean “m” lotes. En el esquema están señalados todos los posibles tipos de lotes
aceptados. Quedan fuera del esquema todos los tipos de lotes rechazados.
El número de lotes aceptados de tipo “0 unidades defectuosas localizadas en n1”
será: m*PBinomial(0;7*0.04;FALSO), que designaremos PB(0,n1), para abreviar. Para este
tipo de lotes se incurre en el gasto de inspección (n1*I) más el gasto de fallo externo ((N-
n1)*pm*A).
El número de lotes aceptados de tipo “1 unidades defectuosas en n1” será:
m*PB(1,n1). Para este tipo de lotes se incurre en el gasto de inspección (n1*I) más el
gasto de fallo externo ((N-n1)*pm*A).
El número de lotes aceptados de tipo “2 unidades defectuosas en n1 y 0 unidades
defectuosas en n2” será: m*PB(2,n1)*PPoisson(0;50*0.04;FALSE), que designaremos
como PP(0;n2), para abreviar. Para este tipo de lotes se incurre en el gasto de inspección
(n1*I+n2*I) más el gasto de fallo externo ((N-n1-n2)*pm*A).
El número de lotes aceptados de tipo “2 unidades defectuosas en n1 y 1 unidades
defectuosas en n2 y 0 unidades defectuosas en n3” será: m*PB(2,n1)*PP(1,n2)* PP(0,n3).
Para este tipo de lotes se incurre en el gasto de inspección (n1*I+n2*I+n3*I) más el
gasto de fallo externo ((N-n1-n2-n3)*pm*A).
El número de lotes aceptados de tipo “3 unidades defectuosas en n1 y 0 unidades
defectuosas en n2 y 0 unidades defectuosas en n3” será: m*PB(3,n1)*PP(0,n2)*PP(0,n3).
Para este tipo de lotes se incurre en el gasto de inspección (n1*I+n2*I+n3*I) más el
gasto de fallo externo ((N-n1-n2-n3)*pm*A).
= (n1*I+(N-n1)*pm*A)*m*(PB(0,n1)+PB(1,n1)) +
+ ((n1+n2)*I+(N-n1-n2)*pm*A)* m*PB(2,n1)*PP(0,n2) +
+ ((n1+n2+n3)*I+(N-n1-n2-n3)*pm*A)*m*(PB(2,n1)*PP(1,n2)*PP(0,n3)+
+PB(3,n1)*PP(0,n2)* PP(0,n3))
Ahora, el número de lotes rechazados de tipo “por no ser de ninguno de los tipos
de lote aceptado indicados antes” será de:
m*(1-(PB(0,n1)+PB(1,n1)+PB(2,n1)*PP(0;n2)+PB(2,n1)*PP(1,n2)* PP(0,n3)+
+PB(3,n1)*PP(0,n2)* PP(0,n3))).
Algunos de estos lotes sólo habrán sufrido la primera inspección (n1) ya que en ésta se
habrán localizado 4 o más elementos defectuosos. Para este tipo de rechazados, el coste
asociado a la inspección y rectificación posterior del lote sería: n1*I+(N-n1)*I.
Ahora bien, algunos otros habrán sufrido la primera y la segunda inspecciones porque se
habrán localizado 2 elementos defectuosos en n1 y 2 o más elementos defectuosos en n2.
Para este tipo de rechazados, el coste asociado a la inspección y rectificación posterior
del lote sería: (n1+n2)*I+(N-n1-n2)*I.
Ahora bien, algunos otros habrán experimentado las tres inspecciones al localizarse o
bien (2 elementos defectuosos en n1, 1 elemento defectuoso en n2 y 1 o más elementos
defectuosos en n3), o bien (3 elementos defectuosos en n1, 0 elementos defectuosos en
n2 y 1 o más elementos defectuosos en n3). Para este tipo de rechazados, el coste
asociado a la inspección y rectificación posterior del lote sería: (n1+n2+n3)*I+(N-n1-n2-
n3)*I.
Si quisiéramos seguir por este camino habría que cuantificar el número de lotes
que se habrían rechazado en la primera inspección, el número de ellos rechazado en la
segunda inspección y el número rechazado en la tercera ya que cada uno de ellos tendría,
a priori (o podría tener) un coste distinto. Ahora bien, si suponemos que el valor de “I” es
el mismo independientemente de la inspección realizada, podemos ver que cualquiera de
los lotes rechazados incurre en el MISMO gasto, igual a N*I con lo que podemos calcular,
PARA TODOS ELLOS CONJUNTAMENTE, el coste asociado:
Coste de los lotes rechazados al inspeccionarlos para rectificarlos=
= m*(1-(PB(0,n1)+PB(1,n1)+PB(2,n1)*PP(0;n2)+PB(2,n1)*PP(1,n2)* PP(0,n3)+
+PB(3,n1)*PP(0,n2)* PP(0,n3))) * N*I
Por lo tanto, la estimación del COSTE GLOBAL (por lote) del plan de muestreo
con rectificación planteado será:
Coste (por lote) global=
= I*[n1*(PB(0,n1)+PB(1,n1)) + (n1+n2)*PB(2,n1)*PP(0,n2) +
+ (n1+n2+n3)*(PB(2,n1)*PP(1,n2)*PP(0,n3)+PB(3,n1)*PP(0,n2)* PP(0,n3))] +
+ pm*A*[(N-n1)* (PB(0,n1)+PB(1,n1)) + (N-n1-n2)*PB(2,n1)*PP(0,n2) +
+ (N-n1-n2-n3)*(PB(2,n1)*PP(1,n2)*PP(0,n3)+PB(3,n1)*PP(0,n2)* PP(0,n3))] +
+ N*I*[(1-(PB(0,n1)+PB(1,n1)+PB(2,n1)*PP(0,n2)+
+PB(2,n1)*PP(1,n2)* PP(0,n3)+PB(3,n1)*PP(0,n2)* PP(0,n3)))]
que llamando
PAcep_1st≡PB(0,n1)+PB(1,n1)
PAcep_2nd≡ PB(2,n1)*PP(0,n2)
PAcep_3rd≡ PB(2,n1)*PP(1,n2)*PP(0,n3)+PB(3,n1)*PP(0,n2)* PP(0,n3)
PAceptación≡PAcep_1st + PAcep_2nd + PAcep_3rd
quedaría
2
MUESTREO
3
PAcep_1st≡PP(0,n)+PP(1,n)+PP(2,n)+PP(3,n)+PP(4,n)
PAceptación≡PAcep_1st
Con ese nuevo plan de muestreo diseñado, el productor sólo se engaña a sí mismo,
ya que lo único que consigue es que un lote de “buena calidad” se caracterice por un
1.8% de defectuosos (frente al 2.5% de defectuosos del otro plan, el AQL ha cambiado
poco) y, sin embargo, el consumidor lo va a acusar mucho ya que un lote de “mala
calidad” le llegará (le llegará pocas veces ya que tiene una probabilidad muy baja de ser
aceptado pero cuando le llegue) con un 19.5% de defectuosos (frente al 10% de
defectuosos del otro plan, el QL a empeorado mucho).
Un aspecto final es que se puede operar de forma análoga para obtener la curva
característica de un muestreo que sea doble o múltiple. Lo único que tenemos que hacer
es sacar la expresión para la probabilidad de aceptación global del muestreo múltiple que
sea e ir calculando sus valores resultantes, para distintos valores de p diferentes al
histórico promedio.
si e
l
(rec lote se
tifica re
ción chaza
com
plet
n a)
0(%) 0(%)
(N-n) + n
100 100
número de defectuosos a la entrada
5
Cuando N>>n, prácticamente se pueden calcular los valores de AOQ(%) como PAceptación*pentrada(%)
Así, para aquellos lotes con un porcentaje de defectuosos a la entrada de la
inspección del aprox. 7%, se obtiene la “peor mejora” pasando, de promedio, a la salida
de la inspección con un 5.08% de defectuosos (AOQL) que tampoco está mal la mejora.
Aquellos lotes con un porcentaje de defectuosos a la entrada del 3%, los deja
prácticamente como están. Sin embargo, aquellos lotes con un 10% de defectuosos, la
mejora llega hasta sacarlos, de promedio, con un 4.4% de defectuosos, que también está
muy bien. Con esta curva podemos conocer cuál será el peor porcentaje de defectuosos
con las que saldrá un lote (AOQL) en el muestreo seleccionado con rectificación incluida.
Cabe señalar que, al igual que en el apartado anterior de la determinación de la
curva característica de operación, se puede determinar de forma análoga, el AOQ (y
AOQL) de un plan de muestreo doble o múltiple que se plantee, simplemente
calculando las distintas PAceptación para cada p(%) a la entrada. La expresión para el
cálculo de AOQ (simplificada cuando N>>ni) es la misma que en la simplificada del
muestreo simple. La expresión no simplificada (que se puede obtener de forma análoga a
lo anteriormente visto) sería:
en primera en segunda en tercera
inspeccción inspeccción inspeccción
pta
p(%) p(%) p(%)
ce
(N-n1) + (N-n1-n2) + (N-n1-n2-n3) +
ea
100 100 100
es
N
lot
0(%) 0(%) 0(%)
+ n1 + (n1+n2) + (n1+n2+n3)
el
si
100 100 100
si ectif
el ica
(r
lot c
e ión
se c
n3 n2 n1
re om
ch pl
0(%) 0(%) 0(%)
az eta
N N N
a
a la entrada
Inspecc. NORMAL
- La producción está en régimen estacionario (no se - La producción se vuelve irregular o existen retrasos o
prevén cambios) y
- 1 lote se rechaza o
- 10 lotes* consecutivos aceptados y
- No se alcanza ni el número de aceptación ni el de
- El número total** de unidades defectuosas (o de rechazo de las Tables II-C, III-C, IV-C
defectos/100 unidades) es menor o igual que el número
límite de la Table VIII
Inspecc. REDUCIDA
Para cada uno de los tipos de inspección se utilizan las tablas correspondientes,
según si buscamos un muestreo simple, doble o triple. Se comienza por la Table I que nos
da un código de letra con la que entrar en las Tables II (muestreo sencillo), III (doble) o
IV (múltiple) para determinar el/los valor/es de (n,c) que determina/n el/los plan/es de
muestreo.
Una vez se ha seleccionado un plan de muestreo para tamaño de lote N, que
cumpla con el AQL requerido se puede ya operar como en los apartados 7.2.3.1 y
siguientes para calcular la curva característica de operación del plan de muestreo, su QL
(si corresponde), la curva de AOQ y determinar su AOQL y la curva de cálculo de costes
y el coste máximo esperado en ese plan.
Ejemplo:
Sea una producción de lotes de N=3000 unidades. Se requiere un nivel de
aceptación límite AQL del 1% (es decir, se desea que aquellos lotes con un 1% de
unidades defectuosas sean aceptados en la inmensa mayoría de las veces (sólo un 5% de
veces rechazados)).
Determinar un plan de muestreo acorde, de tipo sencillo (n,c) para una inspección
de grado Normal y calcular su QL, su AOQL y la curva de costes que cabe esperar si el
coste de inspeccionar una unidad, I, es de 10 euro y el coste de fallo externo, A, es de 100
euro. ¿Cuál sería el coste (por lote) promedio si el p(%) histórico de la producción, pm, ha
sido de un 2.9%?
Empezamos por la Table I, entramos por la fila para N=3000, es decir, “1201 To
3200” y vamos hasta inspección Normal y leemos la Code Letter de “K”.
Ahora nos vamos a la Table II-A correspondiente a inspección Normal y muestreo
sencillo, buscamos en la primera columna la letra “K” y nos sugiere un tamaño de
muestra de n=125. Ahora nos desplazamos en la fila hasta el AQL 1 y leemos el número
de aceptación y de rechazo: aceptación c=3, rechazo “4”. Es decir (n,c)=(125,3).
Ahora, con el plan de muestreo establecido, ya podemos construir las distintas
curvas según hicimos en los apartados 7.2.3.1 a 7.2.3.3 Se dan las condiciones para
emplear la Distribución de Poisson.
El AQL (más exactamente) es de 1.1% (que es muy próximo al deseado del 1%).
La Calidad Límite resultante del plan propuesto es de QL=5.3% (lotes con un
5.3% de defectuosos (o más) van a ser rechazados la mayoría de las veces (sólo un 10%
de veces (o menos) serán aceptados)).
La curva de porcentajes de defectuosos a la salida (AOQ(%)), después de
hacerles la rectificación correspondiente a los lotes rechazados, será6:
Y el AOQ límite, AOQL, será en torno a 1.55%. Es decir, el peor tipo de lote a la
salida, después de la rectificación, contendrá un 1.5% de defectuosos.
En cuanto a la curva de costes, aplicamos la expresión correspondiente según el
apartado 7.2.3.3:
PAcep_1st(p)= PP(0,125*p)+PP(1,125*p)+PP(2,125*p)+PP(3,125*p)
PAceptación(p)≡PAcep_1st(p)
Coste (por lote) global(p)= I*125*PAceptación(p)+ p*A*(3000-125)*PAceptación(p) +
+ 3000*I*(1-PAceptación(p))
6
Nota que aquí se ha aplicado la expresión del AOQ donde se considera N>>n, aunque se podría haber
utilizado la expresión no aproximada considerando el factor ((N-n)/N).
un lote rechazado. Pero más que este valor, nos interesará saber el coste promedio que,
para un histórico pm(%)=2.9% corresponde a un coste promedio (por lote) de 20000 euro.
Empezamos por la Table I, entramos por la fila para N=3000, es decir, “1201 To
3200” y vamos hasta inspección Reduced y leemos la Code Letter de “H”.
Ahora nos vamos a la Table II-C correspondiente a inspección Reducida y
muestreo sencillo, buscamos en la primera columna la letra “H” y nos sugiere un tamaño
de muestra de n=20. Ahora nos desplazamos en la fila hasta el AQL 1 y leemos el
número de aceptación y de rechazo: Ac=0, Rech=2. Como indica la nota, si al muestrear
un lote nos diera “1”, el lote se aceptaría pero se volvería a inspección Normal. Es decir
(n,c)=(20,1).
Ahora, con el plan de muestreo establecido, ya podemos construir las distintas
curvas según hicimos en los apartados 7.2.3.1 a 7.2.3.3 Se dan las condiciones para
emplear la Distribución de Poisson.
Y el AOQ límite,. AOQL, será en torno a 4.17%. Es decir, el peor tipo de lote a la
salida, después de la rectificación, contendrá un 4.17% de defectuosos (de promedio).
En cuanto a la curva de costes, aplicamos la expresión correspondiente según el
apartado 7.2.3.3:
Coste (por lote) global(p)= I*20*PAceptación(p)+ p*A*(3000-20)*PAceptación(p) +
+ 3000*I*(1-PAceptación(p))
Y el coste máximo estimado en que se incurrirá con el plan de muestreo elegido
(que se da a partir de un lote con p(%) defectuosos por encima del 15%) será de unos
33000 euro. Pero más que este valor, nos interesará saber el coste promedio que, para un
histórico pm(%)=2.9% corresponde a un coste promedio (por lote) de 11680 euro. Luego
mientras dure (lo poco que dure) la inspección Reducida nos ahorraremos algo del coste
de muestreo. Lo único que durará poco ya que tiene un valor de AQL muy riguroso y un
pm bastante más alejado de él que en el caso de la inspección Normal.
Podríamos trabajar con las tablas deforma análoga para inspección rigurosa y
también para probar con muestreos dobles o múltiples.
Ejemplo:
Se deja como ejercicio el probar qué plan de muestreo triple (n1,c1), (n2,c2),
(n3,c3) en inspección Normal se podría hacer para este caso: N=3000 y un AQL del 1%.
También, para los costes, I=10 euro, A=100 euro y pm(%) histórico del 2.9%.
Resolución:
- Table I. Leemos el code Setter “K”.
- Table IV-A (Continued). Leemos n1=n2=n3=32. Luego leemos (Ac1,Re1)=(”no
hay”,3); (Ac2,Re2)=(0,3); (Ac3,Re3)=(1,4). Con esos números de aceptación y rechazo,
las únicas posibilidades de aceptación serían:
1er muestreo 2do muestreo 3er muestreo
0 0
1 0
MUESTREO 1 0 0
vemos que en el primer muestreo no hay posibilidad de aceptación directa, luego todos
los elementos que haya en el 1er muestreo deberán estar unidos al siguiente muestreo.
También vemos que en 2do muestreo sólo acepta si el acumulado de defectos es 0, con lo
que sólo tendríamos la unión 0-0 como aceptable con el segundo muestreo. El resto de
opciones ya serían mayor de 0 y, por lo tanto, tendrían que estar unidas al 3er muestreo.
En éste, el acumulado no puede superar el 1, con lo que sólo quedan esas posibilidades
dibujadas. Cabe notar que, por ejemplo, el 0-2-0 no serviría porque el 3er muestreo sólo
admite 1 como acumulado. ¿Por qué aparece entonces que el rechazo es en el 2do
muestreo con 3 si nunca llegaremos a 2 o por qué aparece entonces que el rechazo es en
el 3er muestreo con 4 si nunca llegaremos a 4? Porque está hecha la tabla para poder
llegar si se quiere a 7 niveles. Cuando se establecen más niveles sí que se llega a esos
números que quedan lejos cuando el muestreo múltiple es a tres niveles.
Ahora, con el plan de muestreo establecido, ya podemos construir las distintas
curvas según hicimos en los apartados 7.2.3.1 a 7.2.3.3 Se dan las condiciones para
emplear la Distribución de Poisson.
7.3 Norma general para el control del contenido efectivo de los alimentos envasados.
(BOE 8-7-1988)
Ejemplo: ¿cuántas mujeres será necesario estudiar para estimar la prevalencia de dolor
lumbar en una población de embarazadas con un error máximo aceptado del ±8%? Nivel
de confianza 95%. Según consultas bibliográficas el dolor lumbar en ese tipo de
población afecta en torno al 20%.
Solución:
Z 2 p q 1.962 0.20 0.80
n= = = 96 mujeres
E2 (0.08) 2
Ejemplo: ¿Cuántas mujeres habrá que estudiar para determinar la media de glucemia
basal de una población con un error máximo de ±3 mg/dL? Nivel de confianza 95%.
Bibliografía: S= ±16mg/dL.
Solución:
Z 2 S 2 1.962 (16mg / dL) 2
n= = = 109 mujeres
E2 (3mg / dL) 2
n=
( ) 2 Z
Z 2 p(1 − p) + Z p1 (1 − p1 ) + p 2 (1 − p 2 )
→ valor tabulado para probabilidad β
( p1 − p 2 )2
d e cometer un falso negativo
p1 → valor de la proporción en el grupo de referencia
p 2 → valor de la proporción en el grupo de tratamiento
p = p1 + p 2
2
Ejemplo: Se diseña un ensayo clínico para evaluar si el tratamiento T2 es mejor que otro
T1 para el alivio del dolor. La eficacia del fármaco habitual es del 70% y se considera
clínicamente relevante si el nuevo fármaco alivia el dolor en un 90%. El nivel de riesgo
se fija en 0.05 (es decir, con un 5% de probabilidad de decir que existe diferencia
significativa entre ambos fármacos y que no haya tal efecto diferente entre ellos) y se
desea un poder o potencia del estudio (1-) del 80% (es decir, con una probabilidad de
0.2 de decir que no existe diferencia significativa entre ambos fármacos y que en realidad
sí actúen de forma diferente).
Solución: p=(0.70+0.90)/2 = 0.80; Z = Z =1.645; Z = Z =0.842
n=
(1.645 2 0.80 (1 − 0.80) + 0.842 0.70(1 − 0.70) + 0.90(1 − 0.90) )
2
n=
d2 d e cometer un falso negativo
S → desviaciónestándar del valor medio
d → diferencia mínima del valor medio entre los dos grupos
de tratamiento
Ejemplo: Se desea utilizar un nuevo fármaco antidiabético y se considera clínicamente
eficaz si logra un descenso de 15 mg/dL respecto al tratamiento habitual. Se conoce por
bibliografía que la desviación estándar de la glucemia en pacientes que reciben el
tratamiento habitual es de ±16 mg/dL. Se acepta un riesgo de 0.05 y se desea un poder
estadístico (1- ) del 90% para detectar diferencias si es que existen.
Solución: Z = Z =1.645; Z = Z =1.282
2 (1.645 + 1.282) 2 (16mg / dL) 2
n= = 20 pacientes en cada grupo
(−15mg / dL) 2