Open Florian-Soyka opened 3 years ago
I guess I can extract the time of the action potential from R-SENN. The place is (for a uniform E-field) always the first or the last node. So it is not so useful. But maybe for other exposure settings it might be. So I'll try to extract the node as well. Like you said, G-SENN can deal with more than one AP. We'd have to see how to include this in the GUI.
I added the AP time and node information to the results table in the GUI. That slightly changed the threshold value as well. Previously the last entry in the table was used (see picture in previous comment). Now the last entry for which the NXGT entry is equal to the number of consecutively activated nodes is used, resulting in slightly raised thresholds. I guess that is just a definition issue of which value to use from the titration procedure. Which one does G-SENN use?
For G-SENN we could take the first action potential and put it in the results table - and omit the others. One can still get the detailed results table for G-SENN by specifying an E-field value.
I just put our thoughts from the email here:
I thought it might be useful to also put the time at which the action potential is elicited into the results table.
I agree that also putting the time and location of the action potential in the results table would be very useful: this is indeed also possible with our SENN code (with ActivationTable=1 and SweepActI = 0, SENNfiberCP will calculate the location/position and also the speed with which the AP propagates in the orthodromic and antidromic directions). The activation table in G-SENN will list all locations/times at which action potentials are initiated (in most cases, this is only one location, but for some waveforms it is possible that action potentials start at multiple locations almost simultaneously). We can make an additional table in the GUI or make the existing table scrollable, to include the information of the activation table. For example, the AP location and latency can be calculated after the threshold is determined for a current injection just above the threshold (e.g., 10% higher than the threshold value).
I think R-SENN provides the time at which the maximum voltage is reached and the node at which this happens. I think with the current parameters this is always the first or the last node if I’m not mistaken. I’ll try to extract that information so that we can put it in the EONS results table.