Mi versión modificada puede bajarse de aquí y permite folding y publicar el RSS del status para sindicar. Esta última opción hace pública tu key. Siendo que es obvio que tu status es público y que no podes usarlo para cosas sensibles.. ¿a alguien se le ocurre las consecuencias que esto conlleva?

Para borrar mis huellas
destrocé mi camisa,
confundí con estrellas
las luces de neón.

Hice trampas al poker.
Defraudé a mis amigos.
Sobre el banco de un parque,
dormí como un lirón

Por decir lo que pienso,
sin pensar lo que digo,
más de un beso me dieron
y más de un bofetón.

Lo que sé del olvido
lo aprendí de la luna.
Lo que sé del pecado
lo tuve que buscar

como un ladrón debajo
de las faldas de algunas
de cuyo nombre ahora
no me quiero acordar.

Joaquín

Información útil:

• Qué: Seminario “Touch of Class: How we teach introductory programming”
• Por quién: Bertrand Meyer
• Dónde: Aula Magna de la FRBA-UTN (Medrano 951)
• Cuándo: miércoles 11 de febrero, a las 19hs
• Para quién: Para todos, es abierto y gratuito. No se necesitan conocimientos previos en ninguna materia en particular. Solo hay que ir.

Es otro de los excelentes seminarios organizados por el grupo Athena (aunque su nueva página web sea un insulto a la accesibilidad)

El hecho es que tuve que ingeniármelas en conseguir un tren cuanto antes, para no perder mi vuelo. Por suerte, no tuve mayores complicaciones y pude estar en Frankfurt a las 2pm, con suficiente tiempo para caminar 4 horas por la ciudad. La visita fue algo fugaz, así que no puedo contar mucho.

De la ciudad que sí puedo contar algo es de Amsterdam. Y acá quiero dar por tierra el comentario clásico que uno recibe cuando visita una ciudad donde el vicio es tan legal como guiñar el ojo. Lo digo con mi mayor seriedad posible: el centro de la ciudad cansa al poco tiempo. Las luces rojas, los sex-shops y los coffee-shops son pintorescos al principio, pero después de un rato se vuelven repetitivos y hasta chabacanos. Supongo que lo entretenido que puede tener se pierde con tanto factor for-export y expositivo. Y dado que para consumir (lo que sea) no hay voluntad y/o plata me entretuve más con la linda arquitectura en las casas cercanas a Vondelpark, el museo Van Gogh y Tour a la gorra por la ciudad.

Las fotos (ajenas, porque recuerden que mi cámara murió) estarán en línea esta tarde. O no.. visste’ como es esto.

UPDATE 19Jan: Las fotos, finalmente, de Berlín y 25C3.

Cumplido ello, al foco. Tuve la oportunidad de participar del 25th Chaos Communication Congress (25C3), lo que fue una experiencia realmente increíble. La gente, los anuncios, las mesas de distintos proyectos y todas esas cosas hicieron de CCC uno de mis eventos favoritos.

El año terminó y con él el tema del bug de OpenSSL/Debian. Solo me queda entregar mi informe en el posgrado y defenderlo, pero espero terminarlo cuanto antes y enterrar así el tema. Es hora de cosas nuevas.

En cuanto al resto, pasé un fin de año muy atípico. Con amigos y desconocidos. Raro pero bien. También en estos días me encuentro aplicando para un interesante puesto laboral (como dudo que se concrete, no daré más detalles), con lo que es algo más en mi cabeza. Durante una caída, mi cámara de fotos dejó de funcionar como corresponde, así que pueden ver las fotos de terceros en estos lugares: las de Machu y las de Maxi. Estoy tratando de hacer una subselección para poner en el lugar de siempre, pero no me he hecho del tiempo. En cuanto llegue a casa veré cómo me hago de una nueva cámara.

Como dije, ahora estoy visitando Amsterdam, idea que surgió medio improvisada. El 6 tengo pasaje a casa desde Frankfurt, por lo que tampoco puedo colgarme mucho.

Ahora a dormir, que me estoy cayendo del sueño…

De Buenos Aires a Sao Pablo, en un vuelo repleto de adolescentes que viajaban en grupo. Griteríos y canciones. Espera de 2 horas en el aeropuerto de Sao Pablo. De Sao Pablo a Frankfurt. En Frankfurt 8 horas esperando un tren lechero a Halle. En un cafe de Halle, 4 horas esperando salir a Berlín. En Berlín, 3 horas perdido en el metro para llegar a casa de Martin.

Pero ya estoy acá. Y tengo una cintita en mi muñeca que dice que los próximos días tendré acceso a la 25C3. I’m happy :D.

Las primeras fotos (no mías).

I’ll be out of home from 24th Dec until 7th Jan. It’ll be a non-traditional holidays, be far away from family. But, in the other hand, it gonna be awesome be and talk in a CCC and spend a Christmas Eve in a plane and a New Year’s Eve in Brandenburg Gate with friends.

See you there?

Para los que no hicieron click en el link anterior y dado que la entrada en la Wikipedia en español sobre el tema deja bastante que desear (estoy corrigiéndola), acá está mi breve explicación:

