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 I0 as the median of
the distribution (a more
difficult fit than 2 parameters)
