El siguiente código que hemos desarrollado consta de una clase "Partida", que a su vez contiene una serie de métodos que simulan un juego. El juego consta de 4 participantes que a lo largo de 6 rondas (ver "emparejamientos") muestran y observan las señales de sus parejas. Las señales que son mostradas por cada participante son almacenadas en su memoria de señales mostradas y las señales que son observadas son almacenadas en su memoria de señales observadas. Con una ecuación calculamos la probabilidad de que en cada ronda un participante genere una señal, dependiendo de la memoria y de otras variables.
En el código actual, para calcular la probabilidad de producción de una señal en una ronda determinada (with_b) por un jugador determinado, se tiene en cuenta el cómputo total de señales almacenadas en las memoria desde el inicio del juego. Es decir, las señales que se han ido almacenando desde el inicio del juego determinan en cierta medida lo que se va a producir.
Me gustaría ser flexible en este aspecto. ¿Cómo podría implementar en el código la opción de que en cada ronda solo se tenga en cuenta la memoria almacenada de las últimas x generaciones?
from random import random, sample
from bisect import bisect
import csv
class Partida():
def __init__(self, jugadores, emparejamientos, senales, s, b, x, m):
self.emparejamientos = emparejamientos
self.senales = senales
self.s = s
self.b = b
self.x = x
self.m = m
self.jugadores = {nombre: Partida.Jugador(senales)
for pareja in emparejamientos[0]
for nombre in pareja}
self.memoria = list()
def generar_senales(self):
def with_b(muestra, observa, s, r):
if not (muestra == observa == 0):
result = ((0.98) * (1.0 - self.b) * (1.0 - self.x) * muestra/r) + ((0.98) * (1.0 - self.b) * (self.x) * observa/r) + ((0.98) * self.b * s) + ((self.m / 8))
else:
result = ((0.98) * (1.0 - 0) * (1.0 - self.x) * muestra/r) + ((0.98) * (1.0 - 0) * (self.x) * observa/r) + ((0.98) * 0 * s) + ((self.m / 8))
return result
def choice(opciones, probs):
probAcumuladas = list()
aux = 0
for p in probs:
aux += p
probAcumuladas.append(aux)
r = random() * probAcumuladas[-1]
op = bisect(probAcumuladas, r)
return opciones[op]
yield dict(zip(self.jugadores.keys(), self.senales))
r = 1
while True:
eleccs = dict.fromkeys(self.jugadores.keys())
for nombre, inst in self.jugadores.items():
probs = [with_b(inst.mem_mostradas[op], inst.men_observadas[op], self.s[indx], r)
for indx, op in enumerate(self.senales)]
eleccs[nombre] = choice(self.senales, probs)
r += 1
yield eleccs
def jugar(self):
gen_sens = self.generar_senales()
for n, ronda in enumerate(self.emparejamientos):
senales = next(gen_sens)
self.memoria.append(senales)
for jugador1, jugador2 in ronda:
self.jugadores[jugador1].men_observadas[senales[jugador2]] += 1
self.jugadores[jugador2].men_observadas[senales[jugador1]] += 1
self.jugadores[jugador1].mem_mostradas[senales[jugador1]] += 1
self.jugadores[jugador2].mem_mostradas[senales[jugador2]] += 1
class Jugador():
def __init__(self, senales):
self.mem_mostradas = {senal: 0 for senal in senales}
self.men_observadas = {senal: 0 for senal in senales}
def main():
jugadores = [1, 2, 3, 4]
senales = ['Senal 1', 'Senal 2', 'Senal 3', 'Senal 4']
emparejamientos = [[(1,2),(3,4)],
[(1,3),(2,4)],
[(1,4),(2,3)],
[(1,2),(3,4)],
[(1,3),(2,4)],
[(1,4),(2,3)]]
s=[1,0,0,0]
muestras = [{'b':0.0, 'x':0.5, 'm':0.02} for _ in range(1000)]
simulaciones = 100
estadisticas = {jugador:{muestra:{senal:[0 for ronda in range(1, len(emparejamientos)+1)]
for senal in senales}
for muestra in range(len(muestras))}
for jugador in jugadores}
for mu in range(len(muestras)):
for _ in range(simulaciones):
juego = Partida(jugadores, emparejamientos, senales, s, muestras[mu]['b'],muestras[mu]['x'], muestras[mu]['m'])
juego.jugar()
for n, ronda in enumerate(juego.memoria):
for jugador, senal in ronda.items():
estadisticas[jugador][mu][senal][n] += 1
with open('salidav4.csv', 'wb') as csvfile:
writer =csv.writer(csvfile, delimiter=',',
quotechar='"', quoting=csv.QUOTE_MINIMAL)
writer.writerow(['Muestra' ,'Jugador', 'Ronda', 'b', 'x', 'm'] + senales)
for jugador in jugadores:
for mu in range(len(muestras)):
for ronda in range(1, len(emparejamientos)+1):
aux = [estadisticas[jugador][mu][senal][ronda-1] for senal in senales]
writer.writerow([mu+1, jugador, ronda, muestras[mu]['b'], muestras[mu]['x'], muestras[mu]['m']]+aux)
if __name__ == '__main__':
main()