Articulo de referencia

Pruebas grupales

Ilustración del problema de la bombilla, donde se busca una bombilla rota entre seis. Aquí, las tres primeras están conectadas a la corriente y se encienden (A). Esto indica que...

Este es un buen artículo. Haz clic aquí para obtener más información.

Ilustración del problema de la bombilla, donde se busca una bombilla rota entre seis. Aquí, las tres primeras están conectadas a la corriente y se encienden (A). Esto indica que la bombilla rota debe ser una de las tres últimas (B). Si, en cambio, las bombillas no se encendieran, se podría estar seguro de que la bombilla rota está entre las tres primeras. Siguiendo este procedimiento, se puede localizar la bombilla rota en no más de tres pruebas, en comparación con un máximo de seis pruebas si se revisan las bombillas individualmente.

En estadística y matemáticas combinatorias , las pruebas de grupo son cualquier procedimiento que divide la tarea de identificar objetos en pruebas sobre grupos de elementos, en lugar de probar cada elemento individualmente. Estudiadas por primera vez por Robert Dorfman en 1943, las pruebas de grupo constituyen un campo relativamente nuevo de las matemáticas que puede aplicarse a una amplia gama de aplicaciones prácticas y que actualmente es un área activa de investigación.

Un ejemplo común de prueba grupal consiste en una serie de bombillas conectadas en serie, donde se sabe que una de ellas está rota. El objetivo es encontrar la bombilla rota con el menor número de pruebas posible (donde una prueba consiste en conectar algunas bombillas a una fuente de alimentación). Un método sencillo es probar cada bombilla individualmente. Sin embargo, cuando hay muchas bombillas, resulta mucho más eficiente agruparlas. Por ejemplo, al conectar la primera mitad de las bombillas a la vez, se puede determinar en qué mitad se encuentra la bombilla rota, descartando así la mitad de las bombillas en una sola prueba.

Los esquemas para realizar pruebas grupales pueden ser simples o complejos, y las pruebas involucradas en cada etapa pueden ser diferentes. Los esquemas en los que las pruebas de la siguiente etapa dependen de los resultados de las etapas anteriores se denominan procedimientos adaptativos , mientras que los esquemas diseñados de manera que todas las pruebas se conozcan de antemano se denominan procedimientos no adaptativos . La estructura del esquema de las pruebas involucradas en un procedimiento no adaptativo se conoce como diseño de agrupamiento .

Las pruebas grupales tienen muchas aplicaciones, incluyendo estadística, biología, informática, medicina, ingeniería y ciberseguridad. El interés moderno en estos esquemas de prueba se ha reavivado gracias al Proyecto Genoma Humano . [ 1 ]

Descripción básica y términos

A diferencia de muchas áreas de las matemáticas, los orígenes de las pruebas grupales se remontan a un único informe [ 2 ] escrito por una sola persona: Robert Dorfman . [ 3 ] La motivación surgió durante la Segunda Guerra Mundial , cuando el Servicio de Salud Pública de los Estados Unidos y el Servicio Selectivo emprendieron un proyecto a gran escala para descartar a todos los hombres sifilíticos llamados a filas. La prueba de sífilis consiste en extraer una muestra de sangre y analizarla para determinar la presencia o ausencia de la enfermedad. En aquel entonces, realizar esta prueba era costoso, y analizar a cada soldado individualmente habría sido muy caro e ineficiente. [ 3 ]

Suponiendo que haynorte{\displaystyle n}soldados, este método de prueba conduce anorte{\displaystyle n}Pruebas separadas. Si una gran proporción de las personas están infectadas, este método sería razonable. Sin embargo, en el caso más probable de que solo una proporción muy pequeña de los hombres esté infectada, se puede lograr un esquema de prueba mucho más eficiente. La viabilidad de un esquema de prueba más efectivo depende de la siguiente propiedad: los soldados pueden agruparse y, en cada grupo, se pueden combinar las muestras de sangre. La muestra combinada se puede analizar para comprobar si al menos un soldado del grupo tiene sífilis. Esta es la idea central detrás de las pruebas grupales. Si uno o más soldados de este grupo tienen sífilis, entonces se desperdicia una prueba (se necesitan realizar más pruebas para encontrar qué soldado(s) fue(ron)). Por otro lado, si nadie en el grupo tiene sífilis, entonces se ahorran muchas pruebas, ya que cada soldado de ese grupo puede ser descartado con solo una prueba. [ 3 ]

Los elementos que provocan que un grupo dé positivo en la prueba se denominan generalmente elementos defectuosos (estos son las bombillas rotas, los hombres sifilíticos, etc.). A menudo, el número total de elementos se denota comonorte{\displaystyle n}yd{\displaystyle d}representa el número de defectuosos si se supone que se conoce. [ 3 ]

Clasificación de los problemas de pruebas grupales

Existen dos clasificaciones independientes para los problemas de pruebas grupales; cada problema de prueba grupal es adaptativo o no adaptativo, y probabilístico o combinatorio. [ 3 ]

En los modelos probabilísticos, se supone que los artículos defectuosos siguen una distribución de probabilidad determinada y el objetivo es minimizar el número esperado de pruebas necesarias para identificar la defectuosidad de cada artículo. Por otro lado, con las pruebas de grupo combinatorias, el objetivo es minimizar el número de pruebas necesarias en el "peor escenario posible" —es decir, crear un algoritmo minmax— y no se presupone ningún conocimiento de la distribución de los artículos defectuosos. [ 3 ]

La otra clasificación, la adaptabilidad, se refiere a la información que se puede utilizar al elegir qué elementos agrupar en una prueba. En general, la elección de los elementos a probar puede depender de los resultados de pruebas anteriores, como en el problema de la bombilla mencionado anteriormente. Un algoritmo que procede realizando una prueba y luego utilizando el resultado (y todos los resultados anteriores) para decidir qué prueba realizar a continuación se denomina adaptativo. Por el contrario, en los algoritmos no adaptativos, todas las pruebas se deciden de antemano. Esta idea se puede generalizar a algoritmos multietapa, donde las pruebas se dividen en etapas y cada prueba en la siguiente etapa debe decidirse de antemano, conociendo únicamente los resultados de las pruebas en etapas anteriores. Si bien los algoritmos adaptativos ofrecen mucha más libertad en el diseño, se sabe que los algoritmos adaptativos de pruebas grupales no mejoran a los no adaptativos más que un factor constante en el número de pruebas necesarias para identificar el conjunto de elementos defectuosos. [ 4 ] [ 3 ] Además, los métodos no adaptativos suelen ser útiles en la práctica porque se pueden realizar pruebas sucesivas sin analizar primero los resultados de todas las pruebas anteriores, lo que permite una distribución eficaz del proceso de prueba. [ 5 ]

Variaciones y extensiones

Hay muchas maneras de extender el problema de las pruebas grupales. Una de las más importantes se llama prueba grupal ruidosa y aborda una suposición importante del problema original: que la prueba está libre de errores. Un problema de prueba grupal se llama ruidoso cuando hay alguna posibilidad de que el resultado de una prueba grupal sea erróneo (por ejemplo, que dé positivo cuando la prueba no contenía defectuosos). El modelo de ruido de Bernoulli supone que esta probabilidad es una constante.q{\displaystyle q}, pero en general puede depender del número real de defectuosos en la prueba y del número de elementos probados. [ 6 ] Por ejemplo, el efecto de dilución puede modelarse diciendo que un resultado positivo es más probable cuando hay más defectuosos (o más defectuosos como fracción del número probado) presentes en la prueba. [ 7 ] Un algoritmo ruidoso siempre tendrá una probabilidad no nula de cometer un error (es decir, etiquetar incorrectamente un elemento). [ 6 ]

Las pruebas grupales pueden ampliarse considerando escenarios en los que haya más de dos resultados posibles de una prueba. Por ejemplo, una prueba puede tener los siguientes resultados:0,1{\displaystyle 0,1}y2+{\displaystyle 2^{+}}, lo que corresponde a que no haya defectuosos, un solo defectuoso o un número desconocido de defectuosos mayor que uno. De manera más general, es posible considerar que el conjunto de resultados de una prueba es0,1,,k+{\displaystyle {0,1,\ldots ,k^{+}}}para algunosknorte{\displaystyle k\in \mathbb {N} }. [ 8 ]

Otra extensión consiste en considerar restricciones geométricas sobre qué conjuntos pueden ser probados. El problema de la bombilla mencionado anteriormente es un ejemplo de este tipo de restricción: solo se pueden probar las bombillas que aparecen consecutivamente. De manera similar, los elementos pueden estar dispuestos en un círculo o, en general, en una red, donde las pruebas son caminos disponibles en el grafo. Otro tipo de restricción geométrica sería sobre el número máximo de elementos que pueden ser probados en un grupo, [ a ] o los tamaños de los grupos podrían tener que ser pares, etc. De manera similar, puede ser útil considerar la restricción de que cualquier elemento dado solo puede aparecer en un cierto número de pruebas. [ 8 ]

Existen innumerables maneras de seguir modificando la fórmula básica de las pruebas grupales. Las siguientes explicaciones ilustran algunas de las variantes más singulares. En el modelo «bueno-mediocre-malo», cada elemento es «bueno», «mediocre» o «malo», y el resultado de la prueba es el tipo del elemento «peor» del grupo. En las pruebas grupales con umbral, el resultado es positivo si el número de elementos defectuosos en el grupo supera un valor o proporción umbral. [ 9 ] Las pruebas grupales con inhibidores son una variante con aplicaciones en biología molecular. En este caso, existe una tercera clase de elementos denominados inhibidores, y el resultado es positivo si contiene al menos un elemento defectuoso y ningún inhibidor. [ 10 ]

