code: def _do_lumping(self):
          
          model = LandmarkAgglomerative(linkage='ward',
                                        n_clusters=self.n_macrostates,
                                        metric=self.metric,
                                        n_landmarks=self.n_landmarks,
                                        landmark_strategy=self.landmark_strategy,
                                        random_state=self.random_state)
          model.fit([self.transmat_])
  
          if self.fit_only:
              microstate_mapping_ = model.landmark_labels_
  
          else:
              microstate_mapping_ = model.transform([self.transmat_])[0]
  
          self.microstate_mapping_ = microstate_mapping_
docs: Do the MVCA lumping.