Statistica Sinica

Kim-Hung Li

Abstract:The sampling/importance resampling algorithm is an approximate noniterative sampling method. The algorithm has been used on many occasions to select an approximate random sample of size from a target distribution from input random variates. The selection mechanism is an unequal probability sampling with weights being the importance weights. As the weights are random, sampling without replacement is not always possible and some input variates may have more than one copy in the final sample.Duplication of values in the final output is undesirable as it means dependence among the output variates. In this paper a general and simple determination rule for is proposed. It keeps the duplication problem at a tolerably low level when a tight resampling method is used. We show that (a) if and only if the importance weight is bounded above, (b) if the importance weight has a moment generating function, the suggested is of order , and (c) may need to be as large as if the importance weight has finite -th moment for a . A procedure is suggested to determine numerically. The method is tested on the Pareto, Gamma and Beta distributions, and gives satisfactory results.

Key words and phrases:Importance weight, Monte Carlo sampling, resampling method, sample size, tight resampling algorithm.