Historia y desarrollo

Invención y progreso inicial

El concepto de pruebas grupales fue introducido por primera vez por Robert Dorfman en 1943 en un breve informe [ 2 ] publicado en la sección de Notas de Annals of Mathematical Statistics . [ 8 ] [ b ] El informe de Dorfman, al igual que todos los primeros trabajos sobre pruebas grupales, se centró en el problema probabilístico y tuvo como objetivo utilizar la novedosa idea de las pruebas grupales para reducir el número esperado de pruebas necesarias para descartar a todos los hombres sifilíticos en un grupo determinado de soldados. El método era simple: colocar a los soldados en grupos de un tamaño determinado y utilizar pruebas individuales (pruebas de elementos en grupos de tamaño uno) en los grupos positivos para encontrar cuáles estaban infectados. Dorfman tabuló los tamaños de grupo óptimos para esta estrategia en función de la tasa de prevalencia de la deficiencia en la población. [ 2 ] Stephen Samuels encontró una solución de forma cerrada para el tamaño de grupo óptimo en función de la tasa de prevalencia. [ 12 ]

Después de 1943, las pruebas grupales permanecieron prácticamente sin cambios durante varios años. Luego, en 1957, Sterrett introdujo una mejora al procedimiento de Dorfman. Este nuevo proceso comienza realizando nuevamente pruebas individuales en los grupos positivos, pero deteniéndose tan pronto como se identifica un elemento defectuoso. Luego, los elementos restantes del grupo se prueban juntos, ya que es muy probable que ninguno de ellos sea defectuoso. [ 13 ]

El primer análisis exhaustivo de las pruebas grupales fue realizado por Sobel y Groll en su artículo fundamental de 1959 sobre el tema. Describieron cinco nuevos procedimientos —además de generalizaciones para casos donde se desconoce la tasa de prevalencia— y, para el óptimo, proporcionaron una fórmula explícita para el número esperado de pruebas que requeriría. El artículo también estableció por primera vez la conexión entre las pruebas grupales y la teoría de la información , además de analizar varias generalizaciones del problema de las pruebas grupales y ofrecer algunas aplicaciones nuevas de la teoría. [ 14 ]

El resultado fundamental de Peter Ungar en 1960 muestra que si la tasa de prevalenciapag>pag{\displaystyle p>p_{u}}, dóndepag=(35)/20,38{\displaystyle p_{u}=(3-{\sqrt {5}})/2\approx 0.38}, entonces la prueba individual es el procedimiento óptimo de prueba grupal con respecto al número esperado de pruebas, y sipag<pag{\displaystyle p<p_{u}}, entonces no es óptimo. Sin embargo, es importante señalar que a pesar de 80 años de esfuerzo de investigación, el procedimiento óptimo aún se desconoce parapag<pag{\displaystyle p<p_{u}}y un tamaño de población generalnorte>2{\displaystyle n>2}. [ 15 ]

Pruebas de grupo combinatorias

Las pruebas de grupo fueron estudiadas por primera vez en el contexto combinatorio por Li en 1962, [ 16 ] con la introducción de la de Li.s{\displaystyle s}algoritmo de -etapa . [ 8 ] Li propuso una extensión del 'algoritmo de 2 etapas' de Dorfman a un número arbitrario de etapas que no requería más det=miregistro2(mi)dregistro2(norte){\textstyle t={\frac {e}{\log _{2}(e)}}d\log _{2}(n)}pruebas que se garantizan para encontrard{\displaystyle d}o menos defectuosos entrenorte{\displaystyle n}elementos. La idea era eliminar todos los elementos con resultados negativos en las pruebas y dividir los elementos restantes en grupos, como se hizo con el conjunto inicial. Esto se iba a hacers1{\displaystyle s-1}veces antes de realizar pruebas individuales. [ 16 ]

Las pruebas de grupo combinatorias en general fueron estudiadas más a fondo por Katona en 1973. Katona introdujo la representación matricial de las pruebas de grupo no adaptativas y desarrolló un procedimiento para encontrar el defectuoso en el caso no adaptativo de 1-defectuoso en no más det=registro2(norte){\displaystyle t=\lceil \log _{2}(n)\rceil }pruebas, que también demostró que eran óptimas. [ 17 ]

En general, encontrar algoritmos óptimos para pruebas grupales combinatorias adaptativas es difícil, y aunque no se ha determinado la complejidad computacional de las pruebas grupales, se sospecha que es difícil en alguna clase de complejidad . [ 8 ] Sin embargo, en 1972 se produjo un avance importante con la introducción del algoritmo generalizado de división binaria . Este algoritmo funciona realizando una búsqueda binaria en grupos que dan positivo en las pruebas, y es un algoritmo simple que encuentra un único defectuoso en no más del número de pruebas que el límite inferior de información . [ 18 ]

En escenarios donde hay dos o más defectuosos, el algoritmo generalizado de división binaria aún produce resultados casi óptimos, requiriendo como máximod1{\displaystyle d-1}pruebas por encima del límite inferior de información donded{\displaystyle d}es el número de defectuosos. [ 18 ] Allemann realizó mejoras considerables en esto en 2013, logrando que el número requerido de pruebas fuera menor que0,187d+0,5registro2(d)+5.5{\displaystyle 0.187d+0.5\log _{2}(d)+5.5}por encima del límite inferior de información cuandonorte/d38{\displaystyle n/d\geq 38}yd10{\displaystyle d\geq 10}Esto se logró cambiando la búsqueda binaria en el algoritmo de división binaria a un conjunto complejo de subalgoritmos con grupos de prueba superpuestos. De esta manera, el problema de las pruebas de grupo combinatorias adaptativas —con un número conocido o un límite superior en el número de defectuosos— se ha resuelto esencialmente, con poco margen para futuras mejoras. [ 19 ]

Existe una cuestión abierta sobre cuándo las pruebas individuales son minmax . Hu, Hwang y Wang demostraron en 1981 que las pruebas individuales son minmax cuandonorte(5d+1)/2{\displaystyle n\leq \lfloor (5d+1)/2\rfloor }y que no es minmax cuandonorte>3d{\displaystyle n>3d}. [ 20 ] Actualmente se conjetura que este límite es preciso: es decir, la prueba individual es minmax si y solo sinorte3d{\displaystyle n\leq 3d}. [ 21 ] [ c ] En 2000, Riccio y Colbourn lograron algunos avances al demostrar que para grandesnorte{\displaystyle n}, la prueba individual es minmax cuandodnorte/registro3/2(3)0,369norte{\displaystyle d\geq n/\log _{3/2}(3)\approx 0.369n}. [ 22 ]

Pruebas no adaptativas y probabilísticas

Una de las ideas clave en las pruebas grupales no adaptativas es que se pueden lograr avances significativos eliminando el requisito de que el procedimiento de prueba grupal tenga éxito garantizado (el problema "combinatorio"), y en cambio permitiéndole tener una probabilidad baja pero no nula de etiquetar erróneamente cada elemento (el problema "probabilístico"). Se sabe que a medida que el número de elementos defectuosos se acerca al número total de elementos, las soluciones combinatorias exactas requieren significativamente más pruebas que las soluciones probabilísticas, incluso cuando estas últimas permiten solo una probabilidad de error asintóticamente pequeña. [ 4 ] [ d ]

En este sentido, Chan et al. (2011) introdujeron COMP , un algoritmo probabilístico que no requiere más quet=mid(1+δ)ln(norte){\displaystyle t=ed(1+\delta )\ln(n)}pruebas para encontrar hastad{\displaystyle d}defectuosos ennorte{\displaystyle n}elementos con una probabilidad de error no mayor quenorteδ{\displaystyle n^{-\delta }}. [ 6 ] Esto está dentro de un factor constante de lat=O(dregistro2norte){\displaystyle t=O(d\log _{2}n)}límite inferior. [ 4 ]

Chan et al. (2011) también proporcionaron una generalización de COMP a un modelo ruidoso simple y, de manera similar, produjeron un límite de rendimiento explícito, que nuevamente era solo una constante (dependiente de la probabilidad de una prueba fallida) por encima del límite inferior correspondiente. [ 4 ] [ 6 ] En general, el número de pruebas requeridas en el caso de ruido de Bernoulli es un factor constante mayor que en el caso sin ruido. [ 6 ]

Aldridge, Baldassini y Johnson (2014) produjeron una extensión del algoritmo COMP que agregó pasos de posprocesamiento adicionales. [ 23 ] Demostraron que el rendimiento de este nuevo algoritmo, llamado DD , supera estrictamente el de COMP, y que DD es "esencialmente óptimo" en escenarios donded2norte{\displaystyle d^{2}\geq n}, comparándolo con un algoritmo hipotético que define un óptimo razonable. El rendimiento de este algoritmo hipotético sugiere que hay margen de mejora cuandod2<norte{\displaystyle d^{2}<n}, además de sugerir cuánto podría mejorar esto. [ 23 ]

Formalización de las pruebas de grupos combinatorios

