code: def fit(self, X):
          
          U, V = self.split_matrix(X)
          self.tau = stats.kendalltau(U, V)[0]
          self.theta = self.compute_theta()
          self.check_theta()
docs: Fit a model to the data updating the parameters.
  
          Args:
              X: `np.ndarray` of shape (,2).
  
          Return:
              None