##### Plot Probability Density of an Histogram
 Plot Probability Density of an Histogram Python_TLS Unladen Swallow Posts: 2 Threads: 2 Joined: Jun 2018 Reputation: 0 Jun-27-2018, 11:37 AM Hi, Please =i have a code wich plot a histogram for each iteration and i want to plot the Probability Density of each histogram how can i do this please ? you find attached the code ```import matplotlib.pyplot as plt import numpy as np import pylab importdata=open('Hmz_OT_MULC_Discrete5_TOPO.dens','r+') data=importdata.readlines() importdata.close() tri=[] j=0 for i in range(len(data)): if 'iter' in data[i]: l=i nbre_iter=[] while l+1 < len(data) and 'iter' not in data [l+1]: l+=1 nbre_iter.append(data[l]) tri.append(nbre_iter) j+=1 tri_final=[] for i in range(len(tri)): extract_density=[] for j in range(len(tri[0])): extract_density.append(tri[i][j][12 :]) tri_final.append(extract_density) convert_into_float=[] for i in range(len(tri_final)): convert=[] for j in range(len(tri_final[0])): convert.append(float(tri_final[i][j])) convert_into_float.append(convert) for i in range(len(convert_into_float)): nbre_elements=np.zeros(10) for j in range(len(convert_into_float[0])): if 0<= convert_into_float[i][j]< 0.1: nbre_elements[0]+=1 elif 0.1<= convert_into_float[i][j] <0.2: nbre_elements[1]+=1 elif 0.2<= convert_into_float[i][j] <0.3: nbre_elements[2]+=1 elif 0.3<= convert_into_float[i][j] <0.4: nbre_elements[3]+=1 elif 0.4<= convert_into_float[i][j] <0.5: nbre_elements[4]+=1 elif 0.5<= convert_into_float[i][j] <0.6: nbre_elements[5]+=1 elif 0.6<= convert_into_float[i][j] <0.7: nbre_elements[6]+=1 elif 0.7<= convert_into_float[i][j] <0.8: nbre_elements[7]+=1 elif 0.8<= convert_into_float[i][j] <0.9: nbre_elements[8]+=1 elif 0.9<= convert_into_float[i][j] <=1: nbre_elements[9]+=1 fig, ax = plt.subplots() x=[1,2,3,4,5,6,7,8,9,10] width = 1.0 BarName = ['0.0-0.1','0.1-0.2','0.2-0.3','0.3-0.4','0.4-0.5','0.5-0.6','0.6-0.7','0.7-0.8','0.8-0.9','0.9-1'] rects=plt.bar(x, nbre_elements, width, color='b' ) def autolabel(rects): # attach some text labels for rect in rects: height = rect.get_height() ax.text(rect.get_x()+rect.get_width()/2., 1.0*height, '%d'%int(height), ha='center', va='bottom') plt.xlabel('Densite') plt.ylabel('Nombre d elements') plt.title('NNombre d elements en fonction de leur densite_ITER_'+str(i)) nom_figure='Histogramme'+str(i)+ '.png' pylab.xticks(x, BarName, rotation=40) autolabel(rects) plt.savefig(nom_figure) plt.show() ``` Reply Posts: 818 Threads: 1 Joined: Mar 2018 Reputation: 111 Jun-28-2018, 06:19 AM To plot probability density distribution from empirical data you need to estimate probability density first. One of the ways to do it is to use the kernel density estimation approach (see scipy's kde density estimator). Also, you can use ready-made function from `seaborn` package. ```import seaborn as sns import numpy as np import matplotlib.pyplot as plt x=np.random.randn(10000) sns.distplot(x) #plots histogram and kde-estimation of the pdf. plt.show()``` Reply

 Possibly Related Threads… Thread Author Replies Views Last Post inserting something into probability density SchroedingersLion 1 518 Jan-06-2020, 09:15 AM Last Post: Gribouillis Get underlying function from Kernel Density Estimation jpython 3 635 Dec-05-2019, 11:23 AM Last Post: jpython How to get the probability density function of my data set jpython 1 487 Dec-04-2019, 12:49 PM Last Post: Larz60+ finding the integral of probability density function Staph 3 729 Aug-11-2019, 09:19 AM Last Post: buran How to fit a 2D histogram like in ROOT? Zandar 2 1,123 Aug-06-2019, 03:01 PM Last Post: Zandar MatPlotLib 2d plot of current density? ruben 0 570 May-13-2019, 06:47 AM Last Post: ruben How to check for nested dataframe density? python_newbie09 0 735 Aug-27-2018, 07:34 PM Last Post: python_newbie09 Draw Weibull distribution probability function based on Confidence interval farzadtb 1 1,516 Jul-31-2018, 03:21 PM Last Post: Vysero histogram error alyssaantarctica 1 1,138 Jul-09-2018, 10:46 PM Last Post: Larz60+ histogram with matplotlib vaugirard 10 2,348 Jun-01-2018, 04:13 AM Last Post: vaugirard

Forum Jump:

### User Panel Messages

##### Announcements
Announcement #1 8/1/2020
Announcement #2 8/2/2020
Announcement #3 8/6/2020