code: def uncertainty_timescales(self):
          
          if self.information_ is None:
              self._build_information()
  
          sigma_timescales = _ratematrix.sigma_timescales(
              self.information_, theta=self.theta_, n=self.n_states_)
  
          if self.n_timescales is None:
              return sigma_timescales
          return sigma_timescales[:self.n_timescales]
docs: Estimate of the element-wise asymptotic standard deviation
          in the model relaxation timescales.