ABSTRACT 



Describes a method for pps sampling developed by H. O. Hartley 

 (1966). To facilitate drawing samples using this method, a computer 

 subroutine is provided that identifies the units in the sample, calculates 

 their sampling probability, and summarizes the distribution of this 

 probability over the entire population for use in subsequent estimates 

 of the population total and its variance. An estimate of the error 

 variance is derived in the case that the ratio of measured value to size 

 is a polynomial function of size. This error estimate is a multivariate 

 analog of Hartley's linear model for the superpopulation. To permit 

 sampling without replacement, excessively large units in the population 

 are assigned to a stratum to be measured with absolute certainty. 



