Fit distribución empírica a teóricas con Scipy (Python)?

INTRODUCCIÓN: Tengo una lista de más de 30 000 valores que van de 0 a 47, por ejemplo, [0,0,0,0, .., 1,1,1,1, ..., 2,2,2,2 , ..., 47 etc.] que es la distribución continua.

PROBLEMA: Según mi distribución, me gustaría calcular el valor p (la probabilidad de ver valores mayores) para cualquier valor dado. Por ejemplo, como puede ver, el valor p para 0 se acercará a 1 y el valor p para números más altos tenderá a 0.

No sé si tengo razón, pero para determinar las probabilidades, creo que necesito ajustar mis datos a una distribución teórica que sea la más adecuada para describir mis datos. Supongo que se necesita algún tipo de prueba de bondad de ajuste para determinar el mejor modelo.

¿Hay alguna forma de implementar dicho análisis en Python (Scipy o Numpy)? ¿Podría presentar algún ejemplo?

¡Gracias

Respuestas a la pregunta(6)

Su respuesta a la pregunta