The electronic ground states for two- and three-dimensional Anderson-lattice clusters with large Hubbard repulsion are determined variationally within the Gutzwiller ansatz by use of a stochastic algorithm for exact evaluation of matrix elements of Gutzwiller states. The algorithm is similar to one normally used for finite-temperature simulation and is described in some detail. Variational states for a variety of Hamiltonian parameters and particle densities were considered. Of particular interest is the regime of small interband hopping, which is antiferromagnetic at densities near two particles per unit cell. We found that for some Hamiltonian-parameter regimes, doping the antiferromagnetic cluster with enough particles or holes can destroy the magnetic ordering and produce a state that resembles a heavy-Fermi-liquid state, a metal with a strongly-frequency-dependent self-energy. For a 4\ifmmode\times\else\texttimes\fi{}4 Anderson lattice, we present results of low-temperature quantum Monte Carlo simulations to evaluate the validity of the variational Gutzwiller ansatz.