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