The occurrence frequency distributions (size distributions) are the most important diagnostics for self-organized criticality systems. There are at least three formats for size distributions: (i) the differential size distribution function, (ii) the cumulative size distribution function, and (iii) the rank-order plot. Each of the three formats (or methods) has at least three ranges of event sizes: (i) a range with statistically incomplete sampling; (ii) an inertial range or power law fitting range with statistically complete sampling; and (iii) a range bordering finite system sizes. Only the intermediate range with power law behavior should be used to determine the power law slope from fitting the observed size distributions. The establishment of power law functions in a given observed size distribution depends crucially on the choice of the fitting range, which should have a logarithmic range of at least 2–3 decades. Often the fitted distribution functions exhibit significant deviations from an ideal power law and can be fitted better with alternative functions, such as log-normal distributions, Pareto type-II distributions, and Weibull distributions.