Aims: We study how IMF sampling affects the ionizing flux and emission line spectra of low mass stellar clusters. Methods: We performed 2 × 106 Monte Carlo simulations of zero-age solar-metallicity stellar clusters covering the 20-106 M&sun; mass range. We study the distribution of cluster stellar masses, Mclus, ionizing fluxes, Q(H0), and effective temperatures, T_eff^clus. We compute photoionization models that broadly describe the results of the simulations and compare them with photoionization grids. Results: Our main results are: (a) a large number of low mass clusters (80% for Mclus = 100 M&sun;) are unable to form an H ii region. (b) There are a few overluminous stellar clusters that form H ii regions. These overluminous clusters preserve statistically the mean value of obtained by synthesis models, but the mean value cannot be used as a description of particular clusters. (c) The ionizing continuum of clusters with Mclus ⪉ 104 M&sun; is more accurately described by an individual star with self-consistent effective temperature (T_eff^*) and Q(H0) than by the ensemble of stars (or a cluster T_eff^clus) produced by synthesis models. (d) Photoionization grids of stellar clusters cannot be used to derive the global properties of low mass clusters. Conclusions: Although variations in the upper mass limit, mup, of the IMF would reproduce the effects of IMF sampling, we find that an ad hoc law that relates mup to Mclus in the modeling of stellar clusters is useless, since: (a) it does not cover the whole range of possible cases; and (b) the modeling of stellar clusters with an IMF is motivated by the need to derive the global properties of the cluster: however, in clusters affected by sampling effects we have no access to global information of the cluster but only particular information about a few individual stars.