I have a question regarding handling large amounts of data.
I have a model cell. I am producing an Injected current - firing frequency (I-F) curve. For each value injected current, I would like the time series for membrane potential, so that I can do further analysis. I therefore have saved the results of my protocol to a matrix, and written it to a datafile.
Trouble is, I want to do this same protocol, for various values of a parameter, GKA.
I have some code that looks like this
Code: Select all
for i=1,nvalues {
GKA=value(i)
for j=1,nstimvalues {
stim.amp=ampvalue(i) // this is the injected current amplitude
// set up vectors recording V, etc.
run()
}
// code for saving data for that value of GKA to a file
}
I guess this is equivalent to wanting to save time-series data from the results of a 2dimensional parameter sweep; in GKA and the stimulus amplitude.
Any suggestions would be most appreciated,
Linford