No label defined (Q4527097)
From MaRDI portal
!
This is the item page for this Wikibase entity, intended for internal use and editing purposes. Please use the normal view instead:
scientific article; zbMATH DE number 1560252
| Language | Label | Description | Also known as |
|---|---|---|---|
| default for all languages | No label defined |
||
| English | No label defined |
scientific article; zbMATH DE number 1560252 |
Statements
2 October 2001
0 references
Markov chain Monte Carlo
0 references
law of large numbers
0 references
central limit theorem
0 references
Harris recurrent kernels
0 references
Hilbert space theory
0 references
reversible and non-reversible kernels
0 references
asymptotic efficiency
0 references
0.8314083814620972
0 references
0.8211519122123718
0 references
0.8138890862464905
0 references
0.7990028858184814
0 references
0.7984508275985718
0 references
In the introduction some general remarks are presented. Markov chain Monte Carlo (MCMC) simulates a probability distribution \(\pi\) by running a Markov chain \(X_1,X_2,\dots\) having \(\pi\) as its unique stationary distributions. An expectation \[ \mu=E_\pi \bigl\{f(x) \bigr\}=\int f(x)\pi(dx) \tag{1} \] is approximated by the sample mean over the path of the Markov chain \[ \widehat\mu_n= {1\over n}\bigl(f(x_1) +\cdots+f (x_n)\bigr).\tag{2} \] This is similar to the classical independent sample Monte Carlo (ISMC), which uses a sample mean of the same form as (2) but with \(X_1,X_2,\dots\) independent and identically distributed \(\pi\). The law of large numbers (LLN) and the central limit theorem (CLT) hold. They guarantee that \(\widehat\mu_n\) converges almost surely to \(\mu\) and \(\sqrt n(\widehat\mu_n -\mu)\) converges in law to some normal random variable. This justifies the use of \(\widehat\mu_n\) as an approximation of \(\mu\) and provides an estimate of the approximation error, just as in ISMC. The difference between MCMC and ISMC is that the variance in CLT is typically larger, sometimes much larger in MCMC than in ISMC. An MCMC algorithm is characterized by its transition probability kernel \(P\) defined by \[ P(x,A)= P(X_{n+1}\in A\mid X_n=x). \] A kernel \(P\) has a probability measure \(\pi\) as its stationary distribution if \(\pi(A)=\int P(x,A)\pi(dx)\) holds for all measurable sets \(A\). The authors use Harris recurrent kernels. In the second part of the note preliminary Hilbert space theory is presented and in the third part a necessary and sufficient condition for efficiency is given. In the fourth part the authors construct the inverse. So, they consider Markov chains on finite state spaces and explain how to construct a generalized inverse that represents the inverse Laplacian. In the fifth section the non-reversible Markov chain introduced by P. Diaconis, S. Holmes and R. M. Neal (1997) is analyzed. In the sixth section comparing the performance of reversible and non-reversible kernels is presented. In the last part entitled ``Discussion'' interesting conclusions for non-reversible Markov chains are considered. In this sense a particular example of Diaconis et al. (1997) applying the criterion of asymptotic efficiency is reanalyzed. For other details see the authors' references. For the entire collection see [Zbl 0951.00044].
0 references