La ley de Benford, también conocida como la ley del primer dígito, dice que, en los números que existen en la vida real^TM, la primer cifra tiene muchas más posibilidades de ser 1 que otro valor. Además, según crece este primer dígito, más improbable es que se encuentre en la primera posición.

A por un ejemplo, que seguro es más fácil. Tomemos una tabla con estadísticas cualquiera, como ser el área o población de las provincias argentinas, la capacidad de los estadios del mundo, los muertos por accidentes de tránsito o las estadísticas de visitas en tu blog. Como todo en estadística, mientras más grande la muestra mejor, así que vayan a por tablas realmente grandes.

Intuitivamente, uno podría pensar que, para cualquier número, las posibilidades que empiece con 1 son las mismas que con 9. Es decir, al agarrar una pelota de números cualquiera de la vida real^TM se podría llegar a creer que, si la cantidad es lo suficientemente grande, más o menos 1/9 de la muestra empezarán con 1 (nótese que no tiene sentido que los valores empiezen en 0, por lo que las opciones son 9). Esto es porque creemos que los números que estamos analizando se comportan como si fuesen aleatorios. Si tiramos una moneda al aire una gran cantidad de veces, cerca de la mitad de las oportunidades será seca. Es algo que aprendimos hace mucho y nos parece intutivo que la naturaleza se comporte así. Como si $deity hubiese tirado un dado gigante para decidir el largo de un río, la población de un país o el precio de las acciones en el MERVAL.

Como en mi disco todavía tengo los datos utlizados para el post la chica bajo la curva (dating pool) voy a utilizar la cantidad de casados, por provincia, por edad (puede descargarse desde acá). En mi caso son 30843 regitros. Esta linea de bash cuenta cuantos de los valores empienzan con cada cifra:

$ for i in $(cat provincia_indec.csv | grep años | cut -d ';' -f 4); do echo ${i:0:1}; done | sort | uniq -c
10350 1
5581 2
3550 3
2744 4
2159 5
1903 6
1815 7
1406 8
1335 9
Puede verse que el 33.56% de las cifras empiezan con 1 y que, mientras mayor es el valor del primer dígito, menor es la cantidad de ocurrencias.

En efecto, nuestro nuevo amigo Benford describe este fenómeno y nos dice que, la probabilidad p de que el dígito d aparezca en el primer lugar está dado por la siguiente fórmula:

En el ejemplo anterior, el dígito 1 se encuentra en 10350/30843=0.336 de los casos. Puede verse que la predicción de la fórmula es bastante buena, ya que log₁₀(1+1)=0.301. En la siguiente figura puede verse como se ajusta la fórmula a los datos de la práctica (primer dígito de la cantidad de personas casadas por provincia, por edad según el censo 2001):

De hecho, la misma fórmula puede aplicarse a más de un dígito. Por ejemplo, la probabilidad de que una cifra empiece con 42 (primer dígito 4, segundo 2) es log₁₀(1+1/42)=0.010219. Modificando levemente el script (${i:0:2}) podemos estudiar la cantidad de cifras por la repetición de sus primeros dos números y compararlos con su valor teórico:

Impresionante.. no?

¿Y porqué pasa esto? Ocurre que a las magnitudes del mundo real^TM están distribuidas de forma logarítmica.

Recordemos la fórmula: p_(d)=log₁₀(1+1/d)=log₁₀(d + 1) − log₁₀(d)
Es decir, cuenta cuántos números hay entre d y d+1 dentro de la escala logarítmica.

La mejor explicación que recibí para este fenómeno habla de un cambio de escala. Supongamos por un momento que la distribución de los primeros dígitos de lo largo de los ríos, lo alto de las montañas, lo profundo de los posos es constante. Ahora imaginemos que $deity se levanta una mañana y duplica el tamaño del planeta (o del universo, ustedes elijen). Las medidas que empezaban por 1 hora pasan a empezar por 2 o 3 (160*2=320). Lo que empezaba por 2 ahora lo hace por 4 o por 5 (290*2=580). El 3 se va a 6 o 7 (384*2=768). El 4 a 8 y 9.

¡Pero todos aquellas medidas que empezaban por 5, 6, 7, 8 y 9 ahora empiezan por 1! Si realizamos esta operación varias veces los valores se amontonan rápidamente en los iniciados por 1, generando la escala logarítmica en cuestión. Y acá está el tema. La mayoría de los valores de la vida real^TM son resultados de multiplicaciones.

El siguiente gráfico invita a comparar las superficies rojas (valores que inician con 1) y azules (valores que inician con 8) para estudiar sus probabilidades de ocurrencia en el primer dígito:

Las distribuciones que cubren muchos órdenes de magnitud (que varían mucho entre número y número) cumplen relativamente bien con la ley de Benford. Sin embargo puede no ocurrir así siempre:

