We study $L^p$-$L^q$ estimate for the spectral projection operator $��_��$ associated to the Hermite operator $H=|x|^2-��$ in $\mathbb R^d$. Here $��_��$ denotes the projection to the subspace spanned by the Hermite functions which are the eigenfunctions of $H$ with eigenvalue $��$. Such estimates were previously available only for $q=p'$, equivalently with $p=2$ or $q=2$ (by $TT^*$ argument) except for the estimates which are straightforward consequences of interpolation between those estimates. As shown in the works of Karadzhov, Thangavelu, and Koch and Tataru, the local and global estimates for $��_��$ are of different nature. Especially, $��_��$ exhibits complicated behaviors near the set $\sqrt��\mathbb S^{d-1}$. Compared with the spectral projection operator associated to the Laplacian, $L^p$-$L^q$ estimate for $��_��$ is not so well understood up to now for general $p,q$. In this paper we consider $L^p$--$L^q$ estimate for $��_��$ in a general framework including the local and global estimates with $1\le p\le 2\le q\le \infty$ and undertake the work of characterizing the sharp bounds on $��_��$. We establish various new sharp estimates in extended ranges of $p,q$. First of all, we provide a complete characterization of the local estimate for $��_��$ which was first considered by Thangavelu. Secondly, for $d\ge5$, we prove the endpoint $L^2$--$L^{2(d+3)/(d+1)}$ estimate for $��_��$ which has been left open since the work of Koch and Tataru. Thirdly, we extend the range of $p,q$ for which the operator $��_��$ is uniformly bounded from $L^p$ to $L^q$.

90 pages, 7 figures. The last section of the previous version was taken out to be published separately