'Fit a model using MICE.
Parameters
n_burnin : int
The number of burn-in cycles to skip.
n_imputations : int
The number of data sets to impute'

def fit(self, n_burnin=10, n_imputations=10):

    self.data.update_all(n_burnin)
    for j in range(n_imputations):
        result = self.next_sample()
        self.results_list.append(result)
    self.endog_names = result.model.endog_names
    self.exog_names = result.model.exog_names
    return self.combine()
