Las prioris no informativas están manifiestamente sobrevaloradas

La estadística bayesiana se enseña en cursos de estadística (y, frecuentemente, envuelto en un aparataje matemático tan ofuscante como innecesario). Lo malo es que en los cursos y textos de estadística no existe información previa. La información previa sobre los fenómenos en los que se utilizaría la estadística bayesiana están en las aplicaciones, extramuros del muy agnóstico mundo de la estadística y la matemática. Por eso, a los autores de los libros de estadística bayesiana y quienes enseñan cursos sobre lo mismo, enfrentados al problema de llenar de sentido la problemática distribución a priori, no se les ocurre nada mejor que discutir muy sesudamente la excepción (la priori no informativa) en lugar de la regla (la priori informativa). Reto al lector escéptico a que repase cualquier manual en la materia (que no haya sido escrito por Gelman) y compare el espacio que dedican a la selección de prioris no informativas con el de convenir una priori informativa decente. ...

4 de enero de 2016 · Carlos J. Gil Bellosta

¿Por qué el empate de la CUP es más raro de lo que parece (y de lo que yo mismo digo)?

Menos el de un presunto profesor, La probabilidad de que 3030 votantes en la #ANECUP empaten con 1515 votos es 1 / 3029 = 0.00033014, lo que corresponde a un suceso imposible – Mario Bilbao (@mario_bilbao) December 27, 2015 todos los análisis que he visto al respecto (1, 2, 3), incluido el mío, coinciden en señalar que la probabilidad de empate en el muy manido acto asambleario de la CUP es relativamente alta: alrededor del 1,5%. Y más todavía si se tienen en cuenta los resultados de las votaciones previas. ...

30 de diciembre de 2015 · Carlos J. Gil Bellosta

Empates electorales, sorteos y una inadvertida paradoja

Ayer hubo en España una asamblea en la que 3030 personas votaron sobre un cierto asunto: 1515 votaron que sí y otras tantas, que no. La cosa acabó en empate. Una estimación (de las muchas que pueden hacerse con hipótesis de partida distintas) de la probabilidad de que tal ocurra es dbinom(1515, 2*1515, 0.5), que viene a ser 1.5%. Expresado de otra manera, solo ocurriría en uno de cada 70 congresos de tales características. ...

29 de diciembre de 2015 · Carlos J. Gil Bellosta

La intersección de lo interesante, lo sorprendente, lo cierto y lo basado en datos

Me interesan, obviamente, los problemas interesantes. Me interesan los problemas en que puedo argumentar basándome en datos. Para ello, obviamente, de nuevo, tienen que existir datos con los que tratar de dar respuesta a esas preguntas interesantes del párrafo anterior. Me interesa que los datos revelen respuestas no obvias, que no sepamos ya de antemano. Me interesa que los datos me sorprendan. Me interesa, obviamente, que esas respuestas sorprendentes a preguntas interesantes basadas en datos sean ciertas. O que tengan visos de certeza. Cuando menos, que no se caigan en la primera revisión crítica. ...

29 de diciembre de 2015 · Carlos J. Gil Bellosta

Un poco más sobre el índice de poder de Banzhaf

En el año 2012 escribí esto, que incluye El índice de Banzhaf para un determinado partido político mide su poder en términos del porcentaje de las posibles alianzas mínimas ganadoras en las que participa dentro de su universo total. Una alizanza es ganadora cuando reúne más de la mitad de los votos. Y es mínima cuando todos sus integrantes son necesarios para que sea ganadora; excluye, por ejemplo, la alianza trivial formada por todos los partidos. ...

23 de diciembre de 2015 · Carlos J. Gil Bellosta

Coaliciones de Banzhaf en el 20D

Usando código de una entrada anterior voy a medir el poder de cada partido político de acuerdo con Banzhaf tras las elecciones de diciembre de 2015. escannos <- c(123, 90, 69, 40, 9, 8, 6, 2, 2, 1) names(escannos) <- c( "pp", "psoe", "pod", "c's", "erc", "dl", "pnv", "iu", "bildu", "cc") banzhaf(escannos) da 14 coaliciones mínimas, pp psoe pp pod pp c's erc dl pp c's erc pnv pp c's erc iu bildu pp c's dl pnv pp c's dl iu bildu cc psoe pod c's psoe pod erc dl psoe pod erc pnv iu psoe pod erc pnv bildu psoe pod dl pnv iu bildu psoe pod dl pnv iu cc psoe pod dl pnv bildu cc y un reparto de poder que queda de esta manera: psoe pod pp pnv dl erc c's iu bildu cc 57.1 57.1 50.0 50.0 50.0 42.8 42.8 35.7 35.7 21.4 Gráficamente, ...

22 de diciembre de 2015 · Carlos J. Gil Bellosta

La primera víctima será la sintaxis

No seré yo quien condene el contar palabras en textos para inferir por dónde van los tiros. Lo he hecho y lo hago con frecuencia. Pero lo cuestiono en algunas ocasiones. La principal, esos análisis sumarísimos de los programas electorales —p.e., en la campaña que acabamos de cerrar— que los comparan en términos del número de veces que se mencionan diversos términos o el número de páginas que dedican a diversos temas. ...

21 de diciembre de 2015 · Carlos J. Gil Bellosta

La respuesta es: "Yo no"

La pregunta, a propósito, es esta.

18 de diciembre de 2015 · Carlos J. Gil Bellosta

CRAN, r-devel, GitHub, Travis CI, pruebas automáticas y todo eso

Estoy harto. La gente de CRAN me devolvió (con errores) un paquete que trataba de subir. Había hecho el prescriptivo R CMD check --as-cran etc. y el log era una patena. Pero había un par de NOTES al pasar el test sobre la versión de desarrollo de R, r-devel. No solo hay que probar los paquetes en la versión que hay sino también en la que vendrá (tal y como está docuentado). ...

17 de diciembre de 2015 · Carlos J. Gil Bellosta

La búsqueda de la causa más probable de un efecto

La búsqueda de la causa más probable de un efecto tiene un nombre: razonamiento abductivo. Que el visitante al enlace anterior aprenderá distinto del deductivo y el inductivo. Y que los viejos de estas páginas reconocerán en esta entrada que la formaliza y cuantifica en un caso concreto.

16 de diciembre de 2015 · Carlos J. Gil Bellosta