i am trying to get the the variance of derivative of order 2. i only have the data of an EEG signal sampled at 128 Hz. is the below code correct

Question

I basically want to implement equation 7 present in the added image

interpolation=CubicSpline(x,y)
x_interp = np.linspace(0, 1, 129)
y_interp = interpolation(x_interp,2)
var=np.var(y_interp)
df=((np.sqrt(222/72.7var))/(1+np.sqrt(222/72.7var)))

score 0 · Answer 1 · answered Apr 01 '24 at 16:24

It seems you've implemented Eq. 7, provided that var is indeed the fourth order moment.

However, your approximations for pi and e might not be precise enough. To make them more precise, use the numpy implementations.

L = np.sqrt(2*np.pi*np.exp(1)*var) / (1 + np.sqrt(2*np.pi*np.exp(1)*var))

i am trying to get the the variance of derivative of order 2. i only have the data of an EEG signal sampled at 128 Hz. is the below code correct

I basically want to implement equation 7 present in the added image

1 Answers1