This is not a trivial problem, and the procedure is rife with opportunities for self-delusion. For example, one can take samples of about 50 drops from what is actually an exponential DSD; hypothesize that the DSD is gamma; use the third (LWC), fourth (bracketing R) and sixth (Z) moments to "fit" gamma parameters; find a mean value of about 6 for the gamma shape parameter; and thus convince oneself, quite erroneously, that the DSD is gamma after all. This and other findings have been obtained by Monte Carlo simulations of sampling from known DSDs, which allows:
(1) The "fitted" distribution functions to be compared with the original distribution from which the samples were drawn (something that can never be done in practice) as well as with the sample distributions;
(2) The process to be repeated to assess the effects of sampling variability; and
(3) Aspects (such as the bias) not accessible through error propagation analysis to be investigated.
Some general conclusions from the work thus far:
(a) The bias is greater when higher-order moments are employed (implication: use of Z in the procedure is especially unwise).
(b) The size of the largest drop in the sample is an excellent predictor of the calculated gamma shape parameter.
(c) Samples in the hundreds of drops are required to reduce the bias to any manageable size.