In this paper, we consider gradient estimates for positive solutions to the following equation $$\triangle_V u+au^p\log u=0$$ on complete noncompact Riemannian manifold with $k$-dimensional Bakry-Emery Ricci curvature bounded from below. Using the Bochner formula and the Cauchy inequality, we obtain upper bounds of $ \nabla u $ with respect to the lower bound of the Bakry-Emery Ricci curvature.