diff --git a/neurokit2/eda/eda_intervalrelated.py b/neurokit2/eda/eda_intervalrelated.py
index 8897d0fc6e..506e8034da 100644
--- 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: