  Complex organic molecules (COMs) have been observed towards several low-mass
young stellar objects (LYSOs). Small and heterogeneous samples have so far
precluded conclusions on typical COM abundances, as well as the origin(s) of
abundance variations between sources. We present observations towards 16 deeply
embedded (Class 0/I) low-mass protostars using the IRAM 30m telescope. We
detect CH$_2$CO, CH$_3$CHO, CH$_3$OCH$_3$, CH$_3$OCHO, CH$_3$CN, HNCO, and
HC$_3$N towards 67%, 37%, 13%, 13%, 44%, 81%, and 75% of sources respectively.
Median column densities derived using survival analysis range between
6.0x10$^{10}$ cm$^{-2}$ (CH$_3$CN) and 2.4x10$^{12}$ cm$^{-2}$ (CH$_3$OCH$_3$)
and median abundances range between 0.48% (CH$_3$CN) and 16% (HNCO) with
respect to CH$_3$OH. Column densities for each molecule vary by about one order
of magnitude across the sample. Abundances with respect to CH$_3$OH are more
narrowly distributed, especially for oxygen-bearing species. We compare
observed median abundances with a chemical model for low-mass protostars and
find fair agreement, although some modeling work remains to bring abundances
higher with respect to CH$_3$OH. Median abundances with respect to CH$_3$OH in
LYSOs are also found to be generally comparable to observed abundances in hot
cores, hot corinos, and massive young stellar objects. Compared with comets,
our sample is comparable for all molecules except HC$_3$N and CH$_2$CO, which
likely become depleted at later evolutionary stages.
