code: def posteriori_covariance(self):
          r
          K = GLMM.covariance(self)
          tau = self._ep._posterior.tau
          return pinv(pinv(K) + diag(1 / tau))
docs: r""" Covariance of the estimated posteriori.