¿Extorsión en el Dilema del Prisionero?

La verdad, los tengo encima de la mesa y aun no los he digerido, el PNAS de Press y Dyson y el comentario de Stewart y Plotkin.

Ayer los mencionarón @elaragon y @Ciudadano_Zer0 en el twitter, haciendo referencia a un apunte en technology review:
http://www.technologyreview.com/view/428920/the-emerging-revolution-in-game-theory/
Tambien ha hecho un comentario Lubos Motl:
http://motls.blogspot.com.es/2012/08/freeman-dyson-and-william-press.html

Al parecer los comentarios van sobre  un nuevo preprint (Adami-Hintze) en el que se arguye que estos algoritmos de Press y Dyson no son estables evolutivamente… lo que se me hace un argumento raro, de buenas a primeras, porque se trataba precisamente de que eran algoritmos que vencían a los evolutivos.

Dyson siempre me ha resultado de confianza, incluso tengo sus apuntes de mecanica cuantica relativista del curso de Les Houches (de los sesenta). Asi que habrá que darle un repaso.

Desigualdad en el reparto.

Suponed que tenemos que repartir, digamos, el PIB de España entre todos los españoles, pero que se mete el azar por medio hay que hacerlo sorteando cada euro. En principio no pasa nada, cuando el numero de sorteos sea muchisimo más grande que el de participantes, tocaremos a todos. Por ejemplo, he hecho un intento de repartir 250000 tiradas entre 10000 y sale gente con 1.7 veces la media, pero cuando he repartido 10 millones de tiradas ya apenas habia alguno que otro con 1.1 veces la media.

Ahora, ¿y cuando dinero llama a dinero? ¿qué pasa si el sorteo lo hacemos de forma que cada cual lleve tantos tickets como dinero ha ganado ya?

Esto lo podemos implementar en un programita en python. Primero repartimos una parte a cada uno de los participantes, y para el resto tiramos de loteria de la forma proporcional a la riqueza que lleva cada uno ganada.

from random import randint
from array import array
import sys
total=10000
#reparto=20000
reparto=int(sys.argv[1])
cantidad=array('l',(1 for x in xrange(total)))
for i in xrange(reparto-total):
   concede=randint(1,total)
   x, ticks =-1,0
   while ticks < concede:
      x+=1
      ticks+=cantidad[x]
   cantidad[x]+=1
   total+=1
   #burbuja
   while (x > 0) and (cantidad[x] > cantidad[x-1]):
      cantidad[x],cantidad[x-1]= cantidad[x-1],cantidad[x]
      x=x-1
print cantidad[:20]
ricos=sum(cantidad[:len(cantidad)/100])
print "hemos repartido ", total , " entre ", len(cantidad),
print " (media=",1.0*total/len(cantidad), ")"
print  "el 1% posee el ", 100.0*ricos/total, "% del reparto"
print "el mas rico posee ", 1.0*cantidad[0]*len(cantidad)/total,
print "veces mas que la media"

Y si por ejemplo ejecutamos

python reparte.py 500000

Nos va a salir algo bastante más injusto:

array('l', [595, 469, 450, 440, 403, 401, 358, 355, 355, 352, 347, 343, 339, 334, 334, 328, 327, 321, 319, 317])
hemos repartido 500000 entre 10000 (media= 50.0 )
el 1% posee el 5.6656 % del reparto
el mas rico posee 11.9 veces mas que la media

Para este mismo reparto, el aleatorio sin sesgo saca que el 1% posee el 1.40 % de lo repartido y el más rico posee tan solo vez y media lo que el promedio.

Vamos, que en efecto lo de hacer que dinero llame a dinero genera desigualdades, eso estaba claro… Pero aún asi, no acaban de salirme las típicas brutalidades en las que el 1% posee el 25% de la riqueza nacional y tal. Sacar diez veces la media salarial es brutal, pero no te convierte en Amancio Ortega, y si mirais el array vereis que el que esta en el puesto 15 ya ha bajado a la mitad de beneficios.

Tertulias de Politica

No había previsto en este blog una categoría para política, aunque siempre ha algo en casi todas mis paginas. Incluso en el 95 intenté crear un archivo de documentos políticos. Posiblemente algo de influencia mutua habrá, entre los temas que me han interesado en investigación y las cuestiones políticas. Asi que por lo menos conviene hacer un post para encuadrar el asunto. Más adelante ya veremos.

Y el asunto es que buscando en las news (los foros, para novatos) y listas de email asuntos sobre anarquismo, acabe conociendo a unos cuantos matemáticos españoles próximos a la Idea. En un foro se planeaba una lectura semanal del libro “De Dios”, de Agustin Garcia-Calvo. Y asi, cuando al acabar la tesis me encontré en Madrid haciendo la mili (sí, un poco de paradoja, aquí) en la Marina Española, localizé la Tertulia Politica de Agustín y comenzé a aprovechar algún miercoles en el Ateneo. Un año despues volvería a la Capital de la Esencia y Potencia, esta vez para trabajar en Teléfonica, y pasé a ser asiduo. Continue reading →

Los tiempos de la Lattice

Lamentablemente estoy ya muy desconectado de la investigación de teoria de campos en la lattice. Pero es interesante ver los dos primeros articulos en los que colabore como programador a veces, administrador otras, y hasta soldador y montador si se terciaba. En el Arxiv tenemos dos:

[hep-lat/9210014] The U(1)-Higgs Model: Critical Behaviour in the Confinig-Higgs region

[hep-lat/9302007] Instanton-like Contributions to the Dynamics of Yang-Mills Fields on the Twisted Torus

que quizás son representativos de otras lineas que seguí luego, aunque en el mundo analitico: los instantones con Casahorran y con mi director de tesis, Boya. Y el dichoso Higgs. En esos tiempos J. L. Alonso le estuvo dando vueltas, como especulación, a las consecuencias de la existencia de soluciones antiferromagneticas en el campo de Higgs en la lattice, y si eso podia explicar la diferencia entre las masas entre fermiones y bosones masivos, o quizas entre fermiones cargados y neutrinos, ya no lo recuerdo bien 🙂 . Pero la idea de tener que llevar al limite una solución antiferro conectaba bastante bien, en mi fuero interno, con la separacion del Higgs en dos capas, en una dimensión discreta a la que apuntaban los modelos de Connes-Lott.

Tambien en aquella epoca nos metimos en el asunto de los ordenadores de proposito especifico, colaborado con el grupo APE de Roma. Preparamos varias maquinas basadas en Transputers, en particular una que llamabamos RTN, de 64 procesadores. Mi primera misión en el extranjero fue llegarme a Roma I para traer las placas base ya terminadas y verificadas por nuestro ingeniero, Jarda, que compartiamos con APE. Al pasar la frontera, el guardia civil del scanner me preguntó si llevaba “un piano electronico” en la maleta. Naturalmente le conteste que sí.