Approach 1:
lm.fit <- glm(response ~ 1, offset=log(lam), family="quasipoisson")
summary(lm.fit)
Approach 2: Feed summary.glm with a pearson dispersion.
lm.fit <- glm(response ~ 1, offset=log(lam), family="poisson")
summary(lm.fit,
dispersion=sum(residuals(lm.fit, type="pearson")^2)/df.residual(lm.fit))