To run the simulation start the file init.hoc
This will simulate the network model and reproduce the results
associated with Fig. 6 in:
Conde-Sousa E and Aguiar P, "A working memory model for serial order
that stores information in the intrinsic excitability properties of
neurons", JComputNeuroscience 2013
This is a network model with several neurons: the simulation therefore
takes a while(*)
A raster plot is automatically generated in NEURON after the
simulation. Alternativelly you can run the MATLAB script plotspikes.m
(or plotspikes_for_all_neurons.m) after the NEURON simulation is
completed to generate the raster in MATLAB.
In the raster plot:
- neurons with id from 0 to N_PRINCIPAL_NEURONS-1 are inhibitory
interneurons
- neurons with id from N_PRINCIPAL_NEURONS to 2*N_PRINCIPAL_NEURONS-1
are principal neurons
- all other are gate interneurons
(*)A smaller version of the network can be explored (lines 42-50 in
WMSeqLearn.hoc):
N_PRINCIPAL_NEURONS = 20
GATES_SAMPLING_FRACTION = 1
CONN_RATE = 0
N_PATTERNS = 5
PATTERN_SIZE = 3