The shapes of interspecific range-size distributions at scales finer than the geographic range are highly variable. However, no numerical model has been developed as a basis for understanding this variation. Using self-similarity conditions, we present an occupancy probability transition (OPT) model to investigate the effect of sampling scale (i.e. sample grain) and species saturation (strongly positively correlated with the fractal dimension) on the shape of occupancy frequency distributions (fine scale expression of range-size distributions). In accordance with empirical observations, the model showed that core-modes are likely to be rare in occupancy frequency distributions. The modal occupancy shifted from core to satellite with an increase in sample grain (from coarse scale to fine scale) at a linear rate after log-transformation of occupancy. Saturation coefficients above a particular threshold generated multimodality. Bimodal distributions arose from a combination of different occupancy probability distributions (OPDs), with species-specific saturation coefficients generating occupancy frequency distributions of the shape commonly observed empirically, i.e. bimodal with a dominant satellite mode. This is a consequence of the statistical properties of the OPD, and is also largely insensitive to species richness. The OPT model thus provides a null model for the shape of occupancy frequency distributions. Furthermore, it demonstrates that the sample grain of a study, sampling adequacy (based on a linear sampling assumption) and the distribution of species saturation coefficients in a community are together largely able to explain the patterns observed in empirical occupancy frequency distributions.