Esta sección define formalmente las nociones y los términos relacionados con las pruebas grupales.

  • El vector de entrada ,incógnita=(incógnita1,incógnita2,,incógnitanorte){\displaystyle \mathbf {x} =(x_{1},x_{2},\dots ,x_{n})}, se define como un vector binario de longitudnorte{\displaystyle n}(eso es,incógnita{0,1}norte{\displaystyle \mathbf {x} \in \{0,1\}^{n}}), y el j -ésimo artículo se considera defectuoso si y solo siincógnitaj=1{\displaystyle x_{j}=1}Además, cualquier artículo que no presente defectos se denomina artículo "bueno".

incógnita{\displaystyle \mathbf {x} }tiene como objetivo describir el conjunto (desconocido) de elementos defectuosos. La propiedad clave deincógnita{\displaystyle \mathbf {x} }es que es una entrada implícita . Es decir, no hay conocimiento directo de cuáles son las entradas deincógnita{\displaystyle \mathbf {x} }son, aparte de lo que se puede inferir mediante alguna serie de 'pruebas'. Esto nos lleva a la siguiente definición.

  • Dejarincógnita{\displaystyle \mathbf {x} }sea ​​un vector de entrada. Un conjunto,S{1,2,,norte}{\displaystyle S\subseteq \{1,2,\dots ,n\}}se llama prueba . Cuando la prueba no tiene ruido , el resultado de una prueba es positivo cuando existejS{\displaystyle j\in S}de tal manera queincógnitaj=1{\displaystyle x_{j}=1}y, en caso contrario, el resultado es negativo .

Por lo tanto, el objetivo de las pruebas grupales es encontrar un método para elegir una serie "corta" de pruebas que permitanincógnita{\displaystyle \mathbf {x} }Por determinar, ya sea con exactitud o con un alto grado de certeza.

  • Se dice que un algoritmo de prueba grupal comete un error si etiqueta incorrectamente un elemento (es decir, etiqueta un elemento defectuoso como no defectuoso o viceversa). Esto no es lo mismo que el resultado de una prueba grupal sea incorrecto. Un algoritmo se denomina de error cero si la probabilidad de que cometa un error es cero. [ e ]
  • t(d,norte){\displaystyle t(d,n)}denota el número mínimo de pruebas necesarias para encontrar siempred{\displaystyle d}defectuosos entrenorte{\displaystyle n}elementos con probabilidad cero de error por cualquier algoritmo de prueba grupal. Para la misma cantidad pero con la restricción de que el algoritmo no sea adaptativo, la notaciónt¯(d,norte){\displaystyle {\bar {t}}(d,n)}se utiliza.

límites generales

Dado que siempre es posible recurrir a pruebas individuales mediante la configuraciónSj={j}{\displaystyle S_{j}=\{j\}}para cada1jnorte{\displaystyle 1\leq j\leq n}, debe ser que esot¯(d,norte)norte{\displaystyle {\bar {t}}(d,n)\leq n}Además, dado que cualquier procedimiento de prueba no adaptativo puede escribirse como un algoritmo adaptativo simplemente realizando todas las pruebas sin tener en cuenta su resultado,t(d,norte)t¯(d,norte){\displaystyle t(d,n)\leq {\bar {t}}(d,n)}. Finalmente, cuando0dnorte{\displaystyle 0\neq d\neq n}, hay al menos un artículo cuya defectuosidad debe determinarse (mediante al menos una prueba), y por lo tanto1t(d,norte){\displaystyle 1\leq t(d,n)}.

En resumen (al asumir0dnorte{\displaystyle 0\neq d\neq n}),1t(d,norte)t¯(d,norte)norte{\displaystyle 1\leq t(d,n)\leq {\bar {t}}(d,n)\leq n}. [ f ]

Límite inferior de información

Se puede describir un límite inferior en el número de pruebas necesarias utilizando la noción de espacio muestral , denotadoS{\displaystyle {\mathcal {S}}}, que es simplemente el conjunto de posibles ubicaciones de los defectuosos. Para cualquier problema de prueba de grupo con espacio de muestraS{\displaystyle {\mathcal {S}}}y cualquier algoritmo de prueba grupal, se puede demostrar quetregistro2|S|{\displaystyle t\geq \lceil \log _{2}{|{\mathcal {S}}|}\rceil }, dóndet{\displaystyle t}es el número mínimo de pruebas necesarias para identificar todos los defectuosos con una probabilidad de error cero. Esto se denomina límite inferior de información . [ 8 ] Este límite se deriva del hecho de que después de cada prueba,S{\displaystyle {\mathcal {S}}}se divide en dos subconjuntos disjuntos, cada uno correspondiente a uno de los dos posibles resultados de la prueba.

Sin embargo, el límite inferior de información en sí mismo suele ser inalcanzable, incluso para problemas pequeños. [ 8 ] Esto se debe a que la división deS{\displaystyle {\mathcal {S}}}no es arbitrario, puesto que debe ser realizable mediante alguna prueba.

De hecho, el límite inferior de información se puede generalizar al caso en que existe una probabilidad no nula de que el algoritmo cometa un error. En esta forma, el teorema nos da un límite superior en la probabilidad de éxito basado en el número de pruebas. Para cualquier algoritmo de prueba grupal que realicet{\displaystyle t}pruebas, la probabilidad de éxito,PAG(éxito){\displaystyle \mathbb {P} ({\textrm {success}})}, satisfacePAG(éxito)t/registro2(norted){\displaystyle \mathbb {P} ({\textrm {success}})\leq t/\log _{2}{n \choose d}}Esto se puede reforzar para:PAG(éxito)2t(norted){\displaystyle \mathbb {P} ({\textrm {success}})\leq {\frac {2^{t}}{n \choose d}}}. [ 6 ] [ 24 ]

Representación de algoritmos no adaptativos

Diagrama que muestra una matriz de pruebas grupales junto con los vectores asociados, x e y.
Una configuración típica de prueba grupal. Un algoritmo no adaptativo primero elige la matriz.METRO{\displaystyle M}y luego se le da el vector y . El problema es entonces encontrar una estimación para x .

Los algoritmos para pruebas grupales no adaptativas constan de dos fases distintas. Primero, se decide cuántas pruebas realizar y qué elementos incluir en cada una. En la segunda fase, a menudo denominada etapa de decodificación, se analizan los resultados de cada prueba grupal para determinar qué elementos tienen más probabilidades de ser defectuosos. La primera fase se suele codificar en una matriz como se muestra a continuación. [ 6 ]

  • Supongamos un procedimiento de prueba de grupo no adaptativo paranorte{\displaystyle n}Los elementos consisten en las pruebasS1,S2,,St{\displaystyle S_{1},S_{2},\dots ,S_{t}}para algunostnorte0{\displaystyle t\in \mathbb {N} _{\geq 0}}. La matriz de pruebas para este esquema es lat×norte{\displaystyle t\times n}matriz binaria,METRO{\displaystyle M}, dónde(METRO)ij=1{\displaystyle (M)_{ij}=1}si y solo sijSi{\displaystyle j\in S_{i}}(y es cero en caso contrario).

Así, cada columna deMETRO{\displaystyle M}representa un elemento y cada fila representa una prueba, con un1{\displaystyle 1}en el(i,j)-th{\displaystyle (i,j){\textrm {-th}}}entrada que indica que eli-th{\displaystyle i{\textrm {-th}}}La prueba incluyó laj-th{\displaystyle j{\textrm {-th}}}artículo y un0{\displaystyle 0}indicando lo contrario.

Además del vectorincógnita{\displaystyle \mathbf {x} }(de longitudnorte{\displaystyle n}) que describe el conjunto defectuoso desconocido, es común introducir el vector de resultados, que describe los resultados de cada prueba.

  • Dejart{\displaystyle t}sea ​​el número de pruebas realizadas por un algoritmo no adaptativo. El vector de resultados ,y=(y1,y2,,yt){\displaystyle \mathbf {y} =(y_{1},y_{2},\dots ,y_{t})}, es un vector binario de longitudt{\displaystyle t}(eso es,y{0,1}t{\displaystyle \mathbf {y} \in \{0,1\}^{t}}) tal queyi=1{\displaystyle y_{i}=1}si y solo si el resultado de lai-th{\displaystyle i{\textrm {-th}}}La prueba fue positiva (es decir, contenía al menos un defectuoso). [ g ]

Con estas definiciones, el problema no adaptativo puede reformularse de la siguiente manera: primero se elige una matriz de prueba,METRO{\displaystyle M}, después de lo cual el vectory{\displaystyle \mathbf {y} }se devuelve. Entonces el problema es analizary{\displaystyle \mathbf {y} }para encontrar alguna estimación deincógnita{\displaystyle \mathbf {x} }.

En el caso ruidoso más simple, donde hay una probabilidad constante,q{\displaystyle q}, que una prueba grupal tendrá un resultado erróneo, se considera un vector binario aleatorio,v{\displaystyle \mathbf {v} }donde cada entrada tiene una probabilidadq{\displaystyle q}de ser1{\displaystyle 1}y es0{\displaystyle 0}de lo contrario. El vector que se devuelve es entoncesy^=y+v{\displaystyle {\hat {\mathbf {y} }}=\mathbf {y} +\mathbf {v} }, con la adición habitual en(Z/2Z)norte{\displaystyle (\mathbb {Z} /2\mathbb {Z} )^{n}}(equivalentemente, esta es la operación XOR elemento a elemento ). Un algoritmo ruidoso debe estimarincógnita{\displaystyle \mathbf {x} }usandoy^{\displaystyle {\hat {\mathbf {y} }}}(es decir, sin conocimiento directo dey{\displaystyle \mathbf {y} }). [ 6 ]

Límites para algoritmos no adaptativos

