EJEMPLO DE PRUEBA ESTADISTICA KOLMOGORNOV SMIRNOV

 EJEMPLO DE PRUEBA ESTADISTICA  KOLMOGORNOV SMIRNOV

import numpy as np

from scipy.stats import kstest, norm

import matplotlib.pyplot as plt

np.random.seed(42)

sample = np.random.normal(0, 1, 100)  # Generar una muestra de 100 puntos de una distribución normal

# Calcular ECDF

ecdf = np.linspace(0, 1, len(sample))

# Graficar la ECDF

plt.plot(np.sort(sample), ecdf)

plt.plot(norm.cdf(np.sort(sample)), 'r--', label='CDF teórica (normal)')

plt.title('ECDF vs CDF Teórica (Normal)')

plt.xlabel('Valores de la Muestra')

plt.ylabel('Probabilidad acumulada')

plt.legend()

plt.show()

# Realizar la prueba KS con una distribución normal estándar (media=0, desviación=1)

kstest_result = kstest(sample, 'norm')

print('Estadístico KS:', kstest_result.statistic)

print('Valor p:', kstest_result.pvalue)