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)