Making use of the following:
We can integrate
to get the cumulative distribution function
Which we can fit to the data by maximum-likelihood
methods. Here I fit
q
but fix
I
0
as the median of the
distribution (a more difficult fit than 2 parameters)