In [ ]:

    

# this function code has to be debugged further with matrix dimension matching
import numpy as np
def generate_poisson_spikes(spikeMat, imgPara):
    meanFiringRate = np.mean(spikeMat, 1)
    [numNeuron, NT]= spikeMat.shape # NT is number of frames through 20 trials
    poissSpikes = np.zeros(numNeuron, NT)
    for i in range(numNeuron)
        poissSpikes[i,:] = np.random.poisson(meanFiringRate[i], (1, NT))
    return poissSpikes