{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,17]],"date-time":"2026-08-17T23:14:58Z","timestamp":1787008498625,"version":"build-2736575974"},"reference-count":22,"publisher":"Oxford University Press (OUP)","issue":"3","license":[{"start":{"date-parts":[[2019,4,8]],"date-time":"2019-04-08T00:00:00Z","timestamp":1554681600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100012338","name":"Alan Turing Institute","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100012338","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000266","name":"Engineering and Physical Sciences Research Council","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100000266","id-type":"DOI","asserted-by":"publisher"}]},{"name":"European Research Council"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,9,1]]},"abstract":"<jats:title>Summary<\/jats:title>\n                  <jats:p>We consider the problem of approximating the product of $n$ expectations with respect to a common probability distribution $\\mu$. Such products routinely arise in statistics as values of the likelihood in latent variable models. Motivated by pseudo-marginal Markov chain Monte Carlo schemes, we focus on unbiased estimators of such products. The standard approach is to sample $N$ particles from $\\mu$ and assign each particle to one of the expectations; this is wasteful and typically requires the number of particles to grow quadratically with the number of expectations. We propose an alternative estimator that approximates each expectation using most of the particles while preserving unbiasedness, which is computationally more efficient when the cost of simulations greatly exceeds the cost of likelihood evaluations. We carefully study the properties of our proposed estimator, showing that in latent variable contexts it needs only ${O} (n)$ particles to match the performance of the standard approach with ${O}(n^{2})$ particles. We demonstrate the procedure on two latent variable examples from approximate Bayesian computation and single-cell gene expression analysis, observing computational gains by factors of about 25 and 450, respectively.<\/jats:p>","DOI":"10.1093\/biomet\/asz008","type":"journal-article","created":{"date-parts":[[2019,2,17]],"date-time":"2019-02-17T23:22:18Z","timestamp":1550445738000},"page":"708-715","source":"Crossref","is-referenced-by-count":7,"title":["Unbiased approximations of products of expectations"],"prefix":"10.1093","volume":"106","author":[{"given":"A","family":"Lee","sequence":"first","affiliation":[{"name":"School of Mathematics, University of Bristol, University Walk, Bristol BS8 1TW, UK"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"S","family":"Tiberi","sequence":"additional","affiliation":[{"name":"Institute of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"G","family":"Zanella","sequence":"additional","affiliation":[{"name":"Department of Decision Sciences, BIDSA and IGIER, Bocconi University, Via Roentgen 1, 20136 Milan, Italy"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2019,4,8]]},"reference":[{"key":"2019091012043431500_B1","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1007\/s11222-008-9083-x","article-title":"Bayesian estimation of quantile distributions","volume":"19","author":"Allingham","year":"2009","journal-title":"Statist. Comp."},{"key":"2019091012043431500_B2","doi-asserted-by":"crossref","first-page":"697","DOI":"10.1214\/07-AOS574","article-title":"The pseudo-marginal approach for efficient Monte Carlo computations","volume":"37","author":"Andrieu","year":"2009","journal-title":"Ann. Statist."},{"key":"2019091012043431500_B3","doi-asserted-by":"crossref","first-page":"2661","DOI":"10.1214\/15-AAP1158","article-title":"Establishing some order amongst exact approximations of MCMCs","volume":"26","author":"Andrieu","year":"2016","journal-title":"Ann. Appl. Prob."},{"key":"2019091012043431500_B4","doi-asserted-by":"crossref","first-page":"1139","DOI":"10.1093\/genetics\/164.3.1139","article-title":"Estimation of population growth or decline in genetically monitored populations","volume":"164","author":"Beaumont","year":"2003","journal-title":"Genetics"},{"key":"2019091012043431500_B5","doi-asserted-by":"crossref","first-page":"235","DOI":"10.1111\/rssb.12312","article-title":"Approximate Bayesian computation with the Wasserstein distance","volume":"81","author":"Bernton","year":"2019","journal-title":"J. R. Statist. Soc. B"},{"key":"2019091012043431500_B6","doi-asserted-by":"crossref","first-page":"839","DOI":"10.1111\/rssb.12280","article-title":"The correlated pseudomarginal method","volume":"80","author":"Deligiannidis","year":"2018","journal-title":"J. R. Statist. Soc. B"},{"key":"2019091012043431500_B7","doi-asserted-by":"crossref","first-page":"295","DOI":"10.1093\/biomet\/asu075","article-title":"Efficient implementation of Markov chain Monte Carlo when using an unbiased likelihood estimator","volume":"102","author":"Doucet","year":"2015","journal-title":"Biometrika"},{"key":"2019091012043431500_B8","doi-asserted-by":"crossref","first-page":"419","DOI":"10.1111\/j.1467-9868.2011.01010.x","article-title":"Constructing summary statistics for approximate Bayesian computation: Semi-automatic approximate Bayesian computation","volume":"74","author":"Fearnhead","year":"2012","journal-title":"J. R. Statist. Soc. B"},{"key":"2019091012043431500_B9","doi-asserted-by":"crossref","first-page":"1559","DOI":"10.1214\/105051604000000396","article-title":"Concentration of permanent estimators for certain large matrices","volume":"14","author":"Friedland","year":"2004","journal-title":"Ann. Appl. Prob."},{"key":"2019091012043431500_B10","first-page":"241","article-title":"On the matching polynomial of a graph","volume-title":"Algebraic Methods in Graph Theory","author":"Godsil","year":"1981"},{"key":"2019091012043431500_B11","doi-asserted-by":"crossref","first-page":"223","DOI":"10.2307\/3318737","article-title":"An adaptive Metropolis algorithm","volume":"7","author":"Haario","year":"2001","journal-title":"Bernoulli"},{"key":"2019091012043431500_B12","doi-asserted-by":"crossref","first-page":"635","DOI":"10.1093\/biomet\/asp036","article-title":"Approximating the $\\alpha$-permanent","volume":"96","author":"Kou","year":"2009","journal-title":"Biometrika"},{"key":"2019091012043431500_B13","doi-asserted-by":"crossref","first-page":"749","DOI":"10.1007\/BF02366856","article-title":"Computing the permanent by importance sampling method","volume":"32","author":"Kuznetsov","year":"1996","journal-title":"Cybernet. Syst. Anal."},{"key":"2019091012043431500_B14","doi-asserted-by":"crossref","first-page":"1167","DOI":"10.1007\/s11222-011-9288-2","article-title":"Approximate Bayesian computational methods","volume":"22","author":"Marin","year":"2012","journal-title":"Statist. Comp."},{"key":"2019091012043431500_B15","doi-asserted-by":"crossref","first-page":"222","DOI":"10.1006\/tpbi.1995.1027","article-title":"Markovian modeling of gene-product synthesis","volume":"48","author":"Peccoud","year":"1995","journal-title":"Theor. Pop. Biol."},{"key":"2019091012043431500_B16","doi-asserted-by":"crossref","volume-title":"Combinatorial Mathematics","author":"Ryser","DOI":"10.5948\/UPO9781614440147"},{"key":"2019091012043431500_B17","doi-asserted-by":"crossref","DOI":"10.1007\/978-0-387-34675-5","volume-title":"Stochastic Orders","author":"Shaked","year":"2007"},{"key":"2019091012043431500_B18","doi-asserted-by":"crossref","first-page":"727","DOI":"10.1093\/biomet\/asx031","article-title":"Pseudo-marginal Metropolis\u2013Hastings sampling using averages of unbiased estimators","volume":"104","author":"Sherlock","year":"2017","journal-title":"Biometrika"},{"key":"2019091012043431500_B19","doi-asserted-by":"crossref","first-page":"238","DOI":"10.1214\/14-AOS1278","article-title":"On the efficiency of pseudo-marginal random walk Metropolis algorithms","volume":"43","author":"Sherlock","year":"2015","journal-title":"Ann. Statist."},{"key":"2019091012043431500_B20","doi-asserted-by":"crossref","first-page":"i647","DOI":"10.1093\/bioinformatics\/bty568","article-title":"Bayesian inference on stochastic gene transcription from flow cytometry data","volume":"34","author":"Tiberi","year":"2018","journal-title":"Bioinformatics"},{"key":"2019091012043431500_B21","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1016\/0304-3975(79)90044-6","article-title":"The complexity of computing the permanent","volume":"8","author":"Valiant","year":"1979","journal-title":"Theor. Comp. Sci."},{"key":"2019091012043431500_B22","doi-asserted-by":"crossref","first-page":"231","DOI":"10.1007\/s11222-014-9491-z","article-title":"Monte Carlo algorithms for computing $\\alpha$-permanents","volume":"26","author":"Wang","year":"2016","journal-title":"Statist. Comp."}],"container-title":["Biometrika"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/biomet\/article-pdf\/106\/3\/708\/29133126\/asz008.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,9,11]],"date-time":"2022-09-11T20:24:43Z","timestamp":1662927883000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biomet\/article\/106\/3\/708\/5431312"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2019,4,8]]},"references-count":22,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2019,4,8]]},"published-print":{"date-parts":[[2019,9,1]]}},"URL":"https:\/\/doi.org\/10.1093\/biomet\/asz008","relation":{},"ISSN":["0006-3444","1464-3510"],"issn-type":[{"value":"0006-3444","type":"print"},{"value":"1464-3510","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2019,9]]},"published":{"date-parts":[[2019,4,8]]}}}