next up previous
Next: The Metropolis-Hastings Update Up: Combination of Updates Previous: Reversibility

State-Dependent Mixing

Green (1995) proposed an algorithm that involves state-dependent mixing having mixing probabilities that depend on the current state. There are a finite or infinite set of transition kernels tex2html_wrap_inline2627 , tex2html_wrap_inline2629 , and state-dependent mixing probabilities tex2html_wrap_inline2631 . The overall transition kernel is

displaymath2633

To make a move when at x, we choose kernel tex2html_wrap_inline2547 with probability tex2html_wrap_inline2631 and then simulate the next state with probability tex2html_wrap_inline2641 .

No nice properties are transferred from the tex2html_wrap_inline2547 to P, so we introduce the substochastic kernels tex2html_wrap_inline2647 . Then tex2html_wrap_inline2649 , so if all of the tex2html_wrap_inline2651 are reversible with respect to tex2html_wrap_inline2523 , then so is P, and if P is stochastic, then tex2html_wrap_inline2523 is a stationary distribution of P. Note that each tex2html_wrap_inline2651 determines tex2html_wrap_inline2561 and tex2html_wrap_inline2547 through

  equation158

and

  equation161

so we may consider that we have been given the tex2html_wrap_inline2651 to specify the algorithm.

A simple trick allows us to use the same argument when we are given a set of substochastic kernels tex2html_wrap_inline2671 , tex2html_wrap_inline2629 , that sum to a substochastic kernel, that is

  equation166

Define the defect

displaymath2675

and a new kernel

displaymath2677

where I is the identity kernel, defined by

displaymath2681

Then tex2html_wrap_inline2683 is reversible with respect to any distribution tex2html_wrap_inline2523 since

displaymath2687

is trivially symmetric under the interchange of u and v. If we add tex2html_wrap_inline2683 to our set of kernels, then the sum is stochastic.

Thus we have the following formulation of state-dependent mixing. Suppose we are given a set of substochastic kernels satisfying (1.7). Then the following combined update

  1. Choose with i probability tex2html_wrap_inline2631 defined by (1.5). With probability tex2html_wrap_inline2699 , skip step 2 and stay at the current position.
  2. Simulate a new value of x from the probability distribution tex2html_wrap_inline2641 defined by (1.6).
has the stochastic transition kernel tex2html_wrap_inline2705 and is reversible with respect to tex2html_wrap_inline2523 if each of the tex2html_wrap_inline2651 is reversible with respect to tex2html_wrap_inline2523 .


next up previous
Next: The Metropolis-Hastings Update Up: Combination of Updates Previous: Reversibility

Charles Geyer
Fri Jul 5 15:26:21 CDT 1996