Nótese que la clave está en las grandes magnitudes (recuerden la explicación del universo que se duplica). Existen tablas de números de la vida real^TM que no cumplen la ley dada que estan acotadas en cuento a su rango, por ejemplo los datos del cierre del MERVAL de los últimos 3 años. Si bien son números grandes, su máximo y mínimo es acotado. Imagino (no tengo uno a mano) que los precios unitarios de los productos en un ticket de supermercado tampoco se agustan a la ley por razones parecidas. La tabla de goleadores de un torneo de fútbol padece el mismo trauma. ¿Se les ocurre algún otro ejemplo excepcional a la regla?

Espero hayan aprendido algo nuevo y ahora quieran a la matemática un poquito más :)

PD: Me olvidaba. Tarea para el hogar: Demostrar que la vida real^TM incluye a Fibonacci ;)

Hace un par de meses divagué con π. Este semana terminé de leer el libro La Proporción Áurea de Mario Livio. Este libro, que evidentemente trata sobre la constante φ (phi), lo compré en Madrid, atraído por la promesa de desasnarme con respecto a la historia de este número.

Primero me gustaría hablarles sobre el libro como tal, el cual resultó sumamente entretenido, aunque tal vez le sobren unas 50 o 100 páginas. Es un paseo por la historia del arte plástico, arquitectónico, escultural y matemático con φ como eje central hilador. Además de dar ejemplos de como se expresan los números en la naturaleza, explica con profundo detalle ciertas propiedades matemáticas en sus numerosos apéndices e invita al lector a pensar sobre porque la matemática encaja también como explicación del universo. Por momentos se va por las ramas pero casi siempre de forma afable y atrapante. Son altamente rescatables los pasajes donde el autor trata con ácida ironía los esfuerzos de los numerólogos de meter a φ lugares arbitrarios cualquiera, como en las dimensiones de Panteón. En mi opinión es un libro altamente recomendable para el aficionado a la matemática y mi calificación (hace mucho que no califico libros con estrellas) es .

En segundo lugar quisiera comentar algunas nerdeadas matemáticas inspiradas en el la lectura de este libro. Empezando por lo principio, presentar a quienes no conozcan el número φ, en las palabras de Euclides:

La longitud total a+b es al segmento más largo a como a es al segmento más corto b

El valor en cuestión se puede calcular despejando algebraicamente tomando b=1, y es:

Puede verse que, al igual que π(pi), es un número irracional. Pero, a diferencia de π, es un número algebraico, lo que lo emparenta bastante con teoría de números. Sin embargo, su expresión como fracción continua compuesta de solo 1s (unos), hace que converja muy lentamente, convirtiéndolo en el más irracional entre todos los irracionales. En términos llanos, φ puede expresarse como la siguiente fracción continua:

Nótese que, cada iteración de esta fracción continua puede expresarse como (los primeros 15 valores, esta secuencia fue resultado de este script):

2 / 1 = 2
3 / 2 = 1.5
5 / 3 = 1.666666666666666666666666667
8 / 5 = 1.600000000000000000000000000
13 / 8 = 1.625
21 / 13 = 1.615384615384615384615384615
34 / 21 = 1.619047619047619047619047619
55 / 34 = 1.617647058823529411764705882
89 / 55 = 1.618181818181818181818181818
144 / 89 = 1.617977528089887640449438202
233 / 144 = 1.618055555555555555555555556
377 / 233 = 1.618025751072961373390557940
610 / 377 = 1.618037135278514588859416446
987 / 610 = 1.618032786885245901639344262

Aquel lector atento notará que las fracciones tiene cierto patrón. Por un lado, el numerador de cada reglón pasa a ser el denominador en el siguiente. Los números son:
1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,….
¡Exacto! ¡Es la sucesión de Fibonacci! Cada número es la suma de los dos anteriores. En efecto, siendo F(n) el termino n-ésimo de la serie de Fibonacci:

Así es como puede verse una fuerte relación entre φ y Fibonacci. Y si con Fibonacci se puede calcular la proporción áurea… ¿cómo es la relación a la inversa?

Con esta formula se puede calcular el valor del término n-ésimo de la serie de Fibonacci:

Es fácil reconocer que dentro de esa fórmula está φ.

La propoción áurea φ tiene un montón de raras y divertidas propiedades. Entre las más atractivas se encuentran:

En lo personal aprendí muchas cosas nuevas sobre la proporción áurea. Pero en el apéndice 9 (como dije, es un libro con muchos apéndices) se comenta un concepto que me voló la cabeza: la ley de Benford. El tema se toca tangencialmente, promediando el último capítulo, como un ejemplo de como las matemáticas sorprenden. Y sí que lo hacen.

Pero ese será tema del siguiente post.

Podría extenderme hablando de las nuevas experiencias que viví y de lo bien que la pasé. Incluso intentar mencionar a las fascinantes personas que conocí durante ese mes, las nuevas oportunidades que se abrieron y lo lindo que es estar en otro país. Podría ahondar en comparaciones estériles y generalizaciones innecesarias como el contraste entre el homofóbico Marruecos y la Viena gay friendly. Podría aburrirlos con la pila de ideas ridículas y concepto aprendidos que ahora giran en mi cabeza.

Pero no. Se que reducir todo a cifras abstractas sería la opción más ridícula. Pero la carne es débil. Y pensar en un resumen o balance sin números y estadísticas me da escalofríos. Así que acá están:

