Low Frequency Stimulation model
Posted: Wed Oct 27, 2021 10:36 am
Hello NEURON community,
I am still very new to NEURON (I have only ever created two models with it), so I am very open to criticism and being pointed in any direction with a project I am working on.
I am attempting to replicate a study done in the link provided below:
https://pubmed.ncbi.nlm.nih.gov/23981713/
In short, I am opting to replicate the suppression of of epileptic activity through extracellular Low Frequency Stimulation (LFS). My hypothesis was that by utilizing an existing model (like shown below), I can use extracellular electrodes that output LFS to inhibit the epileptic activity in the program.
https://senselab.med.yale.edu/ModelDB/S ... 20/#tabs-1
My question pertains mostly with how to make extracellular AC electrodes. In the past I simple made a vectors in MATLAB with time and DC values that I input into NEURON with the following code:
objref pfile, pvector
pfile=new File("potentials.txt") //Voltage vector generated in MATLAB
pfile.ropen()
pvector=new Vector(51)
pvector.scanf(pfile)
pfile.close()
objref Ifile, Ivec
Ifile=new File("squarepulse.txt") //time vector generated in MATLAB
Ifile.ropen()
Ivec=new Vector()
Ivec.scanf(Ifile)
Ifile.close()
objref stimvec[51]
fract_dist_of_first_segment=(axon.L/nseg/2)/axon.L
inter_segment_fract_length=1/nsegobjref pfile, pvector
for i=0,50{
stimvec=Ivec.c.mul(pvector.x)
stimvec.play(&axon.e_extracellular(fract_dist_of_first_segment+i*inter_segment_fract_length),dt)}
run()
Would this still be practical using an existing model?
I am still very new to NEURON (I have only ever created two models with it), so I am very open to criticism and being pointed in any direction with a project I am working on.
I am attempting to replicate a study done in the link provided below:
https://pubmed.ncbi.nlm.nih.gov/23981713/
In short, I am opting to replicate the suppression of of epileptic activity through extracellular Low Frequency Stimulation (LFS). My hypothesis was that by utilizing an existing model (like shown below), I can use extracellular electrodes that output LFS to inhibit the epileptic activity in the program.
https://senselab.med.yale.edu/ModelDB/S ... 20/#tabs-1
My question pertains mostly with how to make extracellular AC electrodes. In the past I simple made a vectors in MATLAB with time and DC values that I input into NEURON with the following code:
objref pfile, pvector
pfile=new File("potentials.txt") //Voltage vector generated in MATLAB
pfile.ropen()
pvector=new Vector(51)
pvector.scanf(pfile)
pfile.close()
objref Ifile, Ivec
Ifile=new File("squarepulse.txt") //time vector generated in MATLAB
Ifile.ropen()
Ivec=new Vector()
Ivec.scanf(Ifile)
Ifile.close()
objref stimvec[51]
fract_dist_of_first_segment=(axon.L/nseg/2)/axon.L
inter_segment_fract_length=1/nsegobjref pfile, pvector
for i=0,50{
stimvec=Ivec.c.mul(pvector.x)
stimvec.play(&axon.e_extracellular(fract_dist_of_first_segment+i*inter_segment_fract_length),dt)}
run()
Would this still be practical using an existing model?