code: def _mse_converged(self):
          
  
          prior = self.global_prior_[0:self.prior_size]
          posterior = self.global_posterior_[0:self.prior_size]
          mse = mean_squared_error(prior, posterior,
                                   multioutput='uniform_average')
          if mse > self.threshold:
              return False, mse
          else:
              return True, mse
docs: Check convergence based on mean squared difference between
              prior and posterior
  
          Returns
          -------
  
          converged : boolean
              Whether the parameter estimation converged.
  
          mse : float
              Mean squared error between prior and posterior.