Let $Z$ be a 3$\times$2 matrix, $G$ be a 2$\times$2 positive definite symmetric (covariance) matrix, $\sigma^2$ a positive scalar, and $I$ a 2$\times$2 unit diagonal matrix.
Numerically I always see that
$$Z^{\prime}(ZGZ^{\prime}+\sigma^2I)^{-1}Z < (G^{-1}+Z^{\prime}Z/\sigma^2)$$
Is that true always ?