Update eda_intervalrelated.py

Compute the mean peaks amplitude only on the peaks (and not on the whole signal)
neuropsychology · Oct 2, 2023 · b05a700 · b05a700
1 parent ce33299
commit b05a700
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/neurokit2/eda/eda_intervalrelated.py b/neurokit2/eda/eda_intervalrelated.py
@@ -114,7 +114,7 @@ def _eda_intervalrelated(
         )
         output["SCR_Peaks_Amplitude_Mean"] = np.nan
     else:
-        output["SCR_Peaks_Amplitude_Mean"] = np.nanmean(data["SCR_Amplitude"].values)
+        output["SCR_Peaks_Amplitude_Mean"] = np.nanmean(data[data["SCR_Peaks"] == 1]["SCR_Amplitude"].values)
 
     # Get variability of tonic
     if "EDA_Tonic" in colnames: