Calculadora de IC por bootstrap

Estima IC por remuestreo.

Ejemplos:

Usa directamente los cuantiles de la distribución bootstrap

Use la misma semilla para resultados reproducibles

Ingrese datos y haga clic en 'Ejecutar Bootstrap' para calcular el IC

Usa intervalos de confianza bootstrap cuando no puedas asumir normalidad, tengas muestras pequeñas, o necesites intervalos para estadísticos complejos como medianas, razones o percentiles. El remuestreo bootstrap no hace supuestos distribucionales.

¿Qué es el bootstrapping?

El bootstrapping crea muchos conjuntos de datos remuestreados extrayendo con reemplazo de tus datos originales, luego calcula tu estadístico en cada remuestra:

  1. Muestra original: n observaciones
  2. Extrae n observaciones con reemplazo → Muestra bootstrap 1
  3. Calcula el estadístico en la muestra bootstrap 1
  4. Repite B veces (típicamente 1000-10000)
  5. Usa la distribución de estadísticos bootstrap para inferencia

Métodos de IC bootstrap

Método de percentiles

El enfoque más simple: usa percentiles de la distribución bootstrap.

IC95%=[B0,025,B0,975]IC_{95\%} = [B_{0,025}, B_{0,975}]

Donde BpB_{p} es el percentil p-ésimo de los estadísticos bootstrap.

Pros: Simple, intuitivo Contras: Puede estar sesgado, pobre cobertura para distribuciones sesgadas

Método básico (percentil inverso)

IC95%=[2θ^−B0,975,2θ^−B0,025]IC_{95\%} = [2\hat{\theta} - B_{0,975}, 2\hat{\theta} - B_{0,025}]

Donde θ^\hat{\theta} es la estimación de la muestra original.

Refleja la distribución bootstrap alrededor de la estimación original.

Corregido por sesgo y acelerado (BCa)

El estándar de oro para ICs bootstrap:

ICBCa=[Bα1,Bα2]IC_{BCa} = [B_{\alpha_1}, B_{\alpha_2}]

Donde: α1=Φ(z0+z0+zα/21−a(z0+zα/2))\alpha_1 = \Phi\left(z_0 + \frac{z_0 + z_{\alpha/2}}{1 - a(z_0 + z_{\alpha/2})}\right)

α2=Φ(z0+z0+z1−α/21−a(z0+z1−α/2))\alpha_2 = \Phi\left(z_0 + \frac{z_0 + z_{1-\alpha/2}}{1 - a(z_0 + z_{1-\alpha/2})}\right)

  • z0z_0 = corrección de sesgo (proporción de estimaciones bootstrap por debajo de la original)
  • aa = aceleración (mide la asimetría)

Pros: Mejor cobertura, maneja sesgo y asimetría Contras: Computacionalmente intensivo, requiere más muestras bootstrap

Estudentizado (Bootstrap-t)

IC=[θ^−t1−α/2∗⋅EE,θ^−tα/2∗⋅EE]IC = [\hat{\theta} - t^*_{1-\alpha/2} \cdot EE, \hat{\theta} - t^*_{\alpha/2} \cdot EE]

Usa la distribución de estadísticos estudentizados: t∗=θ^∗−θ^EE∗t^* = \frac{\hat{\theta}^* - \hat{\theta}}{EE^*}

Requiere estimar EE para cada muestra bootstrap.

Ejemplo: IC bootstrap para la mediana

Datos: 12, 15, 18, 22, 25, 28, 45, 52 (n = 8)

Mediana original: 23,5

1000 muestras bootstrap dan medianas:

  • Percentil 2,5: 17,5
  • Percentil 97,5: 38,5

IC percentil 95%: [17,5; 38,5]

Comparar con: ¡No existe fórmula simple para IC de mediana con teoría normal!

¿Cuántas muestras bootstrap?

PropósitoB recomendado
Error estándar200-500
Intervalos de confianza1000-2000
IC percentil1000
IC BCa2000-5000
Prueba de hipótesis1000-10000

Más muestras = más precisión, pero rendimientos decrecientes después de ~10.000.

Cuándo usar ICs bootstrap

Situaciones ideales

  1. Datos no normales: Distribuciones sesgadas, valores atípicos
  2. Muestras pequeñas: n < 30 donde el TLC es cuestionable
  3. Estadísticos complejos: Mediana, media recortada, razones, percentiles
  4. Sin EE de forma cerrada: Diferencias de correlación, predicciones de regresión
  5. Datos dependientes: Series de tiempo (usa bootstrap por bloques)

Cuándo NO usar

  1. Muestras muy pequeñas: n < 10 puede no capturar la variabilidad poblacional
  2. Colas pesadas con n pequeño: Bootstrap puede ser inestable
  3. Estadísticos simples con n grande: La teoría normal es adecuada y más rápida

Bootstrap para diferentes estadísticos

Media

Xˉ∗=1n∑i=1nXi∗\bar{X}^* = \frac{1}{n}\sum_{i=1}^{n} X_i^*

