def rvs_single(self): sample = self.X.sample(weights=self.w).iloc[0] perturbed = (sample + np.random.multivariate_normal( np.zeros(self.cov.shape[0]), self.cov)) return perturbed