Unbiased Estimation of The Reciprocal Mean For Non-Negative Random Variables

In recent years, Monte Carlo estimators have been proposed that can estimate the ratio of two expectations without bias. We investigate the theoretical properties of a Taylor-expansion based estimator of the reciprocal mean of a non-negative random variable. We establish explicit expressions for the computational efficiency of this estimator and obtain optimal choices for its parameters. We also derive corresponding practical confidence intervals and show that they are asymptotically equivalent to the maximum likelihood (biased) ratio estimator as the simulation budget increases.

[1]  Don McLeish,et al.  A general method for debiasing a Monte Carlo estimator , 2010, Monte Carlo Methods Appl..

[2]  Nan Chen,et al.  Unbiased Monte Carlo computation of smooth functions of expectations via Taylor expansions , 2015, 2015 Winter Simulation Conference (WSC).

[3]  Jesper Møller,et al.  An Introduction to Simulation-Based Inference for Spatial Point Processes , 2003 .

[4]  Peter W. Glynn,et al.  Unbiased Estimation with Square Root Convergence for SDE Models , 2015, Oper. Res..

[5]  Sandeep Juneja,et al.  Regenerative Simulation for Queueing Networks with Exponential or Heavier Tail Arrival Distributions , 2013, ACM Trans. Model. Comput. Simul..

[6]  I. Lux Monte Carlo Particle Transport Methods: Neutron and Photon Calculations , 1991 .

[7]  Reuven Y. Rubinstein,et al.  Simulation and the Monte Carlo method , 1981, Wiley series in probability and mathematical statistics.

[8]  W. Richter Limit Theorems for Sequences of Random Variables with Sequences of Random Indeces , 1965 .

[9]  Yves F. Atchad'e,et al.  On Russian Roulette Estimates for Bayesian Inference with Doubly-Intractable Likelihoods , 2013, 1306.4032.

[10]  P. Jacob,et al.  On nonnegative unbiased estimators , 2013, 1309.6473.

[11]  Raymond L. Murray,et al.  Particle-transport simulation with the Monte Carlo method By and . Technical Information Center, ERDA, 1975, (TID-26607) 115 pp. $dollar;5.45 , 1977 .

[12]  Colin Wei,et al.  Markov Chain Truncation for Doubly-Intractable Inference , 2016, AISTATS.

[13]  B. L. Nelson Chapter 16 Simulation Algorithms for Regenerative Processes , 2010 .

[14]  Peter W. Glynn,et al.  Exact estimation for Markov chain equilibrium expectations , 2014, Journal of Applied Probability.

[15]  Aarnout Brombacher,et al.  Probability... , 2009, Qual. Reliab. Eng. Int..

[16]  C. Andrieu,et al.  The pseudo-marginal approach for efficient Monte Carlo computations , 2009, 0903.5480.

[17]  M. Beaumont Estimation of population growth or decline in genetically monitored populations. , 2003, Genetics.

[18]  Peter W. Glynn,et al.  Unbiased Monte Carlo for optimization and functions of expectations via multi-level randomization , 2015, 2015 Winter Simulation Conference (WSC).

[19]  Ward Whitt,et al.  The Asymptotic Efficiency of Simulation Estimators , 1992, Oper. Res..

[20]  Stephen G. Walker,et al.  Posterior Sampling When the Normalizing Constant is Unknown , 2011, Commun. Stat. Simul. Comput..