Using an improved lattice Hamiltonian with massive Wilson quarks a variational method is applied to study the dependence of the vector meson mass Mv on the quark mass m and the Wilson parameter r in in the scaling window 1 ≤ 1/g2 ≤ 2, Mv/g is approximately linear in m, but Mv/g obviously does not depend on r (this differs from the quark condensate). Particularly for m → 0 our numerical results agree very well with Bhattacharya’s analytical strong coupling result in the continuum, and the value of ((e)Mv/(e)m) |mm=0 in two-dimensional SU(NC) lattice gauge theory is very close to that in Schwinger model.