Search found 13 matches

by erhervert
Mon Feb 27, 2012 11:29 pm
Forum: Other questions
Topic: Data analysis
Replies: 8
Views: 3059

Re: Data analysis

I found it addictive, feels great when the things start working good.
by erhervert
Mon Feb 27, 2012 2:29 pm
Forum: Other questions
Topic: Data analysis
Replies: 8
Views: 3059

Re: Data analysis

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") savdata.close() } proc doal...
by erhervert
Sun Feb 26, 2012 5:56 pm
Forum: Other questions
Topic: Data analysis
Replies: 8
Views: 3059

Re: Data analysis

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...
by erhervert
Sat Feb 25, 2012 12:46 pm
Forum: Other questions
Topic: Question about proc batchrun()
Replies: 6
Views: 2868

Re: Question about proc batchrun()

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 (...
by erhervert
Sat Feb 25, 2012 11:51 am
Forum: Other questions
Topic: Data analysis
Replies: 8
Views: 3059

Re: Data analysis

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(...
by erhervert
Sat Feb 25, 2012 9:19 am
Forum: Other questions
Topic: Question about proc batchrun()
Replies: 6
Views: 2868

Re: Question about proc batchrun()

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...
by erhervert
Fri Feb 24, 2012 9:17 pm
Forum: Other questions
Topic: Data analysis
Replies: 8
Views: 3059

Data analysis

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 one...
by erhervert
Fri Feb 24, 2012 5:36 pm
Forum: Other questions
Topic: Question about proc batchrun()
Replies: 6
Views: 2868

Re: Question about proc batchrun()

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 ...
by erhervert
Thu Feb 23, 2012 2:02 pm
Forum: Other questions
Topic: Question about proc batchrun()
Replies: 6
Views: 2868

Question about proc batchrun()

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 with t...
by erhervert
Fri Feb 03, 2012 11:38 pm
Forum: Optimization
Topic: automating exploration of parameter space
Replies: 1
Views: 5822

automating exploration of parameter space

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 in each i...
by erhervert
Thu Feb 02, 2012 1:23 am
Forum: Adding new mechanisms and functions to NEURON
Topic: Modelling conductances with H&H model for ionic currents
Replies: 5
Views: 14197

Re: Modelling conductances with H&H model for ionic currents

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...
by erhervert
Tue Jan 31, 2012 7:01 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Modelling conductances with H&H model for ionic currents
Replies: 5
Views: 14197

Re: Modelling conductances with H&H model for ionic currents

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...
by erhervert
Mon Jan 30, 2012 7:47 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Modelling conductances with H&H model for ionic currents
Replies: 5
Views: 14197

Modelling conductances with H&H model for ionic currents

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...