code: def compute_theta(self):
          
          if self.tau == 1:
              theta = 10000
  
          else:
              theta = 2 * self.tau / (1 - self.tau)
  
          return theta
docs: Compute theta parameter using Kendall's tau.
  
          On Clayton copula this is :math:`τ = θ/(θ + 2) \\implies θ = 2τ/(1-τ)` with
          :math:`θ ∈ (0, ∞)`.
  
          On the corner case of :math:`τ = 1`, a big enough number is returned instead of infinity.