Articulo de referencia

Muestreo local de casos y controles

En aprendizaje automático , el muestreo local de casos y controles [ 1 ] es un algoritmo que se utiliza para reducir la complejidad del entrenamiento de un clasificador de regre...

En aprendizaje automático , el muestreo local de casos y controles [ 1 ] es un algoritmo que se utiliza para reducir la complejidad del entrenamiento de un clasificador de regresión logística . El algoritmo reduce la complejidad del entrenamiento seleccionando una pequeña submuestra del conjunto de datos original para el entrenamiento. Presupone la disponibilidad de una estimación piloto (poco fiable) de los parámetros. A continuación, realiza una única pasada sobre todo el conjunto de datos utilizando la estimación piloto para identificar las muestras más "sorprendentes". En la práctica, la estimación piloto puede provenir del conocimiento previo o del entrenamiento utilizando una submuestra del conjunto de datos. El algoritmo es más eficaz cuando el conjunto de datos subyacente está desequilibrado. Aprovecha las estructuras de los conjuntos de datos desequilibrados condicionalmente de forma más eficiente que otros métodos, como el muestreo de casos y controles y el muestreo ponderado de casos y controles.

Conjuntos de datos desequilibrados

En clasificación , un conjunto de datos es un conjunto de N puntos de datos.(incógnitai,yi)i=1norte{\displaystyle (x_{i},y_{i})_{i=1}^{N}}, dóndeincógnitaiRd{\displaystyle x_{i}\in \mathbb {R} ^{d}}es un vector de características ,yi{0,1}{\displaystyle y_{i}\in \{0,1\}}es una etiqueta. Intuitivamente, un conjunto de datos está desequilibrado cuando ciertos patrones estadísticos importantes son poco frecuentes. La falta de observaciones de ciertos patrones no siempre implica su irrelevancia. Por ejemplo, en estudios médicos de enfermedades raras, el pequeño número de pacientes infectados (casos) proporciona la información más valiosa para el diagnóstico y el tratamiento.

Formalmente, un conjunto de datos desequilibrado presenta una o más de las siguientes propiedades:

  • Desequilibrio marginal . Un conjunto de datos está marginalmente desequilibrado si una clase es rara en comparación con la otra clase. En otras palabras,PAG(Y=1)0{\displaystyle \mathbb {P} (Y=1)\approx 0}.
  • Desequilibrio condicional . Un conjunto de datos está desequilibrado condicionalmente cuando es fácil predecir las etiquetas correctas en la mayoría de los casos. Por ejemplo, siincógnita{0,1}{\displaystyle X\in \{0,1\}}, el conjunto de datos está desequilibrado condicionalmente siPAG(Y=1incógnita=0)0{\displaystyle \mathbb {P} (Y=1\mid X=0)\approx 0}yPAG(Y=1incógnita=1)1{\displaystyle \mathbb {P} (Y=1\mid X=1)\approx 1}.

Esquema del algoritmo