Visité 8 ciudades en 4 países de 3 continentes. En 6 de esas ciudades hubo algún evento, en los que dí 8 conferencia/charla/clases, sumando 12 horas de hablar y hablar.

Viajé en 14 aviones de 5 compañías aereas que me obligaron a pisar 10 aeropuertos distintos. Entré o salí de Barajas 8 veces.

Mi mochila pesó un promedio de 20 kg cuando la facturé, lo que ocurrió 5 veces. Volví con 15 libros, de los cuales había llevado 7. Sumé 18 horas de usar mi portátil en el avión o aeropuerto, siempre off-line, tiempo durante el cual produje unas 305 lineas de Python entre otras cosas. También durante estos ratos empecé a leer 4 libros, de lo cuales terminé solo 1.

Dormí (considero dormir a cerrar los ojos durante 3 horas seguidas o más, independientemente de la horizontalidad o de si es de noche) en 12 camas/sofás/asientos distintos, de los cuales 4 fueron camas de hoteles, 4 camas de hostel y 2 camas/sofás de amigos varios. Me registré en hoteles/hostels 7 veces.

Saqué más de 1200 fotos, 27 con una cámara analógica. BTW, acabo de subir las de Alicante y Viena.

Gracias a todos los que organizaron los eventos, me soportaron y pagaron la cerveza. Un gusto haberlos conocido así como a otro montón de gente sumamente agradable. Fueron muy amables e hicieron que este mes de tour haya sido de las cosas más fuera de lo común e inolvidables que he vivido.

Había estado en Viena en el 2005, pero mi recuerdo estaba algo desdibujado. Igual algo seguía tal como lo estaba en la memoria, y es que tanto respeto por las normas y tanto ciudadanismo me hace sentir que estoy en ese irreal pueblo de Gilmore Girls. Todo es asquerosamente perfecto.

Se sabe que es una de esas ciudades que respira arte plástico y música clásica. También que es la casa de Gustav Klimt, de quien siempre me gustó su obra. Para aquellos a quien el nombre no les dice nada, seguramente reconocerán el particular estilo de su arte. La pintura Der Kuβ (El beso) me parece de un romanticismo único. Así que fue una buena oportunidad para reencontrarme con la particular arquitectura de la ciudad y con Klimt.

El evento, DeepSec, no estuvo mal. Aunque esperaba charlas de mejor nivel y la vida social de entre-pasillo es bastante apagada, estilo germano. De cualquier forma hubo momentos más que rescatables. En el campamento de lockpicking hice mis primeros pasos en este lindo/adictivo arte y pasatiempo, gracias a los muchachos de Toool. También conocí el emblemático Metalab, algo así como el edén hacker en la tierra.

Dado que para ese entonces mi economía estaba en franca decadencia y que el transporte público es atemorizántemente caro, mi hice fan a las CityBike que son bicicletas de uso público para la ciudad a bajo costo (incluso gratis, durante la primer hora). Esto me permitió disfrutar mucho de la ciudad y sus paisajes, ya que podía atravesarla en unos 40 minutos pedaleando.

Comí comida alemana, visité los puestos navideños, miré, saqué muchas fotos, dormí en un cuarto con una multitud de gente, tomé cerveza local en las happy hours, me sorprendí con comparaciones innecesarias y volví a casa.

UPDATE 27NOV 10:56 ARTS: El álbum de fotos completo aquí.

La estadía en la capital española fue particularmente social. La pasé con amigos, amigos de amigos, y otros niveles de indirección. Redescubrí la comida vegetariana gracias a Santiago y además cocinó y hospedó como solo los grandes anfitriones saben hacerlo.

Pero claro, tanta sociabilidad acortó las noches y llegué realmente cansado a la noche del domingo. Así fue como me quedé dormido el lunes y por poco no llego al aeropuerto a tiempo para irme a Viena. Por suerte, llegué con lo justo.

El único traspié fue desencontrarse con los Debianitas ibéricos (Holger, Amaya, Agi, etc) ya que los horarios del transporte de fin de semana me jugaron una mala pasada y llegué muy tarde. Tenía muchísimas ganas de verlos… una verdadera lástima.

BTW, me encontraba intentado reservar una cama en este hostel de Viena cuando, a la hora de pagar el adelanto, algo me dijo que podía ser peligroso hacerlo… adivina qué? :)

La pasé muy bien. Aproveché la oportunidad para ver a Dato, comer hasta el hartazgo y visitar la ciudad. Las fotos, en un par de horas.

UPDATE 27NOV 10:56 ARTS: El álbum de fotos completo aquí.

El hecho es que llegamos el viernes y salimos a recorrer la ciudad. Creo que nunca vi llover tanto en mi vida. A la mañana siguiente visitamos la Mezquita Hassan II. Casablanca no resultó ser tan turística como esperábamos, así que decidimos salir ese mismo día a Fez.

Fez es la primer ciudad imperial y antigua capital del reino. Su medina es la más conocida de Marruecos. Pasamos momentos muy interesantes y cumplió sobradamente con mis expectativas.

