I reimplemented sim_background in brian2.
This was well possible with the information from the publication 'Zenke, Hennequin, Gerstner et al in PlosCB 2013' together with the code 'sim_background' provided online. Thanks for providing the code online!
While the manuscript was clear about the values for inhibitory and excitatory resting potentials, as well as the excitatory membrane time constant, I could not find where these were set in the code, or whether Auryn provides default values for those? The documentation nicely mentions default values for different synaptic time constants (AMPA, GABA, NMDA).
The simulation code revealed that membrane voltages were initialised normally distributed around -60mV with a standard deviation of 10mV.
I get an asynchronous irregular regime, as desired, with an alpha (AMPA to NMDA ratio) of .5 as in the manuscript.
However, with the random initialisation of membrane voltages I expected to get into the regime a bit quicker? So there might still be discrepancies between the two implementations.
As zenke pointed out himself, in Eq. (13) and (14) of the manuscript, only the conductance term (g_inh and g_ampa, respectively) is supposed to be divided by the synaptic time constant (tau_inh and tau_ampa, respectively).
Here's the code:
https://github.com/k47h4/scripts/blob/m ... ced_net.py
Discussions on development, Auryn's inner workings, adding new features and how to change things at the core
1 post • Page 1 of 1
Who is online
Users browsing this forum: No registered users and 1 guest