En la regresión logística, dado el modeloθ=(α,β){\displaystyle \theta =(\alpha ,\beta )}, la predicción se realiza de acuerdo conPAG(Y=1incógnita;θ)=pag~θ(incógnita)=exp(α+βTincógnita)1+exp(α+βTincógnita){\displaystyle \mathbb {P} (Y=1\mid X;\theta )={\tilde {p}}_{\theta }(x)={\frac {\exp(\alpha +\beta ^{T}x)}{1+\exp(\alpha +\beta ^{T}x)}}}El algoritmo de muestreo de casos y controles locales presupone la disponibilidad de un modelo piloto.θ~=(α~,β~){\displaystyle {\tilde {\theta }}=({\tilde {\alpha }},{\tilde {\beta }})}Dado el modelo piloto, el algoritmo realiza una única pasada sobre todo el conjunto de datos para seleccionar el subconjunto de muestras que se incluirán en el entrenamiento del modelo de regresión logística. Para una muestra(incógnita,y){\displaystyle (x,y)}, defina la probabilidad de aceptación comoa(incógnita,y)=|ypag~θ~(incógnita)|{\displaystyle a(x,y)=|y-{\tilde {p}}_{\tilde {\theta }}(x)|}El algoritmo procede de la siguiente manera:

  1. Generar independientesziBernoulli(a(incógnitai,yi)){\displaystyle z_{i}\sim {\text{Bernoulli}}(a(x_{i},y_{i}))}parai{1,,norte}{\displaystyle i\in \{1,\ldots ,N\}}.
  2. Ajustar un modelo de regresión logística a la submuestra.S={(incógnitai,yi):zi=1}{\displaystyle S=\{(x_{i},y_{i}):z_{i}=1\}}obteniendo las estimaciones no ajustadasθ^S=(α^S,β^S){\displaystyle {\hat {\theta }}_{S}=({\hat {\alpha }}_{S},{\hat {\beta }}_{S})}.
  3. El modelo de salida esθ^=(α^,β^){\displaystyle {\hat {\theta }}=({\hat {\alpha }},{\hat {\beta }})}, dóndeα^α^S+α~{\displaystyle {\hat {\alpha }}\leftarrow {\hat {\alpha }}_{S}+{\tilde {\alpha }}}yβ^β^S+β~{\displaystyle {\hat {\beta }}\leftarrow {\hat {\beta }}_{S}+{\tilde {\beta }}}.

El algoritmo puede entenderse como la selección de muestras que sorprenden al modelo piloto. Intuitivamente, estas muestras están más cerca del límite de decisión del clasificador y, por lo tanto, son más informativas.

Obtención del modelo piloto

En la práctica, en los casos en que se dispone de un modelo piloto de forma natural, el algoritmo puede aplicarse directamente para reducir la complejidad del entrenamiento. En los casos en que no existe un piloto natural, se puede utilizar una estimación mediante una submuestra seleccionada a través de otra técnica de muestreo. En el artículo original que describe el algoritmo, los autores proponen utilizar un muestreo ponderado de casos y controles con la mitad del presupuesto de muestreo asignado. Por ejemplo, si el objetivo es utilizar una submuestra con tamañonorte=1000{\displaystyle N=1000}Primero, estima un modelo.θ~{\displaystyle {\tilde {\theta }}}usandonorteh=500{\displaystyle N_{h}=500}muestras del muestreo ponderado de casos y controles, luego recoja otranorteh=500{\displaystyle N_{h}=500}muestras mediante muestreo local de casos y controles.

Tamaño de muestra mayor o menor

Es posible controlar el tamaño de la muestra multiplicando la probabilidad de aceptación por una constante.do{\displaystyle c}Para un tamaño de muestra mayor, seleccionedo>1{\displaystyle c>1}y ajustar la probabilidad de aceptación amin(doa(incógnitai,yi),1){\displaystyle \min(ca(x_{i},y_{i}),1)}Para un tamaño de muestra menor, se aplica la misma estrategia. En los casos en que se desea un número preciso de muestras, un método alternativo conveniente consiste en realizar un submuestreo uniforme a partir de una submuestra mayor seleccionada mediante muestreo local de casos y controles.

Propiedades

El algoritmo tiene las siguientes propiedades. Cuando el piloto es consistente , las estimaciones que utilizan las muestras del muestreo de casos y controles local son consistentes incluso bajo una especificación incorrecta del modelo . Si el modelo es correcto, entonces el algoritmo tiene exactamente el doble de la varianza asintótica de la regresión logística en el conjunto de datos completo. Para un tamaño de muestra mayor condo>1{\displaystyle c>1}, el factor 2 se mejora a1+1do{\displaystyle 1+{\frac {1}{c}}}.

Referencias

  1. Fithian, William; Hastie, Trevor (2014). "Muestreo local de casos y controles: submuestreo eficiente en conjuntos de datos desequilibrados" . The Annals of Statistics . 42 (5): 1693– 1724. arXiv : 1306.3706 . doi : 10.1214/14-aos1220 . PMC 4258397. PMID 25492979 .