Realmente me gustaría dar más detalles (si, si… y fotos) pero en 4 horas tengo que levantarme para irme a Alicante.

UPDATE 4NOV 23:03 UTC: El álbum de fotos completo aquí.

Como dicen por ahí, he vivido a cuerpo de rey en La Coruña. Diego organizó un evento de altísima calidad y se preocupó por todos los detalles. Nos llevó de turismo (recuerden, las fotos en formato analógico tardarán semanas en estar on-line) y a comer a una increíble variedad de lugares. El pulpo, las croquetas y la carne a la piedra están haciendo de este viaje un recorrido culinario de lo más colorido.

El evento se realizó en la Universidad de A Coruña y tanto el rector como el resto del personal a sido sumamente hospitalario.

El único traspié fue el golpe de Chema, que después de todo se lo tenía casi que merecido :P.

UPDATE 27NOV 10:59 ARTS: El álbum de fotos completo aquí.

Después del Asegur@IT IV fuimos con el Chema a correr en kartings. Obviamente me pasó el trapo, pero aún así estuvo de pelos. Después fui a lo Santiago quien cocinó rico y me proveyó un sofá donde pasar la noche además de un lugar donde dejar mis cosas. Hacía mucho que nos veíamos y nos pusimos al día.

Me acosté algo afiebrado y con ciertos dolores estomacales. A la mañana siguiente salía hacia La Coruña. Al levantarme me sentí realmente enfermo. Con chuchos y eso. Me auto-diagnostiqué “sobredosis de trajín“.

Sobre llovido mojado, el viaje que tenía que llevarme una hora para llegar al aeropuerto me tomó dos y perdí el vuelo. Por suerte salía otro al rato y no tuve mayores dificultades. Aunque la gracia costó 50 pavos.

En el ajetreo de la salida dejé mi cámara de fotos olvidada en la casa, por lo que tuve que comprar una descartable analógica. La vieja escuela que le dicen. Así que habrá que esperar el revelado para ver que tal. Una experiencia surrealista esto de no saber como salen después de apretar el botón.

Felizmente estoy en La Coruña, en el noveno piso del Hotel Riazor con vista al Atlántico que impresiona. Después de un ibuprofeno 400 recetado por MSN de la salud me siento muchísimo mejor. Y aunque la idea de tener que levantarme en 4 horas y media para poder participar de las Jornadas de Seguridad Informática no me agrada mucho, no puedo quejarme.

Por cierto, si el viejo Walt se hubiese dedicado al transporte público en vez de a los dibujos animados, Madrid sería Disney World.

UPDATE 4NOV 23:03 UTC: El álbum de fotos completo aquí.

Recién (bueno.. fue hace 15 horas) salí de Puebla y ya estoy en el Centro de Formación en Tecnologías de la Información y Comunicaciones Madrid-Sur, donde después de literalmente haber atravesado la ciudad dí una charla en un evento organizados por MicroSoft (si, leyeron bien).

Los últimos días en Puebla – México fueron espectaculares. La he pasado de lo mejor en el ENLi. Muchas nuevas amistades, comidas y cerveza. Y no hubiese sido posible sin el esfuerzo sobrehumano (casi literalmente) de Clemente, Natalia, Gladis y Tlehuek. Muchísimas gracias :). Claro que también tengo que incluir a Tato, que me llevó y me trajo a todos lados además de darme techo y comida. Señores, ustedes saben que son la neta del planeta!

Mañana salgo para un congreso en La Coruña, que tengo que sacarme los modismos mexicanos a fuerza de gallegos. Pasaré el resto de la tarde vagando por Madrid y yendo a los lugares que me quieran llevar :)

¿Las fotos de Puebla? Tal vez mañana ya estén online, veré que puedo hacer esta noche. Claro, siempre que no vuelva a cualquier hora :P

UPDATE 4NOV 23:05 UTC: El álbum de fotos completo aquí.

¡Feliz día, Vieja!

Ayer fue el Día de la Madre en Argentina y corresponder saludarla en esta tradición de carteles cuando uno está lejos. Normalmente linkearía a su foto, pero a mi jefa (como dicen acá) no le gusta salir fotografiada :). Espero la hayas pasado lindo.

Subí algunas fotos (y con comentarios). Lo bueno de visitar el DF por segunda vez es que dejé bastante de lado lo turístico para dedicarme a lo social. Si, visité lugares clásicos, pero con baja prioridad, ya que pasé la mayor parte del tiempo con amigos y amigos de amigos. Otra cosa notable es que me reconcilié con la gastronomía local. Tato se pasó y me llevó a comer a muy buenos lugares. Supongo que las fotos contarán más historia que este post.

Una de las primeras cosas que hice fue dar una charla en la UNAM. Después fue una seguidilla de cafés, bares y amigos, hasta ahora que en Puebla tendré que dar una clases de SSH y seguridad.

La estoy pasando muy bien. Estoy conociendo a mucha gente y eso me es muy divertido…

UPDATE 4NOV 23:00 UTC: El álbum de fotos completo aquí.

