We present benchmark calculations of the Anderson lattice model based on the recently-developed ghost Gutzwiller approximation. Our analysis shows that, in some parameters regimes, the predictions of the standard Gutzwiller approximation can be incorrect by orders of magnitude for this model. We show that this is caused by the inability of this method to describe simultaneously the Mott physics and the hybridization between correlated and itinerant degrees of freedom (whose interplay often governs the metal-insulator transition in real materials). Finally, we show that the ghost Gutzwiller approximation solves this problem, providing us with results in remarkable agreement with dynamical mean field theory throughout the entire phase diagram, while being much less computationally demanding. We provide an analytical explanation of these findings and discuss their implications within the context of ab-initio computation of strongly-correlated matter.