A method for fitting parameters of the gamma distribution to data containing some zero values using maximum likelihood methods is presented. The procedure is based on a conceptual model of the data ...