NeuroKit icon indicating copy to clipboard operation
NeuroKit copied to clipboard

Event-related ECG/EDA Questions.

Open Suitieee opened this issue 2 years ago • 1 comments

Thanks a lot for this package, but I still got stuck by some problems when I set my time window and detect SCR.

Question and context:

First of all, with the same time window, 4s, I set the same codes. But I didn't figure out why the baseline is ranging from -100 to 0 on 'https://neurokit2.readthedocs.io/en/latest/examples/eventrelated.html', while an epoch is just ranging from -0.1s to 4s on the paper.

Secondly, a 4 seconds time window may fail to detect SCR_peaks, which results in no SCR_Amplitude being saved as a SCR_Magnitude of one epoch.

    # Feature 2 EDA - SCR
    # Sampling rate is 200Hz, and an epoch is ranging from -1s to 4s
    scr_max_4s = epoch["SCR_Amplitude"].loc[0:800].max()  # Maximum SCR peak
    # If no SCR, consider the magnitude, i.e.  that the value is 0
    if np.isnan(scr_max_4s):
        scr_max_4s = 0
    else:
        scr_num_64_4s = scr_num_64_4s + 1
    # Store SCR in df
    df64[epoch_index]["SCR_Magnitude_4s"] = scr_max_4s

After that, 4s is also not enough for the HRV detection, as we all know the frequency domain parameters, LF and HF, are ranging from 0.04Hz to 0.4Hz. So I'm wondering how to detect HRV in an event-related study, with such a short time window.

Suitieee avatar Mar 03 '22 06:03 Suitieee

Hi 👋 Thanks for reaching out and opening your first issue here! We'll try to come back to you as soon as possible. ❤️ kenobi

welcome[bot] avatar Mar 03 '22 06:03 welcome[bot]

Is this still relevant? If so, what is blocking it? Is there anything you can do to help move it forward?

This issue has been automatically marked as stale because it has not had recent activity. It will be closed if no further activity occurs.

stale[bot] avatar Nov 25 '22 04:11 stale[bot]