ecg_process(): Sequence too short to compute heart rate - optimising for short signal about neurokit.py

[DominiqueMakowski]

Hey @dicobrazz.

Short answer is, not really (to my knowledge) ^^. Within such time ranges, I expect you to have around 3 to 4 heart beats per epoch. Let's say 4 for the sake of illustration. Thus, you can only compute 3 different momentary heart rates (b1-b0, b2-b1...). Hence, while you could compute the average (of these several momentary hear rate values) over the whole window, computing the "continuous" heart rate is kind of complicated (and I believed not advised). Indeed, the rationale behind this artificial signal is that the momentary heart rate points (that differs at each new heartbeat) are the reflection of a smooth underlying process.

However, attempting at recovering this underlying process (i.e., the continuous and smooth change of heart rate) heavily rely on interpolation and filtering, which do not work well (and can lead to artefacts) for very short time periods.

Thus, depending on what you need this heart rate for, I would recommend 1) to use the average heart rate, or 2) use the momentary heart rates (made of 2 or 3 points). But again, it depends on what you want to do in the end.