{"id":39949,"date":"2020-03-25T11:53:34","date_gmt":"2020-03-25T11:53:34","guid":{"rendered":"https:\/\/keepler.io\/?p=39949"},"modified":"2023-11-27T12:02:35","modified_gmt":"2023-11-27T12:02:35","slug":"paradigma-bayes-escenarios-incertidumbre","status":"publish","type":"post","link":"https:\/\/keepler.io\/es\/2020\/03\/25\/paradigma-bayes-escenarios-incertidumbre\/","title":{"rendered":"El paradigma Bayes para escenarios de incertidumbre"},"content":{"rendered":"<p>\u00bfAcaso nunca os ha pasado que, paseando por el campo y la monta\u00f1a, nos topamos con un camino que se separa en dos, y no sabemos hacia d\u00f3nde ir? En una situaci\u00f3n as\u00ed, lo \u00fanico que sabemos es que estamos en un camino que se divide en dos, y por tanto:<\/p>\n<ul>\n<li>Tenemos que volver por donde hemos venido.<\/li>\n<li>Continuamos por uno de los dos caminos que no sabemos a d\u00f3nde nos llevar\u00e1.<\/li>\n<\/ul>\n<p>Lo natural es que tomemos decisiones racionales, y que si, por ejemplo, estamos cansados, volvamos por donde hemos venido. O si nos vemos con fuerzas, y hace buen tiempo, sigamos por alguno de los dos caminos.\u00a0<\/p>\n<p>Y es que en nuestro d\u00eda a d\u00eda nos topamos con este tipo de decisiones donde, \u201csabiendo lo que sabemos\u201d, tenemos que elegir en un escenario de incertidumbre.<\/p>\n<h3>El paradigma Bayes<\/h3>\n<p>A fin de cuentas, la incertidumbre no es m\u00e1s que un estado en el que el conocimiento se ve limitado por la falta informaci\u00f3n para describir futuras acciones o eventos.<\/p>\n<p>La ciencia de datos y la modelizaci\u00f3n estad\u00edstica han ido desarrollando a lo largo de los a\u00f1os numerosas formas de poder cuantificar la incertidumbre con el fin de mejorar las predicciones de futuros eventos. Existen numerosas disciplinas que hacen uso de la cuantificaci\u00f3n de la incertidumbre, como la ingenier\u00eda, la f\u00edsica, la meteorolog\u00eda, y la econom\u00eda. Uno de los beneficios de utilizar un <strong>enfoque <\/strong><strong><i>bayesiano<\/i><\/strong> (todo empez\u00f3 con este <a href=\"https:\/\/en.wikipedia.org\/wiki\/Thomas_Bayes\" target=\"_blank\" rel=\"noopener noreferrer\">maestro de la estad\u00edstica<\/a>) es que precisamente estamos incorporando la incertidumbre a nuestra toma de decisiones.<\/p>\n<p>Hoy, venimos a hablar del m\u00e9todo m\u00e1s extendido para cuantificar esa incertidumbre. En la siguiente parte de esta serie de posts, realizaremos alg\u00fan ejemplo aplicado a Machine Learning. De momento, vamos a empezar describiendo un poco los diferentes paradigmas que existen para cuantificar la probabilidad y la incertidumbre de un conjunto de datos con un enfoque Bayesiano. \u00bfPara qu\u00e9 es necesario esto? Bueno, no siempre tenemos una \u00fanica soluci\u00f3n a un problema, si no que por lo general, y m\u00e1s en ciencia de datos, disponemos de muchas. Por tanto, entender las diferencias entre uno y otro nos puede servir de base para empezar con buen pie.\u00a0<\/p>\n<h3><b>Los tipos de enfoque<\/b><\/h3>\n<p>Los enfoques para modelar un conjunto de datos m\u00e1s extendido suelen venir de los \u201c<strong>f\u00edsicos<\/strong>\u201d y los \u201c<strong>evidenciales<\/strong>\u201d. Estos enfoques hacen uso de la probabilidad para realizar estudios sobre un conjunto de datos, observaciones, registros, eventos, etc. Cualquier sistema de inteligencia artificial basado en Machine Learning hace (casi siempre) uso de alguno de estos enfoques para realizar predicciones o inferencias, y es que:<\/p>\n<ul>\n<li>La <strong>inferencia frecuentista<\/strong> extrae conclusiones provenientes de los datos mediante el estudio de las frecuencias o proporciones de los mismos. Cuando estudiamos un fen\u00f3meno con un enfoque frecuentista, estamos asumiendo que nuestros datos con los que hacemos el estudio o experimento, son independientes del siguiente experimento (con los mismos pasos e hip\u00f3tesis) que potencialmente podr\u00edamos repetir. As\u00ed, el logro frecuentista reside en que, si repiti\u00e9ramos este experimento todas las veces que quisi\u00e9ramos, obtendr\u00edamos las mismas conclusiones, verificando as\u00ed que el fen\u00f3meno que estamos estudiando no es aleatorio. En otras palabras, que nuestra hip\u00f3tesis se cumple cada vez que realizamos el experimento.<\/li>\n<li>La <strong>inferencia bayesiana o probabil\u00edstica<\/strong>, por el contrario, consiste en ir actualizando una creencia como consecuencia de la observaci\u00f3n de nuevas evidencias. Se llega al descubrimiento mediante un razonamiento aproximado, donde no existen verdades absolutas, y la incertidumbre nos hace dudar del grado de consistencia de nuestras hip\u00f3tesis. As\u00ed, se pretende encontrar las probabilidades de nuestras hip\u00f3tesis condicionadas a los datos y a las evidencias que tenemos.<\/li>\n<\/ul>\n<p>Hist\u00f3ricamente, siempre ha habido disputas dentro de la comunidad cient\u00edfica por determinar cu\u00e1l es el mejor m\u00e9todo para estudiar el mundo que nos rodea. Al principio, el probabilismo (relacionado con lo que venimos denominando como lo bayesiano), no era m\u00e1s que un paradigma epistemol\u00f3gico. Con el pasar de los a\u00f1os, se convirti\u00f3 en otro arma m\u00e1s para la ciencia de cara a realizar experimentos y contrastar hip\u00f3tesis. <\/p>\n<h3><b>En la pr\u00e1ctica<\/b><\/h3>\n<p>La inferencia bayesiana trata de actualizar una probabilidad a posteriori como consecuencia de dos antecedentes: una probabilidad a priori y una funci\u00f3n de verosimilitud. Expresado matem\u00e1ticamente, tenemos:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30618\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-1.png?resize=275%2C85&#038;ssl=1\" alt=\"\" width=\"275\" height=\"85\" \/><\/p>\n<ul>\n<li><strong>A<\/strong> es el evento o hip\u00f3tesis que estamos estudiando y cuya probabilidad de suceso podr\u00eda verse afectada por los datos o evidencia que tenemos.\u00a0<\/li>\n<li><strong>P (A)<\/strong> es la probabilidad estimada de que el evento A suceda <\/span><b>antes<\/b> de que se observan nuevos registros de B.<\/li>\n<li><strong>B<\/strong> son los datos o evidencia de las que disponemos.<\/li>\n<li><strong>P (A \/ B)<\/strong> es la probabilidad posterior, es decir, aquella de que el evento suceda <b>despu\u00e9s<\/b> de haber visto B.<\/li>\n<li><strong>P (B \/ A)<\/strong> es la probabilidad de observar B dado A, o lo que es lo mismo, es la verosimilitud. Indica la compatibilidad de nuestros datos con el evento que estamos estudiando.<\/span><\/li>\n<li><strong>P (B)<\/strong> es la verosimilitud marginal.\u00a0<\/li>\n<\/ul>\n<p>Este teorema tuvo mucho impacto en la ciencia, y lo sigue teniendo a d\u00eda de hoy. Por ejemplo, podemos resolver algoritmos de inferencia estad\u00edstica aplicando esta l\u00f3gica para incorporar a nuestras predicciones el concepto de incertidumbre. Por ejemplo, la regresi\u00f3n lineal. Es un sencillo algoritmo que intenta explicar o predecir una variable usando una combinaci\u00f3n lineal de otras. El enfoque frecuentista para realizar inferencias es sencillo, y es que en una regresi\u00f3n lineal, simplemente tenemos que pasar los datos (matriz X) que queramos inferir a la siguiente ecuaci\u00f3n para obtener nuestras predicciones Y:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30621\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-2.png?resize=163%2C51&#038;ssl=1\" alt=\"\" width=\"163\" height=\"51\" \/>De entre todas las formas que existen para encontrar los mejores par\u00e1metros , la siguiente forma es la m\u00e1s extendida:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30624\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-3.png?resize=184%2C67&#038;ssl=1\" alt=\"\" width=\"184\" height=\"67\" \/>Sin embargo, existe tambi\u00e9n la regresi\u00f3n lineal con enfoque bayesiano. Aqu\u00ed, asumimos que nuestra variable a predecir sigue una distribuci\u00f3n de probabilidad. La idea de una regresi\u00f3n lineal bayesiana, no es conseguir calcular los mejores coeficientes o par\u00e1metros que minimicen el error de entrenamiento y test, si no determinar cu\u00e1l es la distribuci\u00f3n posterior de los par\u00e1metros o coeficientes:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30628\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-4.png?resize=315%2C78&#038;ssl=1\" alt=\"\" width=\"315\" height=\"78\" \/><\/p>\n<p>Nuestra variable de respuesta tendr\u00e1 la siguiente forma:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30631\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-5.png?resize=227%2C68&#038;ssl=1\" alt=\"\" width=\"227\" height=\"68\" \/>En la pr\u00e1ctica, tendremos que definir primero cu\u00e1les son las distribuciones a priori de nuestros par\u00e1metros, nuestras creencias, lo que conocemos, para ir actualizando nuestro distribuci\u00f3n a posteriori v\u00eda un algoritmo de sampleo como <a href=\"https:\/\/towardsdatascience.com\/markov-chain-monte-carlo-in-python-44f7e609be98\" target=\"_blank\" rel=\"noopener noreferrer\">MCMC<\/a> (que, a lo mejor, podr\u00edamos explicar en otro post). As\u00ed, una de las principales ventajas de usar un enfoque bayesiano es que:<\/p>\n<ul>\n<li>Podremos incorporar esa informaci\u00f3n a priori en el c\u00f3mputo y soluci\u00f3n. Si no tenemos esa informaci\u00f3n valiosa, podemos siempre recurrir a la distribuci\u00f3n normal de probabilidad.<\/span><\/li>\n<li>Si no disponemos de muchos datos, estaremos cuantificando la incertidumbre de la aparici\u00f3n de nuevos datos a todo el proceso, de manera que cuando nuestro set de datos tienda a ser infinito, los par\u00e1metros convergen con los calculados por el enfoque frecuentista.<\/li>\n<\/ul>\n<p>Al final, nuestro resultado ser\u00e1 que, en vez de obtener un n\u00famero que nos defina al 100%, nuestras predicciones, obtendremos un rango, o mejor dicho, una distribuci\u00f3n de probabilidad de nuestros coeficientes. <\/p>\n<h3><b>Un caso de uso real en un escenario de incertidumbre<\/b><\/h3>\n<p>En el d\u00eda en el que se escribi\u00f3 este art\u00edculo, se dieron numerosos casos positivos de COVID-19 en el mundo. En <a href=\"https:\/\/github.com\/twiecki\/covid19\" target=\"_blank\" rel=\"noopener noreferrer\">este magn\u00edfico post<\/a>, el autor nos ilustraba c\u00f3mo en el inicio de la pandemia, la curva de infecciones en la primera etapa de la epidemia era de tipo exponencial. <strong>Usando una regresi\u00f3n bayesiana, en este caso no lineal sino exponencial<\/strong>, pod\u00edamos hacer una tentativa de predicci\u00f3n de c\u00f3mo ser\u00edan los primeros 30 d\u00edas en las comunidades aut\u00f3nomas de nuestro pa\u00eds, con la forma:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-30634\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Paradigma-bayes-formula-6.png?resize=364%2C127&#038;ssl=1\" alt=\"\" width=\"364\" height=\"127\" \/>Y es el n\u00famero de infectados, M es la constante, B es la tasa de crecimiento del n\u00famero de infectados, y es un par\u00e1metro para a\u00f1adir aleatoriedad, sigma, a todo el proceso. En Keepler realizamos el mismo ejercicio propuesto por el autor, y obtuvimos este resultado:<\/p>\n<p><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-full wp-image-30637\" src=\"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/10\/Regresion-bayesiana-prediccion.jpg?resize=1080%2C1589&#038;ssl=1\" alt=\"\" width=\"1080\" height=\"1589\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>La l\u00ednea roja marca los casos reales que ha habido desde el caso n\u00famero 100. El eje de las X representa los d\u00edas que hab\u00edan pasado desde el caso n\u00famero 100. En el eje de las Y, el n\u00famero de casos confirmados posibles en un escenario de incertidumbre, simulado 500 veces. A pesar de ser un modelo que asum\u00eda una falta de pol\u00edticas de control de la pandemia, como la cuarentena o el distanciamiento social, ni tampoco se ten\u00eda en cuenta la exposici\u00f3n al virus, y\/o la densidad de poblaci\u00f3n, nuestro modelo fue capaz de, con muy pocos datos, darnos una visi\u00f3n general de lo que podr\u00eda ser el peor de los escenarios\u2026 Est\u00e1bamos, por tanto, cuantificando la incertidumbre en un escenario muy muy pesimista, \u00a1justo lo que busc\u00e1bamos!<\/p>\n<p>Agradecemos a <a href=\"https:\/\/twitter.com\/twiecki\" target=\"_blank\" rel=\"noopener noreferrer\">Thomas Wiecki<\/a><span style=\"font-weight: 400;\"> su contribuci\u00f3n a la comunidad con este magn\u00edfico ejemplo de modelado bayesiano.\u00a0<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u00bfAcaso nunca os ha pasado que, paseando por el campo y la monta\u00f1a, nos topamos con un camino que se separa en dos, y no sabemos hacia d\u00f3nde ir? En una situaci\u00f3n as\u00ed, lo \u00fanico que sabemos es que estamos en un camino que se divide en dos, y por tanto: Tenemos que volver por [&hellip;]<\/p>\n","protected":false},"author":134360170,"featured_media":33912,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_seopress_titles_title":"El paradigma Bayes para escenarios de incertidumbre","_seopress_titles_desc":"El paradigma Bayes, un enfoque para poder cuantificar la incertidumbre con el fin de mejorar las predicciones de futuros eventos.","_seopress_robots_index":"","_seopress_robots_follow":"","_seopress_robots_imageindex":"","_seopress_robots_snippet":"","_seopress_robots_primary_cat":"none","_seopress_robots_breadcrumbs":"","_seopress_robots_freeze_modified_date":"","_seopress_robots_custom_modified_date":"","_seopress_robots_canonical":"","_seopress_social_fb_title":"","_seopress_social_fb_desc":"","_seopress_social_fb_img":"","_seopress_social_fb_img_attachment_id":0,"_seopress_social_fb_img_width":0,"_seopress_social_fb_img_height":0,"_seopress_social_twitter_title":"","_seopress_social_twitter_desc":"","_seopress_social_twitter_img":"","_seopress_social_twitter_img_attachment_id":0,"_seopress_social_twitter_img_width":0,"_seopress_social_twitter_img_height":0,"_seopress_redirections_value":"","_seopress_redirections_enabled":"","_seopress_redirections_enabled_regex":"","_seopress_redirections_logged_status":"both","_seopress_redirections_param":"","_seopress_redirections_type":301,"_seopress_analysis_target_kw":"","_et_pb_use_builder":"","_et_pb_old_content":"","_et_gb_content_width":"","content-type":"","_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_wpcom_ai_launchpad_first_post":false,"_jetpack_feature_clip_id":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"{title}\n\n{excerpt}\n\n{url}","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2},"_wpas_customize_per_network":false,"jetpack_post_was_ever_published":false},"categories":[227],"tags":[325,299,283],"class_list":["post-39949","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-ciencia-de-datos","tag-ciencia-de-datos-2","tag-datos-es","tag-tecnologia"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/keepler.io\/wp-content\/uploads\/2020\/04\/Keepler-paradigma-bayes-escenarios-incertidumbre-1.jpg?fit=850%2C350&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p9CeZw-aol","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/posts\/39949","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/users\/134360170"}],"replies":[{"embeddable":true,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/comments?post=39949"}],"version-history":[{"count":7,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/posts\/39949\/revisions"}],"predecessor-version":[{"id":39959,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/posts\/39949\/revisions\/39959"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/media\/33912"}],"wp:attachment":[{"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/media?parent=39949"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/categories?post=39949"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/keepler.io\/es\/wp-json\/wp\/v2\/tags?post=39949"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}