phat = distfn.fit2(rvs, **opt)