code: def getEmpTraitCovar(self):
          
          if self.P==1:
              out=self.Y[self.Iok].var()
          else:
              out=SP.cov(self.Y[self.Iok].T)
          return out
docs: Returns the empirical trait covariance matrix