code: def uncertainty_pi(self):
          
          if self.information_ is None:
              self._build_information()
  
          sigma_pi = _ratematrix.sigma_pi(
              self.information_, theta=self.theta_, n=self.n_states_)
          return sigma_pi
docs: Estimate of the element-wise asymptotic standard deviation
          in the stationary distribution.