'Simulate model parameters from fitted sampling distribution.'

def _simulate_params(self, result):

    mn = result.params
    cov = result.cov_params()
    return np.random.multivariate_normal(mn, cov)
