A procedure for interval estimation of the shape parameter of the gamma distribution is proposed which enables the experimenter to obtain upper and lower confidence limits after each observation is sampled. An approximate lower bound for the probability that the resulting sequence of confidence intervals will contain the shape parameter is obtained. The proposed procedure may be useful in estimating the coefficient of variationassociated with non-negative variables whose distribution is skewed to the right. Comparisons with a single sample procedure are discussed.