We present a Bayesian stacking technique to directly measure the H I mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several ...