  We calculate the masses of the $0^{++}$, $0^{--}$ and $1^{+-}$ glueballs from
QCD in 3+1 dimensions using an eigenvalue equation method for Hamiltonian
lattice QCD developed and described elsewhere by the authors. The mass ratios
become approximately constants in the coupling region $6/g^2 \in [6.0,6.4]$,
from which we estimate $M(0^{--})/M(0^{++})=2.44 \pm 0.05 \pm 0.20$ and
$M(1^{+-})/M(0^{++})=1.91 \pm 0.05 \pm 0.12$.
