Short-term persistent spiking buffer model of entorhinal cortex described by Randal A. Koene and Michael E. Hasselmo
octave runSim.m
Cell 1
spiketime(ms) 102 212 322 432 543 658 778 902 1027 1152 1277 1402
ISI(ms) 110 110 110 111 115 120 124 125 125 125 125
Phase(percent cycle) 0.82 0.70 0.58 0.46 0.34 0.26 0.22 0.22 0.22 0.22 0.22 0.22
Theta cycle 1 2 3 4 5 6 7 8 9 10 11 12
Cell 2
spiketime(ms) 227 340 451 562 675 794 917 1041 1166 1291 1416
ISI(ms) 113 111 111 113 119 123 124 125 125 125
Phase(percent cycle) 0.82 0.72 0.61 0.50 0.40 0.35 0.34 0.33 0.33 0.33 0.33
Theta cycle 2 3 4 5 6 7 8 9 10 11 12
Cell 3
spiketime(ms) 357 469 581 693 810 931 1055 1180 1305 1430
ISI(ms) 112 112 112 117 121 124 125 125 125
Phase(percent cycle) 0.86 0.75 0.65 0.54 0.48 0.45 0.44 0.44 0.44 0.44
Theta cycle 3 4 5 6 7 8 9 10 11 12
Cell 4
spiketime(ms) 607 717 829 946 1069 1194 1319 1444
ISI(ms) 110 112 117 123 125 125 125
Phase(percent cycle) 0.86 0.74 0.63 0.57 0.55 0.55 0.55 0.55
Theta cycle 5 6 7 8 9 10 11 12
- initialize
- load currents (own file)
- pyramidal cells and gammaNeuron to v=-60 and spiketime=-Inf
- thetaSpikes = [ 0.*[1x125]; 125 .*[1x125]; ...;4.*125/*[1x125] ]
- setinputs to pyrimdal cells ( 1,100 2,225 3,355 4,605)
- for each step of stimulation
- for each pyramidal cell
- update this pyramidal cell voltage
- update gamma cell voltage
- for each pyramidal cell
- collect stats
- plot gamma and pyramidal cells
- if cell spiked within last millisecond
- v =0
- return %don't do anything else
- for j as each current (leak, theta, ahp, adp, input, gamma)
- get "t"
- if current is theta use theta spikes
- if current is gamma use gammaNueron spike time
- otherwise use this pyramidal cell's spike time
- calculate "g"
- using above t and currents(j).tau_fall ...anrom, G, Erev
- add to numerator num=num+g.*(E-v) and denominator denom=denom+num+g portions of deltaV sum
- get "t"
- calculate detaV using collected sums=> num*dt /(1+dt *denom)
- if voltage is above thres AND timeline(i) > spiketime+2ms
- voltage = 0
- update volt of pyrmidal(p)
- update Vhist(p,i) = voltage