# zbMATH — the first resource for mathematics

##### Examples
 Geometry Search for the term Geometry in any field. Queries are case-independent. Funct* Wildcard queries are specified by * (e.g. functions, functorial, etc.). Otherwise the search is exact. "Topological group" Phrases (multi-words) should be set in "straight quotation marks". au: Bourbaki & ti: Algebra Search for author and title. The and-operator & is default and can be omitted. Chebyshev | Tschebyscheff The or-operator | allows to search for Chebyshev or Tschebyscheff. "Quasi* map*" py: 1989 The resulting documents have publication year 1989. so: Eur* J* Mat* Soc* cc: 14 Search for publications in a particular source with a Mathematics Subject Classification code (cc) in 14. "Partial diff* eq*" ! elliptic The not-operator ! eliminates all results containing the word elliptic. dt: b & au: Hilbert The document type is set to books; alternatively: j for journal articles, a for book articles. py: 2000-2015 cc: (94A | 11T) Number ranges are accepted. Terms can be grouped within (parentheses). la: chinese Find documents in a given language. ISO 639-1 language codes can also be used.

##### Operators
 a & b logic and a | b logic or !ab logic not abc* right wildcard "ab c" phrase (ab c) parentheses
##### Fields
 any anywhere an internal document identifier au author, editor ai internal author identifier ti title la language so source ab review, abstract py publication year rv reviewer cc MSC code ut uncontrolled term dt document type (j: journal article; b: book; a: book article)
Simulating ratios of normalizing constants via a simple identity: A theoretical exploration. (English) Zbl 0857.62017

Summary: Let ${p}_{i}\left(w\right)$, $i=1,2$, be two densities with common support where each density is known up to a normalizing constant: ${p}_{i}\left(w\right)={q}_{i}\left(w\right)/{c}_{i}$. We have draws from each density (e.g., via Markov chain Monte Carlo), and we want to use these draws to simulate the ratio of the normalizing constants, ${c}_{1}/{c}_{2}$. Such a computational problem is often encountered in likelihood and Bayesian inference, and arises in fields such as physics and genetics. Many methods proposed in statistical and other literature (e.g., computational physics) for dealing with this problem are based on various special cases of the following simple identity:

${c}_{1}/{c}_{2}={E}_{2}\left[{q}_{1}\left(w\right)\alpha \left(w\right)\right]/{E}_{1}\left[{q}_{2}\left(w\right)\alpha \left(w\right)\right]·$

Here ${E}_{i}$ denotes the expectation with respect to ${p}_{i}$ $\left(i=1,2\right)$, and $\alpha$ is an arbitrary function such that the denominator is non-zero. A main purpose of this paper is to provide a theoretical study of the usefulness of this identity, with focus on (asymptotically) optimal and practical choices of $\alpha$.

Using a simple but informative example, we demonstrate that with sensible (not necessarily optimal) choices of $\alpha$, we can reduce the simulation error by orders of magnitude, when compared to the conventional importance sampling method, which corresponds to $\alpha =1/{q}_{2}$. We also introduce several generalizations of this identity for handling more complicated settings (e.g., estimating several ratios simultaneously) and pose several open problems that appear to have practical as well as theoretical value. Furthermore, we discuss related theoretical and empirical work.

##### MSC:
 62E99 Statistical distribution theory 65C99 Probabilistic methods, simulation and stochastic differential equations (numerical analysis)