La representación matricial permite demostrar algunos límites en las pruebas grupales no adaptativas. El enfoque refleja el de muchos diseños deterministas, donded{\displaystyle d}Se consideran matrices separables, tal como se define a continuación. [ 8 ]

  • Una matriz binaria,METRO{\displaystyle M}, se llamad{\displaystyle d}-separable si cada suma booleana (OR lógico) de cualquierd{\displaystyle d}de sus columnas es distinto. Además, la notaciónd¯{\displaystyle {\bar {d}}}-separable indica que cada suma de cualquiera de hastad{\displaystyle d}deMETRO{\displaystyle M}Las columnas de 's son distintas. (Esto no es lo mismo queMETRO{\displaystyle M}serk{\displaystyle k}-separable para cadakd{\displaystyle k\leq d}.)

CuandoMETRO{\displaystyle M}es una matriz de prueba, la propiedad de serd{\displaystyle d}-separable (d¯{\displaystyle {\bar {d}}}-separable) es equivalente a poder distinguir entre (hasta)d{\displaystyle d}defectuosos. Sin embargo, esto no garantiza que sea sencillo. Una propiedad más fuerte, llamada disyunción, sí lo garantiza.

  • Una matriz binaria,METRO{\displaystyle M}se llamad{\displaystyle d}-disyunto si la suma booleana de cualquierd{\displaystyle d}Una columna no contiene ninguna otra columna. (En este contexto, se dice que una columna A contiene una columna B si para cada índice donde B tiene un 1, A también tiene un 1).

Una propiedad útil ded{\displaystyle d}-las matrices de prueba disjuntas son que, con hastad{\displaystyle d}En el caso de los artículos defectuosos, cada artículo no defectuoso aparecerá en al menos una prueba cuyo resultado sea negativo. Esto significa que existe un procedimiento sencillo para encontrar los defectuosos: basta con retirar todos los artículos que den negativo en una prueba.

Utilizando las propiedades ded{\displaystyle d}-separable yd{\displaystyle d}-matrices disjuntas lo siguiente se puede demostrar para el problema de identificaciónd{\displaystyle d}defectuosos entrenorte{\displaystyle n}Artículos totales. [ 4 ]

  1. El número de pruebas necesarias para una probabilidad de error promedio asintóticamente pequeña aumenta comoO(dregistro2norte){\displaystyle O(d\log _{2}n)}.
  2. El número de pruebas necesarias para una probabilidad máxima de error asintóticamente pequeña se escala comoO(d2registro2norte){\displaystyle O(d^{2}\log _{2}n)}.
  3. El número de pruebas necesarias para una probabilidad de error cero se escala comoO(d2registro2norteregistro2d){\displaystyle O\left({\frac {d^{2}\log _{2}n}{\log _{2}d}}\right)}.

Algoritmo generalizado de división binaria

Una ilustración del algoritmo generalizado de división binaria donde hay 8 defectuosos y 135 artículos en total. Aquí,2α1=16{\displaystyle 2^{\alpha _{1}}=16}y la primera prueba da un resultado negativo, por lo que todos los artículos se declaran no defectuosos. Por lo tanto, quedan 119 artículos, así que2α2=8{\displaystyle 2^{\alpha _{2}}=8}Este segundo grupo da un resultado positivo, por lo que se utiliza una búsqueda binaria para encontrar un defectuoso. Una vez hecho esto, se repite todo el proceso, calculando un nuevoα{\displaystyle \alpha }Utilizar únicamente aquellos artículos cuyo defecto no se haya determinado.

El algoritmo generalizado de división binaria es un algoritmo de prueba grupal adaptativo esencialmente óptimo que encuentrad{\displaystyle d}o menos defectuosos entrenorte{\displaystyle n}elementos como sigue: [ 8 ] [ 18 ]

  1. Sinorte2d2{\displaystyle n\leq 2d-2}, prueba elnorte{\displaystyle n}artículos individualmente. De lo contrario, configurel=norted+1{\displaystyle l=n-d+1}yα=registro2l/d{\displaystyle \alpha =\lfloor \log _{2}{l/d}\rfloor }.
  2. Prueba un grupo de tamaño2α{\displaystyle 2^{\alpha }}. Si el resultado es negativo, cada elemento del grupo se declara no defectuoso; conjuntonorte:=norte2α{\displaystyle n:=n-2^{\alpha }}y vaya al paso 1. De lo contrario, utilice una búsqueda binaria para identificar un defectuoso y un número no especificado, llamadoincógnita{\displaystyle x}, de artículos no defectuosos; conjuntonorte:=norte1incógnita{\displaystyle n:=n-1-x}yd:=d1{\displaystyle d:=d-1}. Vaya al paso 1.

El algoritmo generalizado de división binaria no requiere más queT{\displaystyle T}pruebas donde T={nortenorte2d2(α+2)d+pag1norte2d1{\displaystyle T={\begin{cases}n&n\leq 2d-2\\(\alpha +2)d+p-1&n\geq 2d-1\end{cases}}}. [ 8 ]

Paranorte/d{\displaystyle n/d}grande, se puede demostrar queTdregistro2(norte/d){\displaystyle T\rightarrow d\log _{2}(n/d)}, [ 8 ] que se compara favorablemente con elt=miregistro2midregistro2(norted){\displaystyle t={\frac {e}{\log _{2}e}}d\log _{2}\left({\frac {n}{d}}\right)}pruebas requeridas para Lis{\displaystyle s}-algoritmo de etapa. De hecho, el algoritmo generalizado de división binaria es casi óptimo en el siguiente sentido. Cuandod2{\displaystyle d\geq 2}Se puede demostrar queTBI(d,norte)(d1){\displaystyle T-B_{I}(d,n)\leq (d-1)}, dóndeBI(d,norte)=registro2i=0d(nortei){\displaystyle B_{I}(d,n)=\left\lceil \log _{2}\sum _{i=0}^{d}{n \choose i}\right\rceil }es el límite inferior de la información. [ 8 ] [ 18 ]

Algoritmos no adaptativos

Los algoritmos de prueba grupal no adaptativos tienden a asumir que se conoce el número de defectuosos, o al menos un buen límite superior para ellos. [ 6 ] Esta cantidad se denotad{\displaystyle d}en esta sección. Si no se conocen los límites, existen algoritmos no adaptativos con baja complejidad de consulta que pueden ayudar a estimard{\displaystyle d}. [ 25 ]

Búsqueda de coincidencia ortogonal combinatoria (COMP)

Ilustración del algoritmo COMP. COMP identifica el elemento a como defectuoso y el elemento b como no defectuoso. Sin embargo, etiqueta erróneamente a c como defectuoso, ya que queda "oculto" por los elementos defectuosos en todas las pruebas en las que aparece.

El algoritmo Combinatorial Orthogonal Matching Pursuit , o COMP, es un algoritmo sencillo de prueba de grupos no adaptativo que constituye la base de los algoritmos más complejos que se describen a continuación en esta sección.

Primero, cada entrada de la matriz de prueba se elige iid para ser1{\displaystyle 1}con probabilidad1/d{\displaystyle 1/d}y0{\displaystyle 0}de lo contrario.

El paso de decodificación procede columna por columna (es decir, por elemento). Si cada prueba en la que aparece un elemento es positiva, entonces el elemento se declara defectuoso; de lo contrario, se asume que el elemento no es defectuoso. O equivalentemente, si un elemento aparece en alguna prueba cuyo resultado es negativo, el elemento se declara no defectuoso; de lo contrario, se asume que el elemento es defectuoso. Una propiedad importante de este algoritmo es que nunca crea falsos negativos , aunque se produce un falso positivo cuando todas las ubicaciones con unos en la j -ésima columna deMETRO{\displaystyle M}(que corresponden a un artículo no defectuoso j ) están "ocultos" por los de otras columnas que corresponden a artículos defectuosos.

El algoritmo COMP no requiere más quemid(1+δ)ln(norte){\displaystyle ed(1+\delta )\ln(n)}pruebas para tener una probabilidad de error menor o igual anorteδ{\displaystyle n^{-\delta }}. [ 6 ] Esto está dentro de un factor constante del límite inferior para la probabilidad promedio de error anterior.

En el caso ruidoso, se relaja el requisito en el algoritmo COMP original de que el conjunto de ubicaciones de unos en cualquier columna deMETRO{\displaystyle M}correspondiente a un elemento positivo estar completamente contenido en el conjunto de ubicaciones de unos en el vector de resultados. En cambio, se permite un cierto número de “desajustes”: este número de desajustes depende tanto del número de unos en cada columna como del parámetro de ruido.q{\displaystyle q}Este algoritmo COMP ruidoso no requiere más que4.36(δ+1+δ)2(12q)2dregistro2norte{\displaystyle 4.36({\sqrt {\delta }}+{\sqrt {1+\delta }})^{2}(1-2q)^{-2}d\log _{2}{n}}pruebas para lograr una probabilidad de error como máximonorteδ{\displaystyle n^{-\delta }}. [ 6 ]

Defectuosos definitivos (DD)

El método de defectuosos definidos (DD) es una extensión del algoritmo COMP que intenta eliminar cualquier falso positivo. Se ha demostrado que las garantías de rendimiento de DD superan con creces las de COMP. [ 23 ]

El paso de decodificación utiliza una propiedad útil del algoritmo COMP: que cada elemento que COMP declara como no defectuoso ciertamente no lo es (es decir, no hay falsos negativos). Procede de la siguiente manera.

  1. Primero se ejecuta el algoritmo COMP y se eliminan los artículos que no presentan defectos. Todos los artículos restantes se consideran "posiblemente defectuosos".
  2. A continuación, el algoritmo analiza todas las pruebas positivas. Si un elemento aparece como el único "posible defectuoso" en una prueba, entonces debe ser defectuoso, por lo que el algoritmo lo declara defectuoso.
  3. Se supone que todos los demás artículos no presentan defectos. La justificación de este último paso radica en la suposición de que el número de artículos defectuosos es mucho menor que el número total de artículos.

Cabe destacar que los pasos 1 y 2 nunca cometen errores, por lo que el algoritmo solo puede equivocarse si declara que un artículo defectuoso no lo es. Por lo tanto, el algoritmo DD solo puede generar falsos negativos.

COMP secuencial (SCOMP)

SCOMP (Sequential COMP) es un algoritmo que utiliza el hecho de que DD no comete errores hasta el último paso, donde se supone que los elementos restantes no son defectuosos. Sea el conjunto de elementos declarados defectuosos.K{\displaystyle K}Una prueba positiva se llama explicada porK{\displaystyle K}si contiene al menos un elemento enK{\displaystyle K}La observación clave con SCOMP es que el conjunto de defectos encontrados por DD puede no explicar todas las pruebas positivas, y que cada prueba inexplicable debe contener un defecto oculto.

El algoritmo procede de la siguiente manera.

  1. Realice los pasos 1 y 2 del algoritmo DD para obtenerK{\displaystyle K}, una estimación inicial para el conjunto de defectuosos.
  2. SiK{\displaystyle K}Explica cada prueba positiva, finaliza el algoritmo:K{\displaystyle K}es la estimación final para el conjunto de defectuosos.
  3. Si hay alguna prueba inexplicada, encuentre el "posible defectuoso" que aparece en la mayor cantidad de pruebas inexplicadas y declárelo como defectuoso (es decir, agréguelo al conjunto).K{\displaystyle K}). Vaya al paso 2.

