phat = ws.norm.fit2(np.log(Hs))