phats = FitDistribution(dist, data, floc=0.5, method='mps')