En simulaciones, se ha demostrado que SCOMP tiene un rendimiento casi óptimo. [ 23 ]

Pozos polinomiales (PP)

Polynomial Pools (PP) es un algoritmo determinista que garantiza identificar con exactitud hastad{\displaystyle d}positivos. [ 26 ] El algoritmo es para la construcción de la matriz de agrupaciónMETRO{\displaystyle M}, que se puede utilizar directamente para decodificar las observaciones eny{\displaystyle y}. De forma similar a COMP, una muestra se decodifica según la relación: incógnitai=1   si   METRO(:,i) . y=METRO(:,i){\displaystyle x_{i}=1~~{\text{ if }}~~M(:,i)~.*~y=M(:,i)}, dónde.{\displaystyle .*}representa la multiplicación elemento a elemento yMETRO(:,i){\displaystyle M(:,i)}es eli{\displaystyle i}columna deMETRO{\displaystyle M}. Dado que el paso de decodificación no es difícil, PP está especializado en generarMETRO{\displaystyle M}.

Formación de grupos

Diseño de un grupo conqdo1=9{\displaystyle q^{c-1}=9}muestras (azul) de un conjunto denorte=qdo=27{\displaystyle n=q^{c}=27}muestras totales utilizando el algoritmo de grupos polinomiales

Un grupo/piscina{\displaystyle \ell }Se genera utilizando una relación polinómica que especifica los índices de las muestras contenidas en cada grupo. Un conjunto de parámetros de entrada determina el algoritmo. Para un número primopag>1{\displaystyle p>1}y un número enteronorte1{\displaystyle n\geq 1}cualquier potencia prima se define porq=pagnorte{\displaystyle q=p^{n}}. Para un parámetro de dimensióndo2{\displaystyle c\geq 2}El número total de muestras esnorte=qdo{\displaystyle n=q^{c}}y el número de muestras por grupo esqdo1{\displaystyle q^{c-1}}. Además, el campo finito de ordenq{\displaystyle q}se denota porFq{\displaystyle \mathbb {F} _{q}} (es decir, los números enteros){0,1,2,,q1}{\displaystyle \{0,1,2,\ldots ,q-1\}}definido por operaciones aritméticas especiales que aseguran que la suma y la multiplicación enFq{\displaystyle \mathbb {F} _{q}}permanece enFq{\displaystyle \mathbb {F} _{q}}). El método dispone cada muestra en una cuadrícula y la representa mediante coordenadasincógnita=(,v){\displaystyle x=(u,v)}Las coordenadas se calculan según una relación polinómica utilizando los números enteros. 1ldo1{\displaystyle 1\leq l\leq c-1},0ilq1{\displaystyle 0\leq u_{i_{l}}\leq q-1}

v = ado1 ido1++a i1+b,a,b,ilFq.{\displaystyle v~=~a^{c-1}~u_{i_{c-1}}+\cdots +a~u_{i_{1}}+b,\quad a,b,u_{i_{l}}\in \mathbb {F} _{q}.}

La combinación de recorrer el bucleil{\displaystyle u_{i_{l}}}los valores están representados por un conjunto conqdo1{\displaystyle q^{c-1}}elementos de una secuencia ded1{\displaystyle d-1}números enteros, es decir, i1××ido1={(i1,,ido1)}{\displaystyle u_{i_{1}}\times \cdots \times u_{i_{c-1}}=\{(i_{1},\ldots ,i_{c-1})\}}, dónde 0ilq1{\displaystyle 0\leq i_{l}\leq q-1}. Sin pérdida de generalidad , la combinación es tal que id1{\displaystyle i_{d-1}}ciclos cadaq{\displaystyle q}veces,id2{\displaystyle i_{d-2}}ciclos cadaq2{\displaystyle q^{2}}tiempos hasta i1{\displaystyle i_{1}}ciclos solo una vez. Fórmulas que calculan los índices de muestra y, por lo tanto, los grupos correspondientes, para valores fijosa{\displaystyle a}yb{\displaystyle b}, son dados por

i=l=1do1 qd1l ilvi=l=1do1 al il+b(calculado en Fq)incógnitaqi+vi=(i,vi){\displaystyle {\begin{aligned}u_{i}&=\sum _{l=1}^{c-1}~q^{d-1-l}~i_{l}\\v_{u_{i}}&=\sum _{l=1}^{c-1}~a^{l}~i_{l}+b\quad ({\text{computed in }}\mathbb {F} _{q})\\x_{qu_{i}+v_{u_{i}}}&=(u_{i},v_{u_{i}})\end{aligned}}}

Los cálculos enFq{\displaystyle \mathbb {F} _{q}}se puede implementar con bibliotecas de software disponibles públicamente para campos finitos, cuandoq{\displaystyle q}es un poder primordial. Cuandoq{\displaystyle q}es un número primo entonces los cálculos enFq{\displaystyle \mathbb {F} _{q}}simplificar a aritmética modular, es decir,vi=(l=1do1alil+b) mod q{\displaystyle v_{u_{i}}=(\sum _{l=1}^{c-1}a^{l}i_{l}+b)~{\text{mod}}~q}Un ejemplo de cómo generar un pool{\displaystyle \ell }cuando a=1,b=0,do=2{\displaystyle a=1,b=0,c=2}En la tabla siguiente se muestra la información, mientras que la selección de muestras correspondiente se muestra en la figura superior.

Este método utilizaq(do1)(d+1){\displaystyle q(c-1)(d+1)}pruebas para identificar con exactitud hastad{\displaystyle d}aspectos positivos entrenorte=qdo{\displaystyle n=q^{c}}muestras. Debido a esto, PP es particularmente eficaz para tamaños de muestra grandes, ya que el número de pruebas crece solo linealmente con respecto ado{\displaystyle c}mientras que las muestras crecen exponencialmente con este parámetro. Sin embargo, PP también puede ser eficaz para tamaños de muestra pequeños. [ 26 ]

Ejemplos de aplicaciones

La generalidad de la teoría de las pruebas grupales le confiere numerosas aplicaciones, entre las que se incluyen la detección de clones, la localización de cortocircuitos eléctricos; [ 8 ] redes informáticas de alta velocidad; [ 27 ] exámenes médicos, búsqueda de cantidades, estadística; [ 20 ] aprendizaje automático, secuenciación de ADN; [ 28 ] criptografía; [ 29 ] [ 30 ] y análisis forense de datos. [ 31 ] Esta sección ofrece una breve descripción general de una pequeña selección de estas aplicaciones.

Canales de acceso múltiple

Ilustración de un canal de acceso múltiple que muestra un mensaje exitoso y una colisión de mensajes.

Un canal de acceso múltiple es un canal de comunicación que conecta a muchos usuarios simultáneamente. Cada usuario puede escuchar y transmitir en el canal, pero si más de un usuario transmite al mismo tiempo, las señales colisionan y se reducen a ruido ininteligible. Los canales de acceso múltiple son importantes para diversas aplicaciones del mundo real, especialmente para redes informáticas inalámbricas y redes telefónicas. [ 32 ]

Un problema importante con los canales de acceso múltiple es cómo asignar tiempos de transmisión a los usuarios para que sus mensajes no colisionen. Un método simple es dar a cada usuario su propio intervalo de tiempo para transmitir, lo que requierenorte{\displaystyle n}ranuras. (Esto se denomina multiplexación por división de tiempo o TDM). Sin embargo, esto es muy ineficiente, ya que asignará ranuras de transmisión a usuarios que quizás no tengan un mensaje, y generalmente se asume que solo unos pocos usuarios querrán transmitir en un momento dado; de lo contrario, un canal de acceso múltiple no sería práctico en primer lugar.

