We report the bivariate
- and
-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for
,
, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the
and
conditional distributions and the galaxy stellar mass function (GSMF). For the conditional distributions we use the results from the compilation presented in Calette et al. [(2018) RMxAA, 54, 443.]. For determining the GSMF from
to
, we combine two spectroscopic samples from the Sloan Digital Sky Survey at the redshift range
. We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is
, consistent with previous determinations. The obtained
MFs agree with radio blind surveys. Similarly, the
MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while
of the
and
mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of
and total neutral H in late-type galaxies of
and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting
into stars and in transforming
gas into
. Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas–stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.