We give estimators with bias exponentially small as the sample size increases for a broad class of multivariate power series families for lattice random variables. Examples include the Poisson, multinomial and negative multinomial distributions. Their use is illustrated by simulations and a real data application.