Comparar dos medias en Python
La prueba t, y la variante que no asume varianzas iguales.
import numpy as np
from scipy import stats
rng = np.random.default_rng(1)
antes = rng.normal(90, 10, size=40)
despues = rng.normal(97, 12, size=40)
welch = stats.ttest_ind(antes, despues, equal_var=False)
print(round(welch.statistic, 3), f"{welch.pvalue:.5f}")
print(welch.confidence_interval(0.95).low.round(2))
pareada = stats.ttest_rel(antes, despues)
print(round(pareada.statistic, 3), f"{pareada.pvalue:.5f}")
print(stats.ttest_1samp(despues, popmean=90).pvalue < 0.05)
Cómo funciona
ttest_indcompara muestras independientes.equal_var=Falsees la prueba de Welch, la opción por defecto más segura.ttest_reles para mediciones pareadas de los mismos sujetos.
Palabras clave y builtins usados aquí
asprintround
El intento, en números
- Líneas
- 14
- Caracteres a escribir
- 477
- Tokens
- 147
- Ritmo de tres estrellas
- 105 tpm
Al ritmo de tres estrellas de 105 tokens por minuto, este intento toma unos 84 segundos.
Paso 1 de 4 en Pruebas de hipótesis; paso 4 de 23 en Computación científica con SciPy.