Rarezas: estadística algebraica

Matemáticas y estadística son peras y manzanas. La una es la ciencia del ser; la otra, del parecer. Se encuentran en la teoría de la probabilidad, pero se miran de reojo y con recelo. Por eso este curso de estadística algebraica es toda una rareza. Contiene resultados, como la proposición 1.1.2 que… bueno, sí, bien, vale: Proposición 1.1.2. Las variables aleatorias discretas X e Y son independientes sí y solo sí la matriz $p = (p{ij})$ tiene rango 1._ ...

14 de enero de 2015 · Carlos J. Gil Bellosta

evtree: árboles globales

Tengo por delante otro proyecto que tiene mucho de análisis exploratorio de datos. Sospecho que más de un árbol construiré. Los árboles son como la Wikipedia: prácticamente nunca el último pero casi siempre el primer recurso. Esta vez, además, por entretenerme un poco, probaré el paquete [evtree](http://cran.r-project.org/web/packages/evtree/index.html). Aunque no porque espere sorprendentes mejoras con respecto a los tradicionales, ctree y rpart. ¿Qué tiene aquel que los diferencie de los otros dos? Que la optimización es global. Tanto ctree como rpart utilizan algoritmos recursivos: al definir un nuevo corte del espacio, el algoritmo solo tiene en cuenta la región definida por los cortes anteriores. La optimización es local. evtree utiliza un algoritmo global de la familia de los evolucionarios (¡qué tufillo a lentorro!). Los detalles están aquí. ...

12 de enero de 2015 · Carlos J. Gil Bellosta

La curtosis de una variable aleatoria constante

Una mañana de hace veinte $\pm \epsilon$ años sufrí mi primera hora de clase de estadística reglada. No la olvidaré: fue un monográfico sobre momentos muestrales de todo orden; los sumatorios se salían por ambos márgenes de las transparencias de acetato. Horrible. Sin embargo, aquel día perdí la ocasión de levantar la mano y preguntar por la curtosis de una variable aleatoria constante. Porque necesito un valor razonable por defecto y no se me ocurre ninguno. ¿Cero acaso? ¿Alguna sugerencia? ...

5 de enero de 2015 · Carlos J. Gil Bellosta

El problema de la estimación inversa

Supongamos que tenemos unos niños de los que sabemos las edades $x_i$ y las alturas $y_i$. Supongamos además que podemos estimar las segundas en función de las primeras con un modelo lineal clásico $$ y_i \sim N(a_0 + a_1 x_1, \sigma).$$ Este modelo nos permite, dada una edad, estimar la altura y los correspondientes intervalos de confianza. Pero, dada una altura, ¿qué nos dice de la edad? Este es el problema conocido como de la estimación inversa. ...

31 de diciembre de 2014 · Carlos J. Gil Bellosta

Cómo no nació el "big data"

En julio anuncié en mi cuenta de Twitter (léase de abajo a arriba): Ya está disponible.

30 de diciembre de 2014 · Carlos J. Gil Bellosta

Modelos mixtos por doquier

Los códigos postales, por ejemplo, son un problema a la hora de crear modelos predictivos: son variables categóricas con demasiados niveles. Así, por ejemplo, los bosques aleatorios de R solo admiten variables categóricas con no más de 32 niveles. Hay trucos de todo tipo para mitigar el problema. Hace un año, Jorge Ayuso me puso sobre la pista de uno de los que tiene más recorrido. Consiste en [su versión más simplificada en]: ...

29 de diciembre de 2014 · Carlos J. Gil Bellosta

¿Hubo alguna vez un millón de palentinas?

En el año 2013 hubo 54 muertes de mujeres por violencia de género. Eso da una tasa nacional de poco más de dos por millón (de mujeres). El Mundo nos lo ha querido mostrar su distribución provincial así: Diríase que la tasa palentina es enorme, cinco veces la nacional. Pero en Palencia viven del orden de cien mil mujeres y hubo un único caso en 2013 (además, ni la mujer ni el agresor, se ve, eran de la provincia sino de un pueblo limítrofe de Cantabria; solo que el cadáver apareció en el sur de la linde). ...

10 de diciembre de 2014 · Carlos J. Gil Bellosta

Ruido en las estadísticas oficiales

Hacía tiempo que no hablaba de este tema. Pero han salido de mi LIFO de artículos potencialmente interesantes dos a los que merece la pena echar un ojo. El primero, este, arranca con Los organismos estadísticos gubernamentales suelen publicar los estimadores puntuales de las estadísticas económicas oficiales. La documentación metodológica de dichos organismos puede hacer constar que estas estimaciones están sujetas a incertidumbre, pero no suelen cuantificarla. Las notas de prensa raramente discuten el error potencial. ...

9 de diciembre de 2014 · Carlos J. Gil Bellosta

La correlación ni siquiera implica "correlación"

Esto es, según Andrew Gelman, la correlación entre dos variables en una muestra ni siquiera implica su “correlación” (entre comillas, por distinguirlas) en la población de interés. El enlace anterior también discute otras variantes del archiconocido “la correlación no implica causalidad”, tales como la causalidad está correlacionada con la correlación, la falta de correlación está correlacionada con la falta de causalidad, etc. que, si yo fuera tú, me apresuraría a consultar en el enlace anterior. ...

8 de diciembre de 2014 · Carlos J. Gil Bellosta

Sí, señor ministro... y encuestas

Cayó en mis manos que son los resultados de una encuesta en la que la misma pregunta (en puridad, una pregunta sobre una cuestión global y otra sobre un asunto particular de la anterior) reciben respuestas manifiestamente contrarias y contradictorias por parte de una muestra del ostentador de la soberanía. Lo cual me recordó que hacía tiempo había dado con https://www.youtube.com/watch?v=G0ZZJXw4MTA extraído de Yes, Minister y que en inglés sin subtitular ilustra muy amenamente los efectos que sobre el público tiene la manera en que se plantean las cuestiones.

4 de diciembre de 2014 · Carlos J. Gil Bellosta