Я собираюсь посмотреть поведение ансамблей данных.Лучший способ визуализировать ансамбль - использовать различные интервалы неопределенности и среднее значение.Мы можем вычислить интервал неопределенности в различных процентилях, используя "st.scoreatpercentile". Но, используя следующий код, я сталкиваюсь с ошибкой, и рисунок не отображается.
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as st
# generate some data
x = 100*np.sin([np.linspace(0,10,100)])
X = np.vstack([x,x,x,x,x,x,x,x,x,x,x,x,x,x,x,x,x,x,x,x])
e = 10*np.random.randn(20,100)
X_err = X+e
ll = st.scoreatpercentile(X_err, 10) # 10th percentile
ml = st.scoreatpercentile(X_err, 50) # 50th percentile
ul = st.scoreatpercentile(X_err, 90) # 90th percentile
plt.plot(ml,'g', lw=2, label= ' Median ' )
plt.plot(ul,'r', label= ' 90% ' )
plt.plot(ll,'b', label= ' 10% ' )
plt.xlabel( ' Time ' )
plt.ylabel( ' X ' )
plt.legend(loc= 'best' )
plt.show()