En el contexto de las pruebas grupales, este problema generalmente se aborda dividiendo el tiempo en "épocas" de la siguiente manera. [ 8 ] Un usuario se considera "activo" si tiene un mensaje al comienzo de una época. (Si se genera un mensaje durante una época, el usuario solo se vuelve activo al comienzo de la siguiente). Una época termina cuando todos los usuarios activos han transmitido con éxito su mensaje. El problema consiste entonces en encontrar a todos los usuarios activos en una época determinada y programar un tiempo para que transmitan (si aún no lo han hecho con éxito). Aquí, una prueba en un conjunto de usuarios corresponde a aquellos usuarios que intentan una transmisión. Los resultados de la prueba son el número de usuarios que intentaron transmitir,0,1,{\displaystyle 0,1,}y2+{\displaystyle 2^{+}}, que corresponden respectivamente a ningún usuario activo, exactamente un usuario activo (mensaje exitoso) o más de un usuario activo (colisión de mensajes). Por lo tanto, utilizando un algoritmo de prueba de grupo adaptativo con resultados{0,1,2+}{\displaystyle \{0,1,2^{+}\}}Se puede determinar qué usuarios desean transmitir en ese momento. Luego, a cualquier usuario que aún no haya realizado una transmisión exitosa se le puede asignar un espacio para transmitir, evitando así asignar innecesariamente tiempos a usuarios inactivos.

Aprendizaje automático y detección comprimida

El aprendizaje automático es un campo de la informática con numerosas aplicaciones de software, como la clasificación de ADN, la detección de fraudes y la publicidad dirigida . Uno de los principales subcampos del aprendizaje automático es el problema del "aprendizaje mediante ejemplos", donde la tarea consiste en aproximar una función desconocida a partir de su valor en varios puntos específicos. [ 8 ] Como se describe en esta sección, este problema de aprendizaje de funciones puede abordarse mediante un enfoque de pruebas grupales.

