I'm a Msc student,
I am trying to implement the extracellular stimulation and recording on a model for retina network. For attempt one, I consider having just one ganglion cell. I took the code from the example provided by Ted for extracellular stimulating (http://www.neuron.yale.edu/ftp/ted/neur ... nd_rec.zip). I created the files parameters.hoc and createcells.hoc instead of cell.ses, which was created by saving a CellBuilder to a ses file. The files parameters.hoc and createcells.hoc are for creating one ganglion cell. So my initxstim.hoc file is as following:
Code: Select all
load_file("nrngui.hoc")
load_file("parameters.hoc")
load_file("createcells.hoc")
load_file("anatscale.hoc") // show xyz scale bars
load_file("interpxyz.hoc") // only interpolates sections that have extracellular
load_file("setpointers.hoc") // automatically calls grindaway() in interpxyz.hoc
load_file("field.hoc") // computes extracellularly recorded potential vrec
load_file("vrecc.ses") // graph of vrec(t)
load_file("calcrxc.hoc") // computes transfer r between segments and recording electrodes
load_file("stim.hoc") // extracellular stimulus
Code: Select all
begintemplate Ganglion
public soma,spiketimes, spikecount
create soma
objref spikecount, spiketimes,vc
proc init() {
create soma
spiketimes = new Vector()
soma {
L=25
diam=25
insert pas
e_pas=-65
g_pas=.00005
Ra=110
insert spike
gkbar_spike = 0.012
gabar_spike = 0.036
gcabar_spike = 0.0022
gkcbar_spike = 0.00005
celsius = 22
ena=35
ek=-75
spikecount = new APCount(0.5)
spikecount.thresh = -30
spikecount.record(spiketimes)
vc = new VClamp(0.5)
vc.dur[0]=100
vc.amp[0]=-65
}
forsec "soma" {
insert extracellular xraxial=1e9 xg=1e9 xc=0 e_extracellular=0
insert xtra
}
forall {
define_shape()
}
} // end init()
endtemplate Ganglion
I read the posts about extracellular stimulating in forum but didn’t get any idea.
Can anyone in here help me, please?
Thanks for every suggestion.
Best Regards