The printed identity needs an expectation on its right-hand side. For fixed , condition on and use independence of :
The expression inside this expectation is generally random and cannot equal the unconditional probability by itself.
For each independent Monte Carlo method replicate, draw from its given distribution. If , draw independent -variables, form and set . For the empty sum is . If may equal zero, use when , since . The law of total expectation gives
Since , the estimator has statistical consistency by the strong law of large numbers. It is the conditional Monte Carlo estimator obtained by Rao-Blackwellization of the direct indicator , and the law of total variance gives .
This version assumes that the cumulative distribution function can be evaluated. An easy sampler alone does not automatically supply an easy cumulative distribution function evaluation. If only sampling is available, use the direct indicator estimator, or replace each conditional cumulative distribution function by the average of several independent indicators ; the nested version remains an unbiased estimator. With such inner draws its variance per outer replicate is .
Use the ergodic theorem for a positive Harris recurrent Markov chain: for an invariant distribution and an integrable function , the time average converges almost surely to under the usual ergodicity conditions. This is the relevant theorem for correlated Metropolis–Hastings algorithm output, rather than the strong law of large numbers. Here is integrable. Indeed, putting
we have , which bounds the unnormalized density by a constant times an integrable Gaussian density and also gives finite first moments.
For a proposal scale , initialize anywhere in . At each step draw independent random variables with the standard normal distribution and an independent uniform , propose , and accept when
Otherwise set . The normal proposal is symmetric, so its density cancels; the target normalizing constant also cancels. Include every state in the average, including repeated states after rejection. The ergodic theorem for a positive Harris recurrent Markov chain then gives
A fixed discarded initial segment does not change this limit.
Small proposal variance gives high acceptance but tiny moves and strong serial dependence. Large variance gives more ambitious moves, but many proposals enter very low-density regions and are rejected, creating long runs at one state. An intermediate scale should be assessed by exploration and effective sample size of a Markov chain, not acceptance rate alone.
This density has two modes, at and their interchange. To see this, stationary points satisfy , implying either or . The unequal solutions have and are local minima of ; the equal solution is a saddle. Mode switching is therefore a material part of proposal-scale selection: a chain confined to one mode can have high acceptance and a misleading finite-run estimate.
Figure 1.
Two modes of the quartic target density, with the diagonal saddle between them, explaining slow mode switching in random-walk Metropolis-Hastings
. Details.
The observations have finite mean and variance. The strong law of large numbers gives almost surely, so almost surely and therefore is a consistent estimator. The central limit theorem gives
A renewal process has independent identically distributed, nonnegative, almost surely finite inter-arrival times , epochs , and count . Allowing zero inter-arrivals is harmless provided . The hypothesis ensures there is with .
The strong law of large numbers for the indicators gives , so is finite for finite . Conversely every fixed is finite, so eventually . On their countable common probability-one event, this proves .
For fixed , put . The hinted comparison process has epochs , hence . Equivalently,
If , for this is bounded by , by summing the first binomial distribution. Thus for . If , deterministically. Therefore a strictly positive exponential moment exists for every fixed finite , with the same allowed for all such . Almost-sure finiteness of the inter-arrivals is implicit in a proper renewal process; allowing an infinite inter-arrival with positive probability would invalidate the divergence claim.
The sigma-algebras decrease with , and (b) gives . The permitted reverse martingale convergence theorem therefore gives an integrable random variable such that almost surely and with convergence in L1. In particular,
It remains to prove that is constant. We must not assume that is the tail sigma-algebra of the ; instead use tail measurability of limits of sample averages directly. Define the random variable valued in the extended real numbers . For every fixed ,
because pointwise. Thus is measurable with respect to for every , and hence to the tail sigma-algebra. Also almost surely.
By the Kolmogorov zero-one law, each tail event , for a rational number , has probability zero or one. Since is finite almost surely, its distribution function can only be that of a constant : taking and using the density of the rational numbers gives almost surely. Its expected value identifies . This proves the strong law of large numbers:
Split the adjacent products into two subsequences:
Within each subsequence the pairs use disjoint coordinates, so each subsequence consists of independent and identically distributed random variables. The two subsequences need not be independent of each other. By independence within a pair,
Apply the strong law of large numbers from (c) separately to obtain and simultaneously almost surely. Since
and both prefactors tend to , the strong law for adjacent products gives
The maximum-likelihood estimator is
Under the score-outer-product convention stated in the question,
This is the sample covariance matrix with divisor . The strong law of large numbers, applied entrywise, gives
Thus the empirical score outer-product information consistently estimates the per-observation Fisher information.
The strong law of large numbers gives
Because , it follows that almost surely. If this martingale were uniformly integrable, almost-sure convergence would imply convergence in , and therefore
But the martingale has constant expectation . This contradiction proves that
Fix and choose
Under , the increments remain independent and identically distributed, with mean . The strong law of large numbers therefore gives
Part (b) implies
Let and then . Since and is continuous on ,
The same argument includes by taking .
Conditional on , the holding time has an exponential distribution of rate , and after that holding time the next state is with transition probability . The memoryless property of exponential random variables and the Markov property of therefore make a continuous-time Markov chain.
Moreover , so . The strong law of large numbers gives almost surely, hence and the process is nonexplosive. Its Q-matrix is
Equivalently, if is the diagonal matrix with entries , then
The strong law of large numbers gives almost surely, and hence in probability. Therefore
by the continuous mapping theorem. Thus the Bayes estimator has consistency.
It differs from the maximum-likelihood estimator by
Since , the numerator is and the denominator is , so this difference is . The limit in part a and Slutsky theorem now give
The strong law of large numbers gives
A continuous one-to-one real function is strictly monotone and has a continuous inverse on its image. Hence
whenever lies in that image; an arbitrary extension outside the image handles the remaining finite-sample values. The continuous mapping theorem gives almost surely, so this is a consistent estimator.
Here the total holding rate is
and therefore the jump chain has state-independent transitions
It is an upward biased random walk on . By the strong law of large numbers, almost surely, so . The jump chain, and hence , is transient.
Counting measure is invariant for this jump chain: for every , its incoming mass is . The invariant-measure transfer between a jump chain and a CTMC therefore gives an invariant measure proportional to . Since
the normalized invariant distribution of a continuous-time Markov chain is
This does not contradict transience because an invariant distribution of an explosive chain need not imply positive recurrence.
It remains to verify explosion. Starting from zero, the probability that the biased walk ever reaches is . Once it reaches any state, its probability of returning after departure is
so its expected total number of visits to that state, conditional on reaching it, is . If is the number of visits to , then
Consequently the expected sum of all holding times is
The nonnegative random variable is therefore finite almost surely. By explosion by summable holding times, makes infinitely many jumps by time , so it is explosive.
The strong law of large numbers states that if are independent and identically distributed random variables with , then
almost surely. Under the stronger fourth-moment hypothesis, it has a short proof. Put , , and . Independence and centering make every term in the expansion of vanish unless each index occurs at least twice, so
For every , the Markov inequality gives
The probabilities are summable, so the first Borel-Cantelli lemma says that occurs only finitely often almost surely. Apply this simultaneously to to obtain almost surely, which is the claimed law.
By the strong law of large numbers,
almost surely. Brownian scaling and a maximal inequality show that changing Brownian time by changes its value by ; explicitly, first restrict to , bound the Brownian maximum over a time interval of length , and then let . Consequently
in probability.
But
for every . Since has the law of , Slutsky theorem proves the Central limit theorem from the Skorokhod embedding:
The strong law of large numbers states that for independent identically distributed integrable random variables,
To prove it, set . The tail-sum formula gives , so the Borel-Cantelli lemmas make the two sequences eventually equal. Also
The Kolmogorov convergence theorem implies that converges almost surely, and Kronecker lemma yields
Finally , so the Cesàro mean of these expectations tends to .
Using Bayes' theorem and the fact that the prior is proper,
Thus is an unbiased estimator of , and the Harmonic mean estimator of Bayesian model evidence is . The reciprocal is not itself generally unbiased, though it is consistent when the strong law of large numbers applies.
Convergence in distribution makes tight. Hence for every there is such that for all sufficiently large . Since ,
Letting gives . Also, the continuous-mapping theorem gives . Since
the Slutsky theorem yields
Now write the requested statistic as
The central limit theorem gives
Since , the strong law of large numbers gives
almost surely, and therefore its reciprocal converges in probability to one. Applying the product result proves
This is a self-normalized central limit theorem with a second-moment denominator.
Assume that wherever on the part of the space relevant to . The importance sampling algorithm draws independent and identically distributed random variables from , forms the weights
and returns
A change of measure gives
Hence the estimator is unbiased whenever the expectation exists, and, if , the strong law of large numbers gives almost surely.
Conditionally on , the cumulative distribution function of is the weighted empirical distribution
For fixed , the strong law of large numbers gives, almost surely,
while
The ratio therefore converges almost surely to . This is precisely the claimed convergence in distribution of the self-normalized importance sampling resample.