Bootstrap es innecesario para n grande (usa teoría normal), pero válido para n pequeño o datos no normales.

Mediana

No hay fórmula simple para EE de mediana. Bootstrap es a menudo el mejor enfoque:

X~∗=mediana(X1∗,...,Xn∗)\tilde{X}^* = \text{mediana}(X_1^*, ..., X_n^*)

Correlación

r∗=cor(X∗,Y∗)r^* = \text{cor}(X^*, Y^*)

Bootstrap maneja la distribución muestral no normal de r.

Coeficientes de regresión

Remuestrea pares (X, Y) juntos:

β^∗=(X∗TX∗)−1X∗TY∗\hat{\beta}^* = (X^{*T}X^*)^{-1}X^{*T}Y^*

Cuenta para heterocedasticidad sin asumirla.

Diferencia de medias

Remuestrea cada grupo por separado:

D∗=Xˉ1∗−Xˉ2∗D^* = \bar{X}_1^* - \bar{X}_2^*

Válido incluso con varianzas desiguales y no normalidad.

Bootstrap paramétrico vs. no paramétrico

Bootstrap no paramétrico

Remuestrea de los datos mismos. No hace supuestos distribucionales.

Bootstrap paramétrico

  1. Ajusta un modelo a los datos (por ejemplo, estima μ, σ)
  2. Genera muestras del modelo ajustado
  3. Calcula el estadístico en cada muestra generada

Usar cuando: Estás seguro de la forma distribucional pero incierto sobre los parámetros.

Bootstrap por bloques para datos dependientes

Para series de tiempo o datos espaciales, las observaciones no son independientes. El bootstrap por bloques preserva la dependencia local:

  1. Divide los datos en bloques superpuestos
  2. Remuestrea bloques (no observaciones individuales)
  3. Concatena para formar la muestra bootstrap

El tamaño del bloque requiere ajuste—muy pequeño ignora dependencia, muy grande reduce variabilidad.

Verificando resultados bootstrap

Convergencia

Aumenta B hasta que el EE del EE bootstrap se estabilice:

EEbootstrap(EE)≈EE2BEE_{bootstrap}(EE) \approx \frac{EE}{\sqrt{2B}}

Sesgo

Sesgo=θ∗ˉ−θ^\text{Sesgo} = \bar{\theta^*} - \hat{\theta}

Si el sesgo es grande relativo al EE, considera el método BCa.

Normalidad de la distribución bootstrap

Si la distribución bootstrap es altamente no normal, el método de percentiles puede tener pobre cobertura. Usa BCa.

Errores comunes

Remuestrear unidades incorrectas

Para datos agrupados, remuestrea clusters, no observaciones. Para datos pareados, remuestrea pares.

Ignorar dependencia

El bootstrap estándar asume independencia. Usa bootstrap por bloques para series de tiempo.

Muy pocas muestras bootstrap

B = 100 es muy poco para ICs. Usa al menos 1000.

Esperar milagros con n pequeño

Bootstrap no puede crear información que no está en tu muestra. Con n = 5, tu IC será amplio de todos modos.

Comparación con otros métodos

MétodoSupuestosMejor para
Teoría normalNormalidad, n grandeMedias con muestras grandes
Intervalo tNormalidad (robusto)Medias con muestra pequeña
BootstrapNinguno (IID)Estadísticos complejos, datos no normales
BayesianoEspecificación del a prioriIncorporar conocimiento previo

Reportando resultados bootstrap

Formato estándar: "El tiempo de respuesta mediano fue 23,5 ms (IC BCa bootstrap 95%: 17,5–38,5 ms, B = 5000)."

Incluir:

  • Estimación original
  • Método de IC (percentil, BCa, etc.)
  • Número de muestras bootstrap
  • Cualquier consideración especial (tamaño de bloque para series de tiempo)

Consideraciones computacionales

Consejos de eficiencia

  1. Vectoriza operaciones
  2. Usa procesamiento paralelo para B grande
  3. Almacena solo los estadísticos necesarios, no remuestras completas
  4. Considera submuestreo para conjuntos de datos enormes

Reproducibilidad

Establece semilla aleatoria antes del bootstrapping para resultados reproducibles.

Contexto histórico

Bradley Efron introdujo el bootstrap en 1979, revolucionando la inferencia estadística. El método BCa (1987) mejoró las propiedades de cobertura. Los avances en computación hicieron el bootstrap práctico, cambiando fundamentalmente cómo los estadísticos abordan problemas de inferencia complejos.

Cálculos relacionados

Resumen

Los intervalos de confianza bootstrap funcionan remuestreando tus datos con reemplazo, calculando el estadístico de interés muchas veces y usando la distribución resultante para inferencia. El método de percentiles es el más simple pero BCa se prefiere por sus mejores propiedades de cobertura. Usa bootstrap cuando la normalidad es cuestionable, las muestras son pequeñas, o para estadísticos sin errores estándar de forma cerrada. Elige al menos 1000 remuestras para ICs, 2000+ para BCa. Reporta el método usado y el número de muestras bootstrap.