En una versión simple del problema, hay alguna función desconocida,F:{0,1}norte{0,1}{\displaystyle f:\{0,1\}^{N}\to \{0,1\}}dóndeF(incógnita)=aincógnita{\displaystyle f({\textbf {x}})={\textbf {a}}\cdot {\textbf {x}}}, ya{0,1}norte{\displaystyle {\textbf {a}}\in \{0,1\}^{N}}(utilizando aritmética lógica: la suma es la operación lógica OR y la multiplicación es la operación lógica AND). Aquía{\displaystyle {\textbf {a}}}es 'd{\displaystyle d}escaso', lo que significa que como máximodnorte{\displaystyle d\ll N}de sus entradas son1{\displaystyle 1}. El objetivo es construir una aproximación aF{\displaystyle f}usandot{\displaystyle t}evaluaciones de puntos, dondet{\displaystyle t}es lo más pequeño posible. [ 4 ] (Recuperación exactaF{\displaystyle f}corresponde a algoritmos de error cero, mientras queF{\displaystyle f}(Se aproxima mediante algoritmos que tienen una probabilidad de error distinta de cero).

En este problema, recuperarF{\displaystyle f}es equivalente a encontrara{\displaystyle {\textbf {a}}}. Además,F(pag)=1{\displaystyle f({\textbf {p}})=1}si y solo si existe algún índice,norte{\displaystyle n}, dóndeanorte=pagnorte=1{\displaystyle {\textbf {a}}_{n}={\textbf {p}}_{n}=1}. Por lo tanto, este problema es análogo a un problema de prueba grupal cond{\displaystyle d}defectuosos ynorte{\displaystyle n}artículos totales. Las entradas dea{\displaystyle {\textbf {a}}}son los artículos, que son defectuosos si están1{\displaystyle 1},pag{\displaystyle {\textbf {p}}}especifica una prueba, y una prueba es positiva si y solo siF(pag)=1{\displaystyle f({\textbf {p}})=1}. [ 4 ]

En realidad, a menudo uno estará interesado en funciones más complicadas, como por ejemplo:F:donortedo{\displaystyle f:\mathbb {C} ^{N}\to \mathbb {C} }, de nuevo dondeF(incógnita)=aincógnita{\displaystyle f({\textbf {x}})={\textbf {a}}\cdot {\textbf {x}}}. La detección comprimida , que está estrechamente relacionada con las pruebas grupales, puede utilizarse para resolver este problema. [ 4 ]

En la detección comprimida, el objetivo es reconstruir una señal,vdonorte{\displaystyle {\textbf {v}}\in \mathbb {C} ^{N}}, tomando una serie de mediciones. Estas mediciones se modelan como el producto escalar dev{\displaystyle {\textbf {v}}}con un vector elegido. [ h ] El objetivo es utilizar un número pequeño de mediciones, aunque esto normalmente no es posible a menos que se asuma algo sobre la señal. Una de esas suposiciones (que es común [ 35 ] [ 36 ] ) es que solo un pequeño número de entradas dev{\displaystyle {\textbf {v}}}son significativas , lo que significa que tienen una gran magnitud. Dado que las mediciones son productos escalares dev{\displaystyle {\textbf {v}}}, la ecuaciónMETROv=q{\displaystyle M{\textbf {v}}={\textbf {q}}}sostiene, dondeMETRO{\displaystyle M}es unt×norte{\displaystyle t\times N}matriz que describe el conjunto de mediciones que se han elegido yq{\displaystyle \mathbf {q} }es el conjunto de resultados de medición. Esta construcción muestra que la detección comprimida es una especie de prueba grupal "continua".

La principal dificultad en la detección comprimida radica en identificar qué entradas son significativas. [ 35 ] Una vez hecho esto, existen diversos métodos para estimar los valores reales de las entradas. [ 37 ] Esta tarea de identificación puede abordarse mediante una sencilla aplicación de pruebas de grupo. En este caso, una prueba de grupo produce un número complejo : la suma de las entradas analizadas. El resultado de una prueba se considera positivo si produce un número complejo de gran magnitud, lo que, bajo el supuesto de que las entradas significativas son escasas, indica que al menos una entrada significativa está presente en la prueba.

Existen construcciones deterministas explícitas para este tipo de algoritmo de búsqueda combinatoria , que requierend2(registro2registro2norte)O(1){\displaystyle d2^{(\log _{2}\log _{2}N)^{O(1)}}}mediciones. [ 38 ] Sin embargo, al igual que con las pruebas grupales, estas son subóptimas y las construcciones aleatorias (como COMP) a menudo pueden recuperarF{\displaystyle f}sublinealmente ennorte{\displaystyle N}. [ 37 ]

Diseño de ensayo multiplex para pruebas de COVID-19

Durante una pandemia como el brote de COVID-19 en 2020, los ensayos de detección de virus a veces se realizan utilizando diseños de pruebas grupales no adaptativos. [ 39 ] [ 40 ] [ 41 ] Un ejemplo fue proporcionado por el proyecto Origami Assays, que publicó diseños de pruebas grupales de código abierto para ejecutarse en una placa de laboratorio estándar de 96 pocillos. [ 42 ]

Plantilla de papel de ensayo de origami para diseño de pruebas grupales

En un entorno de laboratorio, uno de los desafíos de las pruebas grupales es que la preparación de las mezclas puede ser laboriosa y difícil de realizar con precisión a mano. Los ensayos de origami proporcionaron una solución a este problema de preparación al ofrecer plantillas de papel que guían al técnico sobre cómo distribuir las muestras de los pacientes en los pocillos de prueba. [ 43 ]

Utilizando el diseño de prueba para grupos grandes (XL3), fue posible analizar 1120 muestras de pacientes en 94 pocillos de ensayo. Si la tasa de verdaderos positivos era suficientemente baja, no se requerían pruebas adicionales.

Análisis forense de datos

La informática forense es un campo dedicado a encontrar métodos para recopilar evidencia digital de un delito. Estos delitos suelen implicar que un adversario modifique los datos, documentos o bases de datos de una víctima, como por ejemplo la alteración de registros fiscales, un virus que oculta su presencia o un ladrón de identidad que modifica datos personales. [ 31 ]

Una herramienta común en el análisis forense de datos es la función hash criptográfica unidireccional . Esta función toma los datos y, mediante un procedimiento difícil de revertir, produce un número único llamado hash. [ i ] Los hashes, que suelen ser mucho más cortos que los datos, permiten comprobar si estos han sido modificados sin tener que almacenar innecesariamente copias completas de la información: el hash de los datos actuales se puede comparar con un hash anterior para determinar si se han producido cambios. Una desventaja de este método es que, si bien es fácil saber si los datos han sido modificados, no hay forma de determinar cómo: es decir, es imposible recuperar qué parte de los datos ha cambiado. [ 31 ]

Una forma de sortear esta limitación es almacenar más hashes —ahora de subconjuntos de la estructura de datos— para acotar la ubicación del ataque. Sin embargo, para encontrar la ubicación exacta del ataque con un enfoque ingenuo, sería necesario almacenar un hash para cada dato de la estructura, lo que anularía el propósito de los hashes. (Sería mejor almacenar una copia normal de los datos). Las pruebas de grupo pueden utilizarse para reducir drásticamente el número de hashes que deben almacenarse. Una prueba consiste en una comparación entre los hashes almacenados y los actuales, que es positiva cuando hay una discrepancia. Esto indica que al menos un dato editado (que se considera defectuoso en este modelo) está contenido en el grupo que generó el hash actual. [ 31 ]

De hecho, la cantidad de hashes necesarios es tan baja que, junto con la matriz de prueba a la que hacen referencia, pueden incluso almacenarse dentro de la estructura organizativa de los propios datos. Esto significa que, en lo que respecta a la memoria, la prueba puede realizarse "gratis". (Esto es cierto con la excepción de una clave maestra/contraseña que se utiliza para determinar secretamente la función hash). [ 31 ]

Notas

  1. El problema original que estudió Dorfman era de esta naturaleza (aunque no lo tuvo en cuenta), ya que, en la práctica, solo se podía agrupar una cierta cantidad de sueros sanguíneos antes de que el procedimiento de prueba se volviera poco fiable. Esta fue la razón principal por la que el procedimiento de Dorfman no se aplicó en aquel momento. [ 8 ]
  2. Sin embargo, como suele ocurrir en matemáticas, las pruebas grupales se han reinventado varias veces desde entonces, a menudo en el contexto de aplicaciones. Por ejemplo, Hayes ideó de forma independiente la realización de consultas a grupos de usuarios en el contexto de protocolos de comunicación de acceso múltiple en 1978. [ 11 ]
  3. A esto se le conoce a veces como la conjetura de Hu-Hwang-Wang.
  4. El número de pruebas,t{\displaystyle t}debe escalar comot=O(d2registrodnorte){\displaystyle t=O\left(d^{2}\log _{d}n\right)}para diseños deterministas, en comparación cont=O(dregistro2norte){\displaystyle t=O(d\log _{2}n)}para diseños que permiten probabilidades de error arbitrariamente pequeñas (comod{\displaystyle d\to \infty }ynorte{\displaystyle n\to \infty }). [ 4 ]
  5. Es importante distinguir entre el resultado falso de una prueba y el fallo del procedimiento de prueba grupal en su conjunto. Es posible cometer un error sin realizar pruebas incorrectas y, al mismo tiempo, no cometerlo aunque se realicen algunas. La mayoría de los algoritmos combinatorios modernos presentan una probabilidad de error distinta de cero (incluso sin pruebas erróneas), ya que esto reduce significativamente el número de pruebas necesarias.
  6. De hecho, es posible hacerlo mucho mejor. Por ejemplo, el de Lis{\displaystyle s}El algoritmo de -etapa proporciona una construcción explícita.tmiregistro2midregistro2(norte/d){\displaystyle t\leq {\frac {e}{\log _{2}e}}d\log _{2}{(n/d)}}.
  7. Alternativamentey{\displaystyle \mathbf {y} }puede definirse mediante la ecuacióny:=METROincógnita{\displaystyle \mathbf {y} :=M\mathbf {x} } , donde la multiplicación es AND lógico ({\displaystyle \wedge }) y la suma es OR lógico ({\displaystyle \vee }). Aquí,y{\displaystyle \mathbf {y} }tendrá una1{\displaystyle 1}en posicióni{\displaystyle i}si y solo si(METRO)i,j{\displaystyle (M)_{i,j}}yincógnitaj{\displaystyle \mathbf {x} _{j}}son ambos1{\displaystyle 1}para cualquierj{\displaystyle j}. Es decir, si y solo si al menos un artículo defectuoso estaba incluido en eli-th{\displaystyle i{\textrm {-th}}}prueba.
  8. Este tipo de medición aparece en muchas aplicaciones. Por ejemplo, ciertos tipos de cámaras digitales [ 33 ] o máquinas de resonancia magnética [ 34 ] , donde las limitaciones de tiempo requieren que solo se tome un pequeño número de mediciones.
  9. Formalmente, las funciones hash poseen una propiedad denominada resistencia a colisiones, que consiste en que la probabilidad de obtener el mismo hash a partir de diferentes entradas es muy baja para datos de un tamaño adecuado. En la práctica, la posibilidad de que dos entradas diferentes produzcan el mismo hash suele ignorarse.

Referencias

Citas

  1. Colbourn, Charles J.; Dinitz, Jeffrey H. (2007), Handbook of Combinatorial Designs (2.ª  ed.), Boca Raton: Chapman & Hall/ CRC, p. 574, Sección 46: Pooling Designs , ISBN 978-1-58488-506-1
  2. 1 2 3 Dorfman, Robert (diciembre de 1943), "La detección de miembros defectuosos de grandes poblaciones", The Annals of Mathematical Statistics , 14 (4): 436– 440, doi : 10.1214/aoms/1177731363 , JSTOR 2235930 
  3. 1 2 3 4 5 6 7 Ding-Zhu, Du; Hwang, Frank K. (2000). Pruebas de grupos combinatorios y sus aplicaciones (2.ª ed.). Singapur: World Scientific. ISBN  978-9810241070.
  4. 1 2 3 4 5 6 7 8 9 Atia, George Kamal; Saligrama, Venkatesh (marzo de 2012). "Boolean compressed sensing and noisy group testing". IEEE Transactions on Information Theory . 58 (3): 1880– 1901. arXiv : 0907.1061 . Bibcode : 2012ITIT...58.1880A . doi : 10.1109/TIT.2011.2178156 . S2CID 8946216 . 
  5. Knill, E.; Bruno, WJ; Torney, DC (1998). "Pruebas grupales no adaptativas en presencia de errores". Discrete Applied Mathematics . 88 ( 1– 3): 261– 290. doi : 10.1016/S0166-218X(98)00075-4 . MR 1658592 . Para muchas aplicaciones de cribado, es más rentable hacer muchas consultas de subconjuntos en paralelo. Esto lleva a problemas de pruebas grupales no adaptativas. 
  6. 1 2 3 4 5 6 7 8 9 10 11 Chun Lam Chan; Pak Hou Che; Jaggi, Sidharth; Saligrama, Venkatesh (1 de septiembre de 2011). "Pruebas grupales probabilísticas no adaptativas con mediciones ruidosas: límites casi óptimos con algoritmos eficientes". 49.ª Conferencia Anual de Allerton sobre Comunicación, Control y Computación . págs. 1832-1839 . arXiv : 1107.4540 . doi : 10.1109/Allerton.2011.6120391 . ISBN  978-1-4577-1817-5. S2CID 8408114 . 
  7. Hung, M.; Swallow, William H. (marzo de 1999). "Robustez de las pruebas grupales en la estimación de proporciones". Biometrics . 55 ( 1): 231– 7. doi : 10.1111/j.0006-341X.1999.00231.x . PMID 11318160. S2CID 23389365 .  
  8. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Ding-Zhu, Du; Hwang, Frank K. (1993). Pruebas de grupos combinatorios y sus aplicaciones . Singapur: World Scientific. ISBN 978-9810212933.
  9. Chen, Hong-Bin; Fu, Hung-Lin (abril de 2009). "Algoritmos no adaptativos para pruebas de grupos de umbral" . Matemáticas Aplicadas Discretas . 157 (7): 1581– 1585. doi : 10.1016/j.dam.2008.06.003 .
  10. De Bonis, Annalisa (20 de julio de 2007). "Nuevas estructuras combinatorias con aplicaciones a pruebas de grupo eficientes con inhibidores". Journal of Combinatorial Optimization . 15 (1): 77– 94. doi : 10.1007/s10878-007-9085-1 . S2CID 207188798 . 
  11. Hayes, J. (agosto de 1978). "Una técnica adaptativa para la distribución local". IEEE Transactions on Communications . 26 (8): 1178– 86. Bibcode : 1978ITCom..26.1178H . doi : 10.1109/TCOM.1978.1094204 .
  12. Samuels, Stephen (1978). "La solución exacta al problema de las pruebas grupales en dos etapas" . Technometrics . 20 (4): 497– 500. doi : 10.1080/00401706.1978.10489706 .
  13. Sterrett, Andrew (diciembre de 1957). "Sobre la detección de miembros defectuosos de grandes poblaciones" . The Annals of Mathematical Statistics . 28 (4): 1033– 6. doi : 10.1214/aoms/1177706807 .
  14. Sobel, Milton; Groll, Phyllis A. (septiembre de 1959). "Pruebas grupales para eliminar eficientemente todos los defectuosos en una muestra binomial". Bell System Technical Journal . 38 (5): 1179– 1252. Bibcode : 1959BSTJ...38.1179S . doi : 10.1002/j.1538-7305.1959.tb03914.x .
  15. Ungar, Peter (febrero de 1960). "Puntos de corte en pruebas grupales" . Communications on Pure and Applied Mathematics . 13 (1): 49– 54. doi : 10.1002/cpa.3160130105 .
  16. 1 2 Li, Chou Hsiung (junio de 1962). "Un método secuencial para seleccionar variables experimentales". Journal of the American Statistical Association . 57 (298): 455– 477. doi : 10.1080/01621459.1962.10480672 .
  17. Katona, Gyula OH (1973). "Un estudio de la teoría combinatoria". Problemas de búsqueda combinatoria . North-Holland. págs. 285–308 . ISBN  978-0-7204-2262-7.
  18. 1 2 3 4 Hwang, Frank K. (septiembre de 1972). "Un método para detectar todos los miembros defectuosos en una población mediante pruebas grupales". Journal of the American Statistical Association . 67 (339): 605– 608. doi : 10.2307/2284447 . JSTOR 2284447 . 
  19. Allemann, Andreas (2013). "Un algoritmo eficiente para pruebas de grupos combinatorios". Teoría de la información, combinatoria y teoría de la búsqueda . Lecture Notes in Computer Science. Vol. 7777. pp. 569–596 . doi : 10.1007/978-3-642-36899-8_29 . ISBN   978-3-642-36898-1.
  20. 1 2 Hu, MC; Hwang, FK; Wang, Ju Kwei (junio de 1981). "Un problema de frontera para pruebas grupales". SIAM Journal on Algebraic and Discrete Methods . 2 (2): 81– 87. doi : 10.1137/0602011 .
  21. Leu, Ming-Guang (28 de octubre de 2008). "Una nota sobre la conjetura de Hu-Hwang-Wang para pruebas grupales" . The ANZIAM Journal . 49 (4): 561. doi : 10.1017/S1446181108000175 .
  22. Riccio, Laura; Colbourn, Charles J. (1 de enero de 2000). "Límites más precisos en las pruebas grupales adaptativas" . Taiwanese Journal of Mathematics . 4 (4): 669– 673. doi : 10.11650/twjm/1500407300 .
  23. 1 2 3 4 Aldridge, Matthew; Baldassini, Leonardo; Johnson, Oliver (junio de 2014). "Algoritmos de pruebas grupales: límites y simulaciones". IEEE Transactions on Information Theory . 60 (6): 3671– 3687. arXiv : 1306.6438 . Bibcode : 2014ITIT...60.3671A . doi : 10.1109/TIT.2014.2314472 . S2CID 8885619 . 
  24. Baldassini, L.; Johnson, O.; Aldridge, M. (1 de julio de 2013), "La capacidad de las pruebas grupales adaptativas", 2013 IEEE International Symposium on Information Theory , pp. 2676–2680 , arXiv : 1301.7023 , CiteSeerX 10.1.1.768.8924 , doi : 10.1109/ISIT.2013.6620712 , ISBN   978-1-4799-0446-4, S2CID 9987210 
  25. Sobel, Milton; Elashoff, RM (1975). "Pruebas grupales con un nuevo objetivo, la estimación". Biometrika . 62 (1): 181– 193. doi : 10.1093/biomet/62.1.181 . hdl : 11299/199154 ​​.
  26. 1 2 Brust, D.; Brust, JJ (enero de 2023). " Diseños de matriz eficaces para pruebas grupales de COVID-19" . BMC Bioinformatics . 24 (26): 26. doi : 10.1186/s12859-023-05145-y . PMC 9872308. PMID 36694117 .  
  27. Bar-Noy, A.; Hwang, FK; Kessler, I.; Kutten, S. (1 de mayo de 1992). «Un nuevo algoritmo competitivo para pruebas grupales». [ Actas ] IEEE INFOCOM '92: Conferencia sobre Comunicaciones Informáticas . Vol. 2. págs. 786–793 . doi : 10.1109/INFCOM.1992.263516 . ISBN   978-0-7803-0602-8. S2CID 16131063 . 
  28. Damaschke, Peter (2000). "Aprendizaje adaptativo versus no adaptativo eficiente en atributos" . Machine Learning . 41 (2): 197– 215. doi : 10.1023/A:1007616604496 .
  29. Stinson, DR; van Trung, Tran; Wei, R (mayo de 2000). "Códigos seguros a prueba de marcos, patrones de distribución de claves, algoritmos de prueba de grupos y estructuras relacionadas". Journal of Statistical Planning and Inference . 86 (2): 595– 617. CiteSeerX 10.1.1.54.6212 . doi : 10.1016/S0378-3758(99)00131-7 . 
  30. Colbourn, CJ; Dinitz, JH; Stinson, DR (1999). "Comunicaciones, criptografía y redes" . Surveys in Combinatorics . 3 (267): 37– 41. doi : 10.1007/BF01609873 . S2CID 10128581 . 
  31. 1 2 3 4 5 Goodrich, Michael T.; Atallah, Mikhail J.; Tamassia, Roberto (2005). "Indexación de información para análisis forense de datos". Criptografía aplicada y seguridad de redes . Notas de clase en informática. Vol. 3531. págs. 206–221 . CiteSeerX 10.1.1.158.6036 . doi : 10.1007/11496137_15 . ISBN    978-3-540-26223-7.
  32. Chlebus, BS (2001). "Comunicación aleatoria en redes de radio" . En Pardalos, PM; Rajasekaran, S.; Reif, J.; Rolim, JDP (eds.). Manual de computación aleatoria . Kluwer Academic. pp. 401–456 . ISBN  978-0-7923-6957-8.
  33. Takhar, D.; Laska, JN; Wakin, MB; Duarte, MF; Baron, D.; Sarvotham, S.; Kelly, KF; Baraniuk, RG (febrero de 2006). Bouman, Charles A.; Miller, Eric L.; Pollak, Ilya (eds.). "Una nueva arquitectura de cámara de imágenes compresivas que utiliza compresión en el dominio óptico". Electronic Imaging . Computational Imaging IV. 6065 : 606509–606509–10. Bibcode : 2006SPIE.6065...43T . CiteSeerX 10.1.1.114.7872 . doi : 10.1117/12.659602 . S2CID 7513433 .  
  34. Candès, EJ (2014). "Matemáticas de la escasez (y algunas otras cosas)". Actas del Congreso Internacional de Matemáticos. Seúl, Corea del Sur .
  35. 1 2 Gilbert, AC; Iwen, MA; Strauss, MJ (octubre de 2008). "Pruebas de grupo y recuperación de señales dispersas". 42.ª Conferencia de Asilomar sobre Señales, Sistemas y Computadoras . Instituto de Ingenieros Eléctricos y Electrónicos. págs. 1059–63 . doi : 10.1109/ACSSC.2008.5074574 . ISBN  978-1-4244-2940-0.
  36. Wright, SJ; Nowak, RD; Figueiredo, MAT (julio de 2009). "Reconstrucción dispersa mediante aproximación separable". IEEE Transactions on Signal Processing . 57 (7): 2479– 2493. Bibcode : 2009ITSP...57.2479W . CiteSeerX 10.1.1.142.749 . doi : 10.1109/TSP.2009.2016892 . S2CID 7399917 .  
  37. 1 2 Berinde, R.; Gilbert, AC; Indyk, P.; Karloff, H.; Strauss, MJ (septiembre de 2008). "Combinando geometría y combinatoria: un enfoque unificado para la recuperación de señales dispersas". 46.ª Conferencia Anual de Allerton sobre Comunicación, Control y Computación de 2008. págs. 798–805 . arXiv : 0804.4666 . doi : 10.1109/ALLERTON.2008.4797639 . ISBN  978-1-4244-2925-7. S2CID 8301134 . 
  38. Indyk, Piotr (1 de enero de 2008). "Construcciones explícitas para la detección comprimida de señales dispersas". Actas del decimonoveno simposio anual ACM-SIAM sobre algoritmos discretos : 30–33 .
  39. Austin, David. "Columna destacada de la AMS: estrategias de agrupación para las pruebas de COVID-19" . Sociedad Matemática Estadounidense . Consultado el 3 de octubre de 2020 .
  40. ^ Prasanna, Dheeraj. "Agrupación de tapices" . tapestry-pooling.herokuapp.com . Consultado el 3 de octubre de 2020 .
  41. Chiani, M.; Liva, G.; Paolini, E. (febrero de 2022), "Protocolos de pruebas grupales de identificación-detección para COVID-19 en zonas de alta prevalencia", Scientific Reports , 12 (1), Springer Nature: 3250, arXiv : 2104.11305 , Bibcode : 2022NatSR..12.3250C , doi : 10.1038/s41598-022-07205-4 , PMC 8885674 , PMID 35228579 , S2CID 233387831   
  42. "Ensayos de origami" . Ensayos de origami. 2 de abril de 2020. Consultado el 7 de abril de 2020 .
  43. "Ensayos de origami" . Ensayos de origami. 2 de abril de 2020. Consultado el 7 de abril de 2020 .

Referencias generales

  • Ding-Zhu, Du; Hwang, Frank K. (2000). Pruebas de grupos combinatorios y sus aplicaciones (2.ª  ed.). Singapur: World Scientific. ISBN 978-9810241070.
  • Curso de Atri Rudra sobre Códigos Correctores de Errores: Combinatoria, Algoritmos y Aplicaciones (Primavera de 2007), Lección 7 .
  • Curso de Atri Rudra sobre Códigos Correctores de Errores: Combinatoria, Algoritmos y Aplicaciones (Primavera de 2010), Lecciones 10 , 11 , 28 , 29
  • Du, D.; Hwang, F. (2006). Diseños de agrupamiento y pruebas grupales no adaptativas . World Scientific. ISBN 9789814477864.
  • Aldridge, M.; Johnson, O.; Scarlett, J. (2019). "Pruebas grupales: una perspectiva de la teoría de la información" (PDF) . Fundamentos y tendencias en la teoría de la comunicación y la información . 15 ( 3–4 ): 196–392 . arXiv : 1902.06002 . doi : 10.1561/0100000099 . S2CID 62841593 . 
  • Porat, E.; Rothschild, A. (2011). "Esquemas explícitos de prueba de grupos combinatorios no adaptativos". IEEE Transactions on Information Theory . 57 (12): 7982– 89. arXiv : 0712.3876 . Bibcode : 2011ITIT...57.7982P . doi : 10.1109/TIT.2011.2163296 . S2CID 8815474 . 
  • Kagan, Eugene; Ben-gal, Irad (2014), "Un algoritmo de pruebas grupales con aprendizaje informacional en línea", IIE Transactions , 46 (2): 164–184 , doi : 10.1080/0740817X.2013.803639 , ISSN 0740-817X , S2CID 18588494  

Véase también