In code: instead of 'S' , I believe it should be S_obs
#Implement markov chain
for t in range(1,df.shape[0]):
alpha[t]=alpha_fn(DOY[t],c,phi)
S_pred[t]=S_min + alpha[t]*(S_obs[t-1]-S_min)+P[t]
S_pred[t]=np.min([S_pred[t],S_max])
I am also not getting the figures
I appreciate any help.
Thanks
Amitava Chatterjee
In code: instead of 'S' , I believe it should be S_obs
#Implement markov chain
for t in range(1,df.shape[0]):
alpha[t]=alpha_fn(DOY[t],c,phi)
S_pred[t]=S_min + alpha[t]*(S_obs[t-1]-S_min)+P[t]
S_pred[t]=np.min([S_pred[t],S_max])
I am also not getting the figures
I appreciate any help.
Thanks
Amitava Chatterjee