Network of HVC Neurons with Heterogeneous Dopamine Modulations

Prerequisites

The simulations were run with the following packages and versions:

numpy 2.3.3
matplotlib 3.10.1
brian2 2.9.0

Running the Simulations

You can check out the demo by running python3 network_rand_conn_demo.py which contains a few pre-specified parameters and will generate raster plots and firing rate plots.

To reproduce the figure in the paper, run

python3 HVC_single_neurons_DA.py
run_sims.sh

The first line will generate the single neuron results. The second line which will iterate over different parameters and automatically generate .svg images. Simulation data will be stored in ./results.

To try other values of parameters, run

python3 network_rand_conn.py <DA level> <HVC_RA D1/D2 ratio> <HVC_RA D2 / HVC_Int D2 ratio>

with <...> replaced by your specified values for these three parameters.