  Monte Carlo simulations alone could not clarify the corrections to scaling
for the size-dependent p_c(L) above the upper critical dimension. Including the
previous series estimate for the bulk threshold $p_c(\infty)$ gives preference
for the complicated corrections predicted by renormalization group and against
the simple 1/L extrapolation. Additional Monte-Carlo simulations using the
Leath method corroborate the series result for p_c.
