code: def uncertainty_K(self):
          
          if self.information_ is None:
              self._build_information()
  
          sigma_K = _ratematrix.sigma_K(
              self.information_, theta=self.theta_, n=self.n_states_)
          return sigma_K
docs: Estimate of the element-wise asymptotic standard deviation
          in the rate matrix