This matlab code was contributed by S Zannone.
Instructions to run the model:
The file seq_neuromod_plasticity.m is needed to run the plasticity rules called +ACh and -ACh (the presence or absence of acetylcholine can be decided with the flag ACh_flag). This is needed for Figures 2, 3 and 4. This is the bit of the code that partially overlaps with the elife paper:
Brzosko Z, Zannone S, Schultz W, Clopath C, Paulsen O (2017) Sequential neuromodulation of Hebbian plasticity offers mechanism for effective reward-based navigation. Elife
http://dx.doi.org/10.7554/eLife.27756
Note: This model accession 245018 archive extends ModelDB
accession number 233396
by including the new files dynamic_reward_signal.m and
negative_feedback_signal.m. The seq_neuromod_plasticity.m scripts are
essentially the same between the two modeldb entries, with color and
content edits to the plotting commands.
The same file (seq_neuromod_plasticity.m), however, can be also used to run simulations of the r-STDP rule by setting ACh_flag = 0 and changing the STDP parameters (A_pre_post = pre-post window amplitude, A_post_pre = post-pre window amplitude, tau_pre_post = pre-post window time constant, tau_post_pre = post-pre window time constan). This is needed to reproduce Figure 5.
The file dynamic_reward_signal.m was used for Figure 6 and negative_feedback_signal.m is for Figure 7.
Here is a sample screenshot from running negative_feedback_signal which displays trials 1-40 results as they run: