bluecellulab 2.6.38__tar.gz → 2.6.39__tar.gz

{bluecellulab-2.6.38 → bluecellulab-2.6.39}/.gitignore

@@ -1,5 +1,6 @@
 /modlib
 **/x86_64
+**/arm64
 /.tox
 *.pyc
 *.egg-info
@@ -26,3 +27,4 @@ nrnivmodl.log
 *.btr
 *.whl
 .ipynb_checkpoints
+cADpyr_L2TPC_dendrogram.png

{bluecellulab-2.6.38 → bluecellulab-2.6.39}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: bluecellulab
-Version: 2.6.38
+Version: 2.6.39
 Summary: Biologically detailed neural network simulations and analysis.
 Author: Blue Brain Project, EPFL
 License: Apache2.0

bluecellulab-2.6.39/bluecellulab/analysis/analysis.py

@@ -0,0 +1,63 @@
+"""Module for analyzing cell simulation results."""
+try:
+    import efel
+except ImportError:
+    efel = None
+import numpy as np
+
+from bluecellulab.stimulus import StimulusFactory
+from bluecellulab.tools import calculate_rheobase
+from bluecellulab.analysis.inject_sequence import run_stimulus
+from bluecellulab.analysis.plotting import plot_iv_curve
+
+
+def compute_plot_iv_curve(cell, stim_start=100.0, duration=500.0, post_delay=100.0, threshold_voltage=-30, nb_bins=11):
+    """Compute and plot the IV curve from a given cell by injecting a
+    predefined range of currents.
+
+    Args:
+        cell (bluecellulab.cell.Cell): The initialized cell model.
+        stim_start (float): Start time for current injection (in ms). Default is 100.0 ms.
+        duration (float): Duration of current injection (in ms). Default is 500.0 ms.
+        post_delay (float): Delay after the stimulation ends (in ms). Default is 100.0 ms.
+        nb_bins (int): Number of current injection levels. Default is 11.
+
+    Returns:
+        tuple: A tuple containing:
+            - list_amp (np.ndarray): The predefined injected step current levels (nA).
+            - steady_states (np.ndarray): The corresponding steady-state voltages (mV).
+    """
+    rheobase = calculate_rheobase(cell)
+
+    list_amp = np.linspace(rheobase - 2, rheobase - 0.1, nb_bins)  # [nA]
+
+    steps = []
+    times = []
+    voltages = []
+    # inject step current and record voltage response
+    for amp in list_amp:
+        stim_factory = StimulusFactory(dt=0.1)
+        step_stimulus = stim_factory.step(pre_delay=stim_start, duration=duration, post_delay=post_delay, amplitude=amp)
+        recording = run_stimulus(cell.template_params, step_stimulus, section="soma[0]", segment=0.5)
+        steps.append(step_stimulus)
+        times.append(recording.time)
+        voltages.append(recording.voltage)
+
+    steady_states = []
+    # compute steady state response
+    efel.set_setting('Threshold', threshold_voltage)
+    for voltage, t in zip(voltages, times):
+        trace = {
+            'T': t,
+            'V': voltage,
+            'stim_start': [stim_start],
+            'stim_end': [stim_start + duration]
+        }
+        features_results = efel.get_feature_values([trace], ['steady_state_voltage_stimend'])
+        steady_state = features_results[0]['steady_state_voltage_stimend']
+        steady_states.append(steady_state)
+
+    # plot I-V curve
+    plot_iv_curve(list_amp, steady_states)
+
+    return np.array(list_amp), np.array(steady_states)

bluecellulab-2.6.39/bluecellulab/analysis/plotting.py

@@ -0,0 +1,25 @@
+"""Module for plotting analysis results of cell simulations."""
+
+import matplotlib.pyplot as plt
+
+
+def plot_iv_curve(currents, voltages):
+    """Plots the IV curve.
+
+    Args:
+        currents (iterable): The injected current levels (nA).
+        voltages (iterable): The corresponding steady-state voltages (mV).
+    Raises:
+        ValueError: If the lengths of currents and voltages do not match.
+    """
+    if len(currents) != len(voltages):
+        raise ValueError("currents and voltages must have the same length")
+
+    # Plotting
+    plt.figure(figsize=(10, 6))
+    plt.plot(voltages, currents, marker='o', linestyle='-', color='b')
+    plt.title("I-V Curve")
+    plt.ylabel("Injected current [nA]")
+    plt.xlabel("Steady state voltage [mV]")
+    plt.tight_layout()
+    plt.show()

