code: def get_anchor(self):
          
          temp = np.empty([self.n_nodes, 2])
          temp[:, 0] = np.arange(self.n_nodes, dtype=int)
          temp[:, 1] = np.sum(abs(self.tau_matrix), 1)
          anchor = int(temp[0, 0])
          return anchor
docs: Find anchor variable with highest sum of dependence with the rest.