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