phat = FitDistribution(dist, data, floc=0.5, method='ml')