Hace algunos días me pedían en los comentarios de otra pregunta que subiera mi código y reformulara la pregunta. Aquí va.
El siguiente código contiene una serie de variables y funciones, que definen un juego de señales según un modelo determinado. El output nos da las señales que muestra cada jugador en cada ronda.
Necesito saber cómo:
Simular este juego x veces (por ejemplo 1000). He conseguido hacerlo con execfile("Mimodulo.py"), pero no es útil para ir almacenando los resultados.
Almacenar en un archivo csv (filas (señales) y columnas (jugadores)), los resultados. Es decir, que el archivo csv contenga al final de las 1000 simulaciones la cantidad de cada tipo de señal que cada jugador ha producido en cada ronda. Esto es útil para hacer luego el análisis estadístico. Tan solo he conseguido abrir y cerrar un archivo CSV, pero me temo que estoy bloqueado.
Aquí está el código:
import random
##### Variables ######
emparejamientos= ([[1,2],[3,4], #ronda 1 (participante 1 juega con 2 y 3 con 4)
[1,3],[2,4], #ronda 2 (1 con 3 y 2 con 4)
[1,4],[2,3]]) #ronda 3 (1 con 4 y 2 con 3)
#Parejas de cada participante en cada ronda (empezando por la ronda 1).
pareja_part1= [2,3,4]
pareja_part2= [1,4,3]
pareja_part3= [4,1,2]
pareja_part4= [3,2,1]
#Parametros del modelo
s1=1
s2=0
s3=0
s4=0
b=0.5
x=0.5
m=0.02
round=[0,1,2,3,4,5,6] #rounds
#Senales
S1="S1"
S2="S2"
S3="S3"
S4="S4"
opciones = [S1, S2, S3, S4]
#Diccionarios de memoria
Muestra_part1 = {"S1":0,"S2":0,"S3":0,"S4":0}
Observa_part1 = {"S1":0,"S2":0,"S3":0,"S4":0}
Muestra_part2 = {"S1":0,"S2":0,"S3":0,"S4":0}
Observa_part2 = {"S1":0,"S2":0,"S3":0,"S4":0}
Muestra_part3 = {"S1":0,"S2":0,"S3":0,"S4":0}
Observa_part3 = {"S1":0,"S2":0,"S3":0,"S4":0}
Muestra_part4 = {"S1":0,"S2":0,"S3":0,"S4":0}
Observa_part4 = {"S1":0,"S2":0,"S3":0,"S4":0}
##### Funciones #####
#Almacen de senales mostradas en el diccionario de memoria Muestra_partx
def mem_mostradas(e, memoria):
if e == S1:
memoria[S1] +=1
if e == S2:
memoria[S2] +=1
if e == S3:
memoria[S3] +=1
if e == S4:
memoria[S4] +=1
return memoria
#Almacen de senales observadas en el diccionario de memoria Observa_partx
def mem_observadas(pareja_part, memoria):
if pareja_part == 1:
memoria[eleccion1] +=1
if pareja_part == 2:
memoria[eleccion2] +=1
if pareja_part == 3:
memoria[eleccion3] +=1
if pareja_part == 4:
memoria[eleccion4] +=1
return memoria
#Probabilidad para cada participante de mostrar una senal determinada en una ronda
def with_b(muestra, observa, s, r):
result = ((0.98) * (1.0 - 0) * (1.0 - x) * muestra/r) + ((0.98) * (1.0 - 0) * (x) * observa/r) + ((0.98) * 0 * s) + ((m / 8))
if not (muestra == observa == 0):
result = ((0.98) * (1.0 - b) * (1.0 - x) * muestra/r) + ((0.98) * (1.0 - b) * (x) * observa/r) + ((0.98) * b * s) + ((m / 8))
return result
#Funcion que arroja la opcion que cada participante elige para mostrar
from random import random
from bisect import bisect
def choice(opciones, probs):
probAcumuladas=[]
aux = 0
for p in probs:
aux += p
probAcumuladas.append(aux)
r = random() * probAcumuladas[-1]
op = bisect(probAcumuladas, r)
return opciones[op]
##### Juego #####
#RONDA 1
#Senales mostradas por cada participante en ronda 1
eleccion1=S1
eleccion2=S2
eleccion3=S3
eleccion4=S4
print "Senales mostradas en la generacion 1 \n jugador1: {0}, jugador2: {1}, jugador3: {2}, jugador4: {3}".format(S1, S2, S3, S4)
#Almacen de senales mostradas en el diccionario (memoria de senales mostradas)
Muestra_part1=mem_mostradas(eleccion1,Muestra_part1)
Muestra_part2=mem_mostradas(eleccion2,Muestra_part2)
Muestra_part3=mem_mostradas(eleccion3,Muestra_part3)
Muestra_part4=mem_mostradas(eleccion4,Muestra_part4)
#Almacen de senales observadas en el diccionario (memoria de senales observadas)
Observa_part1=mem_observadas(pareja_part1[0], Observa_part1)
Observa_part2=mem_observadas(pareja_part2[0], Observa_part2)
Observa_part3=mem_observadas(pareja_part3[0], Observa_part3)
Observa_part4=mem_observadas(pareja_part4[0], Observa_part4)
#RONDA 2
#Probabilidad para cada participante de mostrar una senal determinada en 2a ronda
Prob_S1_part1 = with_b(Muestra_part1[S1], Observa_part1[S1], s1, round[1])
Prob_S2_part1 = with_b(Muestra_part1[S2], Observa_part1[S2], s2, round[1])
Prob_S3_part1 = with_b(Muestra_part1[S3], Observa_part1[S3], s3, round[1])
Prob_S4_part1 = with_b(Muestra_part1[S4], Observa_part1[S4], s4, round[1])
Prob_S1_part2 = with_b(Muestra_part2[S1], Observa_part2[S1], s1, round[1])
Prob_S2_part2 = with_b(Muestra_part2[S2], Observa_part2[S2], s2, round[1])
Prob_S3_part2 = with_b(Muestra_part2[S3], Observa_part2[S3], s3, round[1])
Prob_S4_part2 = with_b(Muestra_part2[S4], Observa_part2[S4], s4, round[1])
Prob_S1_part3 = with_b(Muestra_part3[S1], Observa_part3[S1], s1, round[1])
Prob_S2_part3 = with_b(Muestra_part3[S2], Observa_part3[S2], s2, round[1])
Prob_S3_part3 = with_b(Muestra_part3[S3], Observa_part3[S3], s3, round[1])
Prob_S4_part3 = with_b(Muestra_part3[S4], Observa_part3[S4], s4, round[1])
Prob_S1_part4 = with_b(Muestra_part4[S1], Observa_part4[S1], s1, round[1])
Prob_S2_part4 = with_b(Muestra_part4[S2], Observa_part4[S2], s2, round[1])
Prob_S3_part4 = with_b(Muestra_part4[S3], Observa_part4[S3], s3, round[1])
Prob_S4_part4 = with_b(Muestra_part4[S4], Observa_part4[S4], s4, round[1])
probabilidades1 = [Prob_S1_part1, Prob_S2_part1, Prob_S3_part1, Prob_S4_part1]
probabilidades2 = [Prob_S1_part2, Prob_S2_part2, Prob_S3_part2, Prob_S4_part2]
probabilidades3 = [Prob_S1_part3, Prob_S2_part3, Prob_S3_part3, Prob_S4_part3]
probabilidades4 = [Prob_S1_part4, Prob_S2_part4, Prob_S3_part4, Prob_S4_part4]
#Senales que muestran los participantes en la RONDA 2. Asignadas segun las probabilidades dadas por la ecuacion de memorias
eleccion1 = choice(opciones, probabilidades1) #senal mostrada por participante 1 en ronda 2
eleccion2 = choice(opciones, probabilidades2) #senal mostrada por participante 2 en ronda 2
eleccion3 = choice(opciones, probabilidades3) #senal mostrada por participante 3 en ronda 2
eleccion4 = choice(opciones, probabilidades4) #senal mostrada por participante 4 en ronda 2
print "Senales mostradas en la generacion 2 \n jugador1: {0}, jugador2: {1}, jugador3: {2}, jugador4: {3}".format(eleccion1, eleccion2, eleccion3, eleccion4)
#Almacen de senales mostradas en el diccionario (memoria de senales mostradas)
Muestra_part1=mem_mostradas(eleccion1,Muestra_part1)
Muestra_part2=mem_mostradas(eleccion2,Muestra_part2)
Muestra_part3=mem_mostradas(eleccion3,Muestra_part3)
Muestra_part4=mem_mostradas(eleccion4,Muestra_part4)
#Almacen de senales observadas en el diccionario (memoria de senales observadas)
Observa_part1=mem_observadas(pareja_part1[1], Observa_part1)
Observa_part2=mem_observadas(pareja_part2[1], Observa_part2)
Observa_part3=mem_observadas(pareja_part3[1], Observa_part3)
Observa_part4=mem_observadas(pareja_part4[1], Observa_part4)
#RONDA 3
#Probabilidad para cada participante de mostrar una senal determinada en 3a ronda
Prob_S1_part1 = with_b(Muestra_part1[S1], Observa_part1[S1], s1, round[2])
Prob_S2_part1 = with_b(Muestra_part1[S2], Observa_part1[S2], s2, round[2])
Prob_S3_part1 = with_b(Muestra_part1[S3], Observa_part1[S3], s3, round[2])
Prob_S4_part1 = with_b(Muestra_part1[S4], Observa_part1[S4], s4, round[2])
Prob_S1_part2 = with_b(Muestra_part2[S1], Observa_part2[S1], s1, round[2])
Prob_S2_part2 = with_b(Muestra_part2[S2], Observa_part2[S2], s2, round[2])
Prob_S3_part2 = with_b(Muestra_part2[S3], Observa_part2[S3], s3, round[2])
Prob_S4_part2 = with_b(Muestra_part2[S4], Observa_part2[S4], s4, round[2])
Prob_S1_part3 = with_b(Muestra_part3[S1], Observa_part3[S1], s1, round[2])
Prob_S2_part3 = with_b(Muestra_part3[S2], Observa_part3[S2], s2, round[2])
Prob_S3_part3 = with_b(Muestra_part3[S3], Observa_part3[S3], s3, round[2])
Prob_S4_part3 = with_b(Muestra_part3[S4], Observa_part3[S4], s4, round[2])
Prob_S1_part4 = with_b(Muestra_part4[S1], Observa_part4[S1], s1, round[2])
Prob_S2_part4 = with_b(Muestra_part4[S2], Observa_part4[S2], s2, round[2])
Prob_S3_part4 = with_b(Muestra_part4[S3], Observa_part4[S3], s3, round[2])
Prob_S4_part4 = with_b(Muestra_part4[S4], Observa_part4[S4], s4, round[2])
probabilidades1 = [Prob_S1_part1, Prob_S2_part1, Prob_S3_part1, Prob_S4_part1]
probabilidades2 = [Prob_S1_part2, Prob_S2_part2, Prob_S3_part2, Prob_S4_part2]
probabilidades3 = [Prob_S1_part3, Prob_S2_part3, Prob_S3_part3, Prob_S4_part3]
probabilidades4 = [Prob_S1_part4, Prob_S2_part4, Prob_S3_part4, Prob_S4_part4]
#Senales que muestran los participantes en la RONDA 3. Asignadas segun las probabilidades dadas por la ecuacion de memorias
eleccion1 = choice(opciones, probabilidades1) #senal mostrada por participante 1 en ronda 3
eleccion2 = choice(opciones, probabilidades2) #senal mostrada por participante 2 en ronda 3
eleccion3 = choice(opciones, probabilidades3) #senal mostrada por participante 3 en ronda 3
eleccion4 = choice(opciones, probabilidades4) #senal mostrada por participante 4 en ronda 3
print "Senales mostradas en la generacion 3 \n jugador1: {0}, jugador2: {1}, jugador3: {2}, jugador4: {3}".format(eleccion1, eleccion2, eleccion3, eleccion4)
#Almacen de senales mostradas en el diccionario (memoria de senales mostradas)
Muestra_part1=mem_mostradas(eleccion1,Muestra_part1)
Muestra_part2=mem_mostradas(eleccion2,Muestra_part2)
Muestra_part3=mem_mostradas(eleccion3,Muestra_part3)
Muestra_part4=mem_mostradas(eleccion4,Muestra_part4)
#Almacen de senales observadas en el diccionario (memoria de senales observadas)
Observa_part1=mem_observadas(pareja_part1[2], Observa_part1)
Observa_part2=mem_observadas(pareja_part2[2], Observa_part2)
Observa_part3=mem_observadas(pareja_part3[2], Observa_part3)
Observa_part4=mem_observadas(pareja_part4[2], Observa_part4)
El output en una simulación ejemplo es el siguiente:
Senales mostradas en la generacion 1
jugador1: S1, jugador2: S2, jugador3: S3, jugador4: S4
Senales mostradas en la generacion 2
jugador1: S1, jugador2: S1, jugador3: S3, jugador4: S4
Senales mostradas en la generacion 3
jugador1: S2, jugador2: S1, jugador3: S1, jugador4: S1
Por lo tanto, se busca ir almacenando en un archivo CSV, durante 1000 simulaciones, las señales que cada jugador ha generando en cada ronda.
Un ejemplo de cómo sería útil almacenar los resultados en un CSV para el posterior análisis estadístico es el siguiente. En el ejemplo, las columnas RONDAX muestran el número de veces que se ha elegido cada señal (los números de esta imagen no corresponden a ninguna simulación con el programa, y por eso los sumatorios son inconsistentes, pues se trata solo de un ejemplo para ilustrar el formato de la tabla).
Edición (11.01.2016):
Partiendo de la última versión del código compartido por @FJSevilla, me gustaría poder ir un poco más allá. De cara al análisis estadístico de los datos sería interesante obtener no ya 1000 simulaciones, sino x muestras de 1000 simulaciones.
Utilizando el mismo código, la tabla de resultados que en este caso se pretende obtener es como la siguiente (ejemplificando cómo sería un CSV que almacenararía los resultados de 3 muestras de 1000 simulaciones del juego completo):
Donde en "muestra" el número 1 hace referencia a las primeras 1000 simulaciones, 2 hace referencia a las segundas 1000 simulaciones y así sucesivamente. Donde "ronda" hace referencia al número de ronda dentro del juego. Donde "jugador" hace referencia al número de jugador.
Nótese que en este caso se incluyen los parámetros b, x y m en la tabla, y que los números en las columnas "Senalx" hacen referencia a las veces que cada señal ha sido elegida en cada ronda.
Es una configuración del orden de los datos diferente a la que teníamos, pero que facilitaría la comparación estadística. ¿Alguna mano amiga con esto?
while
ofor
acompañado de unrange()
y cada vez que termine una iteración guardas los resultados en el csv, que puedes hacer conwith open('archivo.csv', 'a') as escribir
y luego almacenas tu resultado en un string formateado con.format()
y ya únicamente lo escribes así:escribir.write(string)
eso me parece podría funcionar.