{bluecellulab-2.6.38 → bluecellulab-2.6.39}/bluecellulab/tools.py

@@ -23,6 +23,7 @@ import neuron
 import numpy as np
 
 import bluecellulab
+from bluecellulab.cell import Cell
 from bluecellulab.circuit.circuit_access import EmodelProperties
 from bluecellulab.exceptions import UnsteadyCellError
 from bluecellulab.simulation.parallel import IsolatedProcess
@@ -351,3 +352,76 @@ def check_empty_topology() -> bool:
         neuron.h.topology()
 
     return stdout == ['', '']
+
+
+def calculate_max_thresh_current(cell: Cell, threshold_voltage: float = -30.0) -> float:
+    """Calculate the upper bound threshold current.
+
+    Args:
+        cell (bluecellulab.cell.Cell): The initialized cell model.
+        threshold_voltage (float, optional): Voltage threshold for spike detection. Default is -30.0 mV.
+
+    Returns:
+        float: The upper bound threshold current.
+    """
+    # Calculate resting membrane potential (rmp)
+    rmp = calculate_SS_voltage(
+        template_path=cell.template_params.template_filepath,
+        morphology_path=cell.template_params.morph_filepath,
+        template_format=cell.template_params.template_format,
+        emodel_properties=cell.template_params.emodel_properties,
+        step_level=0.0,
+    )
+
+    # Calculate input resistance (rin)
+    rin = calculate_input_resistance(
+        template_path=cell.template_params.template_filepath,
+        morphology_path=cell.template_params.morph_filepath,
+        template_format=cell.template_params.template_format,
+        emodel_properties=cell.template_params.emodel_properties,
+    )
+
+    # Calculate upperbound threshold current
+    upperbound_threshold_current = (threshold_voltage - rmp) / rin
+    upperbound_threshold_current = np.min([upperbound_threshold_current, 2.0])
+
+    return upperbound_threshold_current
+
+
+def calculate_rheobase(cell: Cell,
+                       threshold_voltage: float = -30.0,
+                       threshold_search_stim_start: float = 300.0,
+                       threshold_search_stim_stop: float = 1000.0) -> float:
+    """Calculate the rheobase by first computing the upper bound threshold
+    current.
+
+    Args:
+        cell (bluecellulab.cell.Cell): The initialized cell model.
+        threshold_voltage (float, optional): Voltage threshold for spike detection. Default is -30.0 mV.
+        threshold_search_stim_start (float, optional): Start time for threshold search stimulation (in ms). Default is 300.0 ms.
+        threshold_search_stim_stop (float, optional): Stop time for threshold search stimulation (in ms). Default is 1000.0 ms.
+
+    Returns:
+        float: The rheobase current.
+    """
+    if cell.template_params.emodel_properties is None:
+        raise ValueError("emodel_properties cannot be None")
+
+    # Calculate upper bound threshold current
+    upperbound_threshold_current = calculate_max_thresh_current(cell, threshold_voltage)
+
+    # Compute rheobase
+    rheobase = search_threshold_current(
+        template_name=cell.template_params.template_filepath,
+        morphology_path=cell.template_params.morph_filepath,
+        template_format=cell.template_params.template_format,
+        emodel_properties=cell.template_params.emodel_properties,
+        hyp_level=cell.template_params.emodel_properties.holding_current,
+        inj_start=threshold_search_stim_start,
+        inj_stop=threshold_search_stim_stop,
+        min_current=cell.template_params.emodel_properties.holding_current,
+        max_current=upperbound_threshold_current,
+        current_precision=0.005,
+    )
+
+    return rheobase

{bluecellulab-2.6.38 → bluecellulab-2.6.39}/bluecellulab.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: bluecellulab
-Version: 2.6.38
+Version: 2.6.39
 Summary: Biologically detailed neural network simulations and analysis.
 Author: Blue Brain Project, EPFL
 License: Apache2.0

{bluecellulab-2.6.38 → bluecellulab-2.6.39}/bluecellulab.egg-info/SOURCES.txt

@@ -39,7 +39,9 @@ bluecellulab.egg-info/dependency_links.txt
 bluecellulab.egg-info/requires.txt
 bluecellulab.egg-info/top_level.txt
 bluecellulab/analysis/__init__.py
+bluecellulab/analysis/analysis.py
 bluecellulab/analysis/inject_sequence.py
+bluecellulab/analysis/plotting.py
 bluecellulab/cell/__init__.py
 bluecellulab/cell/cell_dict.py
 bluecellulab/cell/core.py