I ran ICA on the same dataset using two different parameters:
(1) band-pass 3-55Hz, notch filter 60, 120 Hz, 25 ICA components: 2 heartbeats components, when removing them the data looks cleaner
(2) band-pass 0.3-150Hz, notch filter 60, 120Hz, 40 ICA components: 2 heartbeats components, when removing them the data looks more contaminated with heartbeats.
Do you think there is a bug in the code?