code: def getTraitCorrCoef(self,term_i=None):
          
          cov = self.getTraitCovar(term_i)
          stds = sp.sqrt(cov.diagonal())[:,sp.newaxis]
          RV = cov / stds / stds.T
          return RV
docs: Return the estimated trait correlation coefficient matrix for term_i (or the total if term_i is None)
              To retrieve the trait covariance matrix use \see getTraitCovar
  
          Args:
              term_i:     index of the random effect term we want to retrieve the correlation coefficients
          Returns:
              estimated trait correlation coefficient matrix