Un par de semanas atrás, y por razones que no viene al caso, volví al tema desde cero. En este caso llegué a cosas medianamente concretas (dentro del delirio, claro). Así que aquí están algunos de los resultados que me gustaría compartir con ustedes.

Introducción:

Para aquellos que no tienen idea de los que estoy hablando, una breve introducción. En la famosa tira cómica se plantea una sencilla hipótesis que puede resumirse en la siguiente pregunta: “¿a partir de qué edad tus posibilidades de tener pareja empiezan a bajar?”

Para esto tengamos en cuenta dos elementos:

• A mayor edad, menores las posibilidades de soltería. Es decir, si una persona tiene 20 años es más probable de que esté soltera a que ocurra lo mismo con alguien de 40 años.
• A mayor edad, mayor es el rango de edad de los candidatos. Para poner un ejemplo, no es lo mismo que una persona de 40 salga con otra de 55 que un joven de 18 salga con alguien de 33. Es claro no hay edades para el amor, pero estamos hablando de posibles candidatos. Incluso hablamos de las posibilidades de tener cosas en común, de fijarse en el otro. Randall Munroe, el autor de la tira, propone una fórmula a la que llama “standard creepiness rule” (algo así como “regla estándar de viejo-verdez“, asquerosidad, o algo del estilo). Según él (y muchos parecen coincidir), las personas no salen con alguien que tenga la mitad de su edad más siete. Por ejemplo, en mi caso y con 27 años, mi cota inferior está dada por (27/2)+7=20.5, lo que parece bastante razonable o, al menos, políticamente correcto. Esta misma fórmula también sirve como cota superior, ya que uno es el límite inferior de alguien mayor. Para mi edad corresponde a 40 años porque 40=(27-7)*2. Nótese que para una persona de 20 años sus candidatos van entre 17 y 26, mientras que para alguien de 30 sus opciones van entre 22 y 46. Se supone que hay más personas entre el segundo rango que el primero, al menos hasta que la taza de mortalidad empieza a hacer de las suyas.

Primera aproximación:

Lo interesante de esta situación es la posibilidad de analizarla con datos reales. Son datos que existen, que se preguntan en los censos, pero que no siempre están disponibles para cualquier mortal. Por suerte el INDEC (Instituto Nacional de Estadística y Censos) provee la información de una forma que, aunque críptica para el recién llegado, resulta interesante cuando uno le encuentra la vuelta. Así que fue allí a donde fui, a buscar el universo de la población argentina.

Tomando toda la población mayor a 14 años, por edades simples y situación conyugal, calculé la cantidad de personas candidatas para cada edad, tiendo en cuenta los rangos de edades de esos posibles candidatos. Consideremos candidatos a las personas no casadas (solteros, divorciados, separados legales y viudos). Aquí el gráfico:

Nótese la forma de W-invertida de la gráfica. El máximo absoluto se encuentra a los 27 años. Es decir, en este momento me encuentro en el cenit de la soltería. Es interesante observar, que después del punto de inflexión de los 41, a los 47 vuelve a haber un máximo. Llamemos a este punto El club de la divorciadas ;). Las estadísticas demuestran que la frase “la vida da revancha” tiene fundamento estadístico.

Pero los valor absolutos siempre me parecieron lejanos. Soy de los que disfruta hablar en términos de porcentaje, de probabilidades, es decir, en términos relativos. Así fue como generé este segundo gráfico, con la intención de responder a la pregunta “Si conozco a 100 personas dentro de mi rango de edad, ¿cuántas estarán no-casadas?”:

Así, la cantidad de personas sin compromisos maritales dentro de un rango, sobre la cantidad total de personas de ese rango da un porcentaje de disponibilidad para cada edad. Acá la cosa se pone algo más pesimista, asegurando una clara caída hasta los 44 años, edad desde la cual se inicia un tímido y poco constante ascenso.

Segunda aproximación:

Esta visión es más egocentrista y está pensada en términos totalmente propios. En lo personal, estoy interesado por personas de género femenino, que no solo no estén casadas, sino que además no convivan en pareja (las considero personas en relaciones lo suficientemente fuertes como para quitarlas del conjunto de candidatas). A estas limitaciones le sumaré la restricción geográfica de que tenga residencia dentro de mi ciudad natal y alrededores, por lo que solo tendré en cuenta personas dentro de Buenos Aires y conurbano bonaerence. La ubicación geográfica puede acarrear cambios culturales que tengan efectos en la situación conyugal de las personas, por lo que considero importante tenerlo en cuenta.

Felizmente INDEC tiene una forma rápida de filtrar por estas restricciones, así que seleccioné las regiones geográficas y construí un filtro donde solo se refiera a personas de sexo femenino que no estén conviviendo con una pareja.

Nuevamente, las mismas gráficas:


Esta vez la W-invertida refleja un máximo absoluto muy diferente al inicial, 53 años. Tal vez se parezca más al modelo que, creo que intuitivamente, proponía Munroe. En el gráfico de porcentajes, la caída es mucho más brusca y profunda, perdiendo más de una decena de puntos en su mínimo, que además se da bastante antes, a los 38 años. La pendiente de ascenso es empinada y termina por provocar mayores valores hacia el final de la esperanza de vida que la aproximación anterior.

