powered by
discrpareto2.mlekestimate(x, s, kmin=1e-04, kmax=100)
Y=floor(X)
ppareto2
discrpareto2.mleksestimate
discrpareto2.goftest