How the program currently works is that there are 0's in the clustered data. The uncertainty function takes in that clustered data, even though the uncertainty will be skewed due to the 0's. Maybe what needs to be fixed is you do something like
for data in clustered_data[:,~np.all(clustered_data, axis=3)]:
# calculate the covariance of data
How the program currently works is that there are 0's in the clustered data. The uncertainty function takes in that clustered data, even though the uncertainty will be skewed due to the 0's. Maybe what needs to be fixed is you do something like
when taking the covariance of all the data.