Statistical Distributions of Areas of Voronoi Cells

Reference: https://community.wolfram.com/groups/-/m/t/1620190
In[]:=
pts=RandomPoint[Disk[],5000];​​mesh=VoronoiMesh[pts,AxesTrue]
In[]:=
vor=MeshPrimitives[mesh,2];​​vor//Length
In[]:=
vorInner=Select[vor,RegionWithin[Disk[],#]&];​​vorInner//Length
In[]:=
Graphics[{FaceForm[None],EdgeForm[Gray],vorInner}]
In[]:=
areas=Area/@vorInner;​​hist=Histogram[areas,Automatic,"PDF",PlotTheme"Detailed"]
In[]:=
dis=FindDistribution[areas]
In[]:=
Show[hist,Plot[PDF[dis,x],{x,0,.0015}]]
In[]:=
{Mean[dis],Variance[dis],Kurtosis[dis]}
In[]:=
Probability[x>.0005,xdis]