We show that sampling effects in the initial mass function are very important in the low-mass cluster case. To this aim, we compute photoionization models ionized by realistic clusters made up of various combinations of individual stars and clusters built with a synthesis model. We discuss the differences in the position on diagnostic diagrams, their implications, and future and ongoing applications of the present work.