Suposiciones y limitaciones:

Estos análisis son un emporio de suposiciones y limitaciones. Después de todo, dicen que la estadística es la ciencia del prejuicio. Acá una lista de suposiciones que pueden ser ampliamente discutibles:

• Los homosexuales y religiosos existen en cantidades despreciables. Así que todo hombre soltero busca una mujer y viceversa.
• Las personas con alguna relación de pareja no convivientes son consideradas disponibles. El INDEC no pregunta acerca de estado amoroso de las personas.
• La cantidad de personas con n y medio edad es la mitad de la cantidad de personas con n edad.
• La fórmula “standard creepiness rule” es válida. Aunque posiblemente pueda repetirse la experiencia con fómulas que mejor se ajusten a cada cultura. En lo personal, creo tener poco en común con alguien de 40 años. De hecho, formulas independientes para el hombre y la mujer creo que serían más aplicables. Tal vez los hombres tiene un target más orientado a la izquierda del eje. Pero es medio que gusto de cada uno.
• La información del censo corresponde a datos de 2001. Extrapolarlos no tiene sentido, ya que insertaría muchísimo error. Lo dicho, lo interesante son los números en proporciones.
• La definición de “candidato” esta dada por edad, género, posición geográfica y estado conyugal. No se contempló gustos y esas cosas de esas que hace que las parejas funcionen… o eso dicen :P.

Todo lo necesario para repetir la experiencia se encuentra acá. Enjoy it … que yo voy a estar buscando a esa chica bajo la curva :)

click en la imagen para agrandar

• México DF, mx: 16 a 20 de Octubre
• Puebla, mx: 20 a 25 de Octubre
• Madrid, es: 26 a 28 de Octubre
• La Coruña, es: 29 a 31 de Octubre
• Casablanca, ma: 31 de Octubre a 3 de Noviembre
• Madrid, es: 3 a 5 de Noviembre
• Alicante, es: 5 a 7 de Noviembre
• Madrid, es: 8 a 10 de Noviembre
• Viena, at: 10 a 16 de Noviembre

Si estas en alguna de estas ciudades, no dudes en enviarme por un correo así nos vemos, intercambiamos firmas y tomamos una cervezas (en ese estricto orden).

Some weeks ago I started suspecting about the attack to PFS in SSL with EDH. The key point is: the key space is dependent of the PRNG state. The bug affects the initialization of the PRNG, but the random string has not a pattern by it self. If you ask for many random numbers to the PRNG, you gonna get numbers that differ among them, since they are the output of a hash function of them self. So each random number depends on, besides the PID, the state of the PRNG pool in the moment (in other words, amount of bytes that you already pull from the PRNG pool before)

The explained attack was based in a fixed list of private exponents (which are selected randomly during the DHE handshake), presupposing that all the application call RAND_bytes() the same number of times before get it. To make the list of exponent I ran the openssl s_client with all the possible PIDs, hoping that all the applications behaves the same way.

After more tests I notice that that was an overgeneralization. The proof is in the pudding: wget and cURL, two simple CLI file retrievers, gets different exponent between them, even running with the same PID.

I was working on this when I accidentally found a really nice Eric Rescorla’s post which is deeply related with this. The post goes further and analyzes the interaction between how Apache forks off and how it generates SSL handshakes.

So, I made lists of secret exponents for wget, curl, openssl s_client and openssl s_server with a modification version of libssl (appling this messy patch) and running scripts like this:

for i in $(seq $((2**15)));
do
  export MAGICPID=$i;
  LD_LIBRARY_PATH="openssl.broken/" LD_PRELOAD="./getpid.so" \
     wget --no-check-certificate https://localhost/ -q  -O /dev/null ;
  echo $i ;
done

As you can see, I used the HD Moore’s GetPID faker shared library and a normal local Apache with mod_ssl. The broken libssl (which is in .openssl.broken/) store up in /tmp/data.key a csv with command name, PID and all the DH components (g, x, y and p).

But this way is farly unconfortable for others SSL deamon servers. Have you got any better idea?

Maybe the LHC is robing the planet of angular momentum. Maybe having played with quantum gravity has unpredictable consequences. Who cares the reason, it’s my dream becoming true.

And you, haven’t got any plans for the extra day in October yet? Luckily it’s weekend.

Desde ya hace tiempo que tenía intensiones de irme de LiveJournal. No es que funcione mal. Es que simplemente tiene cosas que no me cuadran. Me la paso adaptándome a lo que puede darme (como el caso de hack para el bloguear en planet.debian.org) y tiene limitaciones de diseño. La publicidad que empezó a surgir a la derecha de la pantalla es la gota que derramó el vaso. No es solo antiestética, sino que si decido tener publicidad es porque espero cobrar por ella.

Dado que RaqLink puede proveerme un hosting gratuito y que otros amigos han ofrecido espacio y ancho de banda, decidí mudarme y tener mi propio Blog que dependa de mi mismo.

