Variational Monte Carlo calculations for Bose-Hubbard model based on projected wavefunctions