Кривая гистограммы с биномиальным распределением - PullRequest
0 голосов
/ 26 апреля 2019

У меня есть гистограмма, которая представляет степень распределения сети для меня.Мне нужно приспособить его к биномиальному распределению, но поскольку в Scipy нет метода .fit для дискретных распределений, я не знаю, как получить параметры, необходимые для биномиальной функции.

Похоже, яНе получается получить правильные параметры из гистограммы, так как биномиальный график не соответствует форме гистограммы.что я делаю не так?

G = nx.gnm_random_graph(5,5)# My Graph

d=np.array(list(dict(nx.degree(G,weight='weight')).values()))#Calculating 
degree distribution
c = np.array(range(0,len(d),1))# The number of nodes

h = plt.hist(d,density= True ,bins=25, stacked = True)#histogram

bionom_func= st.binom.pmf(k,n,p)# defining the binomial function


popt, pcov = optimize.curve_fit(f=bionom_func, xdata=c, ydata=d)

estimated_n = popt[0]
estimated_p = popt[1]
modeled_y = st.binom.pmf(c, estimated_n, estimated_p)
plt.plot(c, modeled_y, 'r-')# plotting the fit curve


plt.show()
...