Así fue como me decidí por Wordpress. Es lindo, sencillo y flexible. Por otro lado, dudo de su seguridad. Y este último punto no es menor. Veremos que tal anda durante los próximos meses. Si da mucho problema… volará por otra opción. Escucho opciones.

Una cosa es segura. No más LiveJournal. Este blog (es decir, lbello.livejournal.com) deja de existir como tal. Puedes acceder al nuevo en www.lucianobello.com.ar. Todos los post antiguos está migrados. Incluso los comentarios. En la pasada se han perdido los tags y el threading de los comentarios. Los primeros irán emergiendo con el correr del tiempo y de mi ratos libros. El segundo está definitivamente perdido. Aunque los nuevos comentarios si pueden anidarse, los viejos han quedado planos.

Anoche me hallé frente a la biblioteca buscando nada en espacial. Necesitaba despejar la cabeza, olvidarme de los asuntos terrenales asociados con las complejas interacciones y relaciones humanas.

Fue en esa circunstancia que me encuentro con De los números y su historia de Isaac Asimov (que puede ser descargado desde aquí). Mientras lo hojeaba vi la formula de Leibniz, promediando el capítulo 6 e inmediatamente atrajo mi atención.

Después de expresar la serie, Asimov explica:

<<Ustedes podrán condenar mi falta de perseverancia, pero los invito a calcular la serie de Leibniz simplemente hasta donde la hemos escrito más arriba, es decir hasta 4/15. Incluso pueden enviarme una postal para darme el resultado. Si al terminar se sienten desilusionados al descubrir que su respuesta no está tan cerca de π como lo está el valor 355/113, no se den por vencidos. Sigan sumando términos. Sumen 4/17 al resultado anterior, luego resten 4/19, después sumen 4/21 y resten 4/23, etcétera. Pueden seguir hasta donde lo deseen, y si alguno de ustedes descubre cuántos términos se requieren para mejorar el valor 355/113, escríbanme unas líneas y no dejen de decírmelo. >>
En efecto, la fracción 355/113 parece ser la forma que mejor balancea precisión y simpleza a la hora de arrimarse a π desde lo números racionales, quedando a solo 2.66764189 × 10^-7 de distancia. Mucha gente también utiliza la relación 22/7, aunque con menor precisión.

El desafío propuesto parecía interesante. Resultaba una entretenida oportunidad de jugar con floats en python, que siempre a resultado ser bastante tricky (además de ser un excusa para despejar la cabeza).

Empecé por escribir una versión relativamente obvia de la solución, que me dijo que al llegar al cociente 7497257 (n=3748629 en la sumatoria de Leibniz) el valor de la sumatoria estaría en 3.14159292035, quedando a 2.66764081935 x 10^-7 de π. Dado que los decimales (floats, representaciones de punto flotante) en Python no se comportan de forma intuitiva, decidí hacer una segunda versión.

En esta versión traté de mudarme al predecible mundo de los enteros. Y obtuve valores levemente distintos. En el cociente 7497247 (n=3748624) se obtiene la primera aproximación mejor que 355/113, que es 3.141592386825… (a 2.66764162 x 10^-7 de π)

Tampoco estoy seguro de que la segunda solución sea la correcta. Me pregunto cual será el promedio de valores en las postales que la viuda Asimov seguramente guarda en una caja de zapatos.

El hecho es que fue entretenido intentarlo y tal vez continúe con la experiencia la próxima vez que quiera desconectarme del mundo. Después de todo, como dijo el mismo Asimov: “Los hombres que se acostumbran a preocuparse por las necesidades de unas máquinas, se vuelven insensibles respecto a las necesidades de los hombres”, y hay veces que volverse insensible se ve, erróneamente, como una propuesta seductora. Hay gente que ahoga penas en alcohol, mucho más sano es ahogarlas matemáticas…

BTW, ya que tenía los datos armé un plot de los primeros 1000 valores de la serie para mostrar de forma gráfica la convergencia.

Ahora es momento de volver a la realidad, que buena o mala, irracional o racional (seguro que real al menos), inexacta o precisa, después de este recreo metal ya no se ve tan mal :-)

UPDATE: Sat, 20 Sep 2008 22:18:18 -0300: La solución estaba realmente cerca (en mi misma ciudad) y se llama Facundo, quien dejó un comentario. Bah.. se llama decimal, pero existe gracias a Facundo :). La biblioteca decimal permite manejar de forma exacta una cantidad arbitraria de decimales. La respuesta correcta es…….. (redoblantes de suspenso) con n=(7497259+1)/2! aproxima a Π en 3.141592386825668744985771256 a solo 2.66764125 × 10^-7 del número irracional. El nuevo script.. aquí.

Gracias Facundo, gracias decimal, gracias Asimov, gracias Leibniz….

Obviamente lo único que soy capaz de entender es mi propio nombre, pero lo que cuenta es la intención :). Gracias Koen por el detalle, se que entenderás este post tanto como yo la nota :P.

Y hablando de correo físicamente transportado (aunque sin relación con lo anterior, que fue totalmente gratis), uno se da cuenta que vive en las antípodas del mundo donde le gustaría estar cuando recibe facturas de encomienda como esta: