www.neuron.yale.edu

I found it addictive, feels great when the things start working good.

Done objref savdata, tempmatrix proc savetofile() { savdata = new File() savdata.wopen("results.dat") savdata.printf("t avg\n") tempmatrix = new Matrix() tempmatrix.resize(avg.size(),2) tempmatrix.setcol(0, time) tempmatrix.setcol(1, avg) tempmatrix.fprint(savdata, " %g"...

First a comment: if you don't prevent negative values from being assigned to your synapse's start time, you're going to have some runs that produce nonuseful results. Sorry, I did it already using the tsynon() function as you recommend but I posted the uncorrected code (by mistake), this is the goo...

You are right. I implement this in my program and it works fine. objref r r = new Random() {r.normal(.5, .05)} objref avg func tsynon() { ton = -1 while (ton<0) { ton = (((10-5)*r.repick())+5) } return ton } proc batchavg() { local i for i=0, $1 { syn.onset = tsynon() run() } avg = new Vector() if (...

Ok, I incorporate that into my program and this is what I got: /////////////////////////// //* simulation control *// /////////////////////////// dt = 0.025 tstop = 50 v_init = -65 objref r r = new Random() r.normal(.5, .05) objref avg proc batchavg() { local i, index for i=0, $1 { index = r.repick(...

Well, if I call batchrun(2) I only see one run, which corresponds to i == 1. Changing for i=1, $1-1 { to for i=1, $1 { makes batchrun(2) generate two runs, which seems to be an improvement. That's true. I also added a number generator to vary the value of syn.onset randomly. This is what I got: obj...

Hello again, I'd like to know if there is a way to customize the output ASCII file obtained from a simulation for further analysis. For example, if I execute my model n times and then I used the print & file window manager to select the V(t) plot, what I get is a file with two columns, one x and...

Hello Ted! To give an answer to the first question I'm working with a model published in the following paper --> Hines ML, Davison AP, Muller E. NEURON and Python. Front Neuroinform. 2009;3:1. Epub 2009 Jan 28 and I really don't know why they used a modified version of the standard run system but I ...

I'm working along with a cell model taken from Muller et al.,2009. I'm trying to execute multiple simulations varying the syn.onset parameter using the batch run procedure according with http://www.neuron.yale.edu/phpbb/viewtopic.php?f=28&t=2333 . The problem is that when I run the simulation wi...

I'd like to ask you a couple of things more regarding the control of simulations. I was wondering if it's possible to repeat a simulation "n" times varying parameters in each run, like for example the onset of the action potential or the stimulus intensity without have to change it by hand...

I think you are right Ted. In fact, what I'm trying to do is to reproduce the action potential of a sympathetic neuron and I was thinking in use that paper as a basis but I'm not constricted to it. Now, I'm considering to look for more recent papers with more detailed descriptions or source code for...

Hi Ted, Well, thanks for reply. The complete citation of the model that I'm trying to reproduce is the next: Belluzzi O., Sacchi O. (1991) A five-conductance model of the action potential in the rat sympathetic neurone. Prog Biophys Mol Biol. 1991;55(1):1-30. Review. I checked the paper again lookin...

Hi there, Well, first of all I'd like to say that this is my first approach to NMODL and I'm currently in the learning process... Now, the issue is that I'm trying to insert some conductances in a model but most of the available examples on the internet are wrote in terms of the alpha and beta funct...