self.sigmay = np.sqrt(((self.rlevels + 1 - muy) ** 2 * py).sum(1))