Radio luminosity function at 3 GHz of our RS-NIR-faint galaxies. The points show the values computed as described in Section 4.1, while the dashed vertical lines indicate the minimum radio luminosity observable in a given redshift bin given the sensitivity of our survey. All the points corresponding to fainter luminosities (shaded points) are not included in the fitting. The shaded areas report the fitting of the modified Schechter function, with its 1 σ uncertainty. For reference, we also report the radio luminosity function (converted to 3 GHz by assuming a standard radio slope of α = −0.7) estimated by Novak et al. (2017), Enia et al. (2022), and van der Vlugt et al. (2022) on their samples of radio NIR-bright galaxies.