www.neuron.yale.edu

Hi there, After a great tutorial of NetPyNE done by Salvador Dura I wanted to extend my question from the chat :) I would like to make a simulation where I will stimulate a network in a closed-loop with external stimulation similar to Reinforcement Learning paradigm. So let's say I have: The network...

Sorry Ted, I tried to prepare the code for you and make sure that there is no interference from the core I wrote, but then I switched to a different task, and left it for later, but I will prepare the code next week and let you know.

The function h.parent_connection() has wrong description and behavior: child.connect(parent(0.7), 0.0) h.parent_connection(child) # returns 0.0 child.parentseg().x # returns 0.7 Reference description of h.parent_connection() indeed points to the y value, however from the description you may expect t...

Thanks, that's an interesting comment. The problem arises with the Python object management and garbage collector. With simple objects, the removal is almost on demand, that's why you may see this scenario works well: class Experiment: def __init__(self): self.cell = Cell(name="cell") self...

Thanks, I will check the thread you mentioned. According to your own response from the last year (the first response to this thread): https://www.neuron.yale.edu/phpBB/viewtopic.php?f=2&t=1484&p=17633&hilit=delete_section#p17633 delete_section doesn't work in regular Python usage: You mu...

It doesn't work for me. I created a custom object-oriented framework to manage NEURON's Python wrapper. So I have objects such as Cell, Simulation or Record. In a simplistic view, what I done looks like this: class Environment: def run(self): cell = Cell(name="cell") cell.add_sec(name=&quo...

Hi there, I'm using Python wrapper. Before each run I would like to setup a fresh environment, so I want to remove all objects (sections, iclamps, netstims, point processes). Is there any way to do that? As indicated in the reference the function h.nrninit() is not helpful, so I see 3 ways: 1. Run e...

Dear Ted, I checked once more the code and I discovered that I use cvode(True), which means that there is a variable timestep. My current readouts for the same model but different compartment are: constant timestep with shape: ms: 50.00000000000341 mV: -8.280244904579009 no shape: ms: 50.00000000000...

Hi there! I use a model of Combe et al. 2018: https://doi.org/10.1523/JNEUROSCI.0449-18.2018 https://senselab.med.yale.edu/modeldb/ShowModel?model=244416 When I removed all tmpo2 variable usage from cell_setup.hoc: link to the file: https://senselab.med.yale.edu/modeldb/ShowModel?model=244416&fi...

Thanks for the fast and detailed answer! However I work with Hay et al. 2011 model. In one file, they call distribute_channels() function: https://senselab.med.yale.edu/modeldb/showModel.cshtml?model=180373&file=/ShaiEtAl2015/models/L5PCbiophys5.hoc#tabs-2 $o1.distribute_channels("apic"...

Hi! When I iterate over segments in a section using Python: [x for x in sec] I got the following values: [0.166667, 0.5, 0.833333] However if I iterate over the same section in HOC: for(x) {print x} I got: [0, 0.166667, 0.5, 0.833333, 1] So as you can see - the beginning and the end (0 and 1) are al...

Hi! I have a STDP-like POINT_PROCESS which I want to modulate with additional Dopamine stimulation (another POINT_PROCESS). I created a dummy da.mod POINT_PROCESS which only receives weight which is read by the STDP-like POINT_PROCESS (through POINTER). Then the weight on the dummy da.mod is zero-ou...

I want to perform multiple runs for the same simulation. Between runs I want to give new stimulations to my neurons, but do not reset the simulation - the network (and its state) must remain the same. The motivation for this comes from the artificial neural network field (ANN). Just like for ANN you...

Thank you for the fast and very detailed response! :) Currently I don't care so much about extracellular space so I wonder if there is any other way to pump out Ca2+? Because I assume that the extracellular computation maybe much more expensive. Also I would like to specify explicitly 2 reactions: c...

Hi! I would like to model PMCA pump on the cell membrane. In the Calcium wave example there is a serca pump model: cyt_er_membrane = rxd.Region(h.allsec(), geometry=rxd.DistributedBoundary(1)) serca = rxd.MultiCompartmentReaction(ca[cyt], ca[er], gserca / ((kserca / (1000. * ca[cyt])) ** 2 + 1), mem...