En la última lección, realizaste una distribución normal multivariante utilizando la media y la matriz de covarianza de dia. Ahora, ¡responderás a preguntas de interés utilizando los resultados simulados!

Te preguntarás: ¿por qué hacemos simulaciones si tenemos datos históricos? ¿No podemos utilizar los propios datos para responder a preguntas de interés?

Esta es una gran pregunta. Las simulaciones de Montecarlo se basan en la modelización mediante distribuciones de probabilidad, que ofrecen toda la distribución de probabilidad para su inspección (un gran número de muestras), en lugar del número limitado de puntos de datos disponibles en los datos históricos.

Por ejemplo, puedes hacer preguntas como ¿cuál es el cuantil 0,1 de la variable age para los pacientes diabéticos de nuestra simulación? No podemos responder a esta pregunta con los propios datos históricos dia: como sólo tiene 442 registros, no podemos calcular cuál es el valor de la milésima parte. En su lugar, puedes aprovechar los resultados de una simulación de Montecarlo, ¡lo que harás ahora!

El conjunto de datos de la diabetes se ha cargado como un DataFrame, dia, y se han importado para ti las siguientes bibliotecas: pandas como pd, numpy como np, y scipy.stats como st.

¿Por qué hay que explorar los datos?

Simulaciones Montecarlo en Python

Ejercicio interactivo práctico