Write , where has the standard normal distribution. Apply Gaussian integration by parts to :
Since ,
The polynomial growth assumption and decay of the normal distribution density justify the integration by parts.
Condition on and use the characteristic function of a standard normal distribution:
The integrand has modulus one, so Fubini's theorem is immediate. Taking expectation over yields
The intended monotone rearrangement is
Here is the quantile function of , agreeing with the ordinary inverse when is continuous and strictly increasing. With an atomless measure , its cumulative distribution function is continuous, and the probability integral transform makes uniform on for . Thus . The one-dimensional monotone rearrangement theorem says this transport map minimizes the cost for convex continuous , whenever the cost integrals are well defined. Values at exceptional endpoints may be chosen arbitrarily.
The printed assumptions omit an essential source condition. Invertibility of alone does not ensure an admissible transport map. For example, and a standard normal distribution satisfy the stated condition on , but is always a Dirac measure. There is no solution to the Monge optimal transport problem in this example. The boxed answer therefore requires the additional assumption that is an atomless measure, or an equivalent condition making the displayed map admissible. For arbitrary sources the always admissible monotone transport plan is , which need not be induced by a map.
Since and a standard normal variable has second and fourth moments and ,
Thus has the chi-squared distribution with degrees of freedom and
Convergence in distribution means that the cumulative distribution functions satisfy at every continuity point of . The Lévy continuity theorem states that this is equivalent to pointwise convergence of characteristic functions to the characteristic function of , provided the limiting function is continuous at zero.
For fixed , independence and a Gaussian integral give
Using ,
Therefore , the characteristic function of the standard normal distribution. Lévy's theorem now gives
This derives the limit directly, without assuming the central limit theorem.
Let be the likelihood function. Since the prior distribution is , the posterior density has the form
where is the normalizing constant.
Choose and use the Gaussian autoregressive proposal reversible with respect to a standard normal distribution
Thus
whose covariance matrix is the required . Reversibility with respect to the standard normal distribution says
so
The Metropolis–Hastings acceptance probability therefore reduces to
This is the Preconditioned Crank–Nicolson algorithm with proposal scale .
Let . Since has an everywhere positive density, it is continuous and strictly increasing, so . Membership in the Kolmogorov neighborhood of a distribution gives
By the symmetry of the standard normal distribution,
The stated asymptotic-bias formula for the sample median therefore yields
If , then has a scaled half-normal distribution with
Put . The population median is , and the asymptotic distribution of a sample median gives
Under , positivity of the scale means . If is the -quantile of the standard normal distribution, an asymptotically level- test rejects for
Larger scale makes the population median larger, so this is the appropriate one-sided rejection region.
Let be the corresponding standard normal distribution quantile. Slutsky theorem permits replacing by the consistent estimator , giving the confidence interval
After standardization, the preceding convergence in distribution and continuity of the normal distribution show that its coverage tends to .
The commands choose observations and predictors; generate a matrix of independent standard-normal values; generate independent binary random variables with success probability ; and report the observed number of successes, . The first glm call then fits the Bernoulli logistic-regression model
The prediction command returns the fitted response probabilities and sums them.
The likelihood is
At the maximum-likelihood estimator, the intercept component of the score function is
Therefore
exactly, by the fitted-mean balance for logistic regression with an intercept.
The second fit uses probit regression,
where is the standard normal distribution function. Its intercept score is
a weighted residual equation. It does not imply , so exact equality is not expected. The output should nevertheless be close to : the data were generated with constant success probability independently of , so the fitted slopes should be small, the linear predictors should cluster near a common intercept, and the score weights should be nearly constant.
Put . The standardized efficient score function for scale is
An optimal B-robust scale M-estimator solves
where the optimal bounded score has the clipped-score form
The constants enforce Fisher consistency, the chosen normalization, and the clipping bound on the influence function.
For the standard normal distribution, , hence . The score therefore simplifies to
with chosen so that for .
Let be independent random variables with the standard normal distribution. The unconstrained Milstein method step is
The noise is additive, so its Milstein correction vanishes. The diagonal, independent noise fields commute, so no cross iterated stochastic integral is required.
Implement each reflecting boundary condition for a diffusion by folding the proposed coordinate back into . One formula that also handles multiple overshoots is
The complete update is
Under the usual smoothness assumptions, the Milstein discretization has strong order of convergence one and weak order of convergence one. The reflection enforces the boundary pathwise.
The probability density function must integrate to one, so
Also, linearity of expectation and the mean and variance of the standard normal distribution give
Using the probability density function to compute this expected value gives
These are the values and of the gamma function.
Since is nonnegative, its cumulative distribution function is zero for . For ,
The last equality uses the symmetry of the standard normal distribution. Thus
Apply the central limit theorem to the independent and identically distributed random variables , whose mean and variance are and . It gives the approximate standard normal distribution
The supplied standard normal quantile leaves probability in each tail. Hence
Since , this event is equivalent to
For , equivalently integer , both denominators below are positive. Inverting the inequalities gives the approximate confidence interval
The approximation is justified for large ; positivity at by itself does not guarantee accurate coverage. This construction is a confidence interval from a relative normal pivot: the unknown scale remains inside the pivotal ratio until the inequalities are inverted. Its meaning is that the random interval covers the fixed statistical parameter with probability approximately under repeated sampling.
The sample mean is exactly , so has the standard normal distribution. The given quantile gives probability to the interval . Inverting this pivot yields the exact confidence interval
Let and . The standard normal distribution is symmetric and assigns no probability mass to zero. Thus , and its equal two-sided tails give
Therefore
so the two events are independent. The permitted approximation gives and , but independence is exact.
An expected value of a signed random variable is finite precisely when its absolute value is integrable. With the standard normal distribution probability density function,
For , this integral equals . For , its integrand does not decay sufficiently and the integral diverges. Therefore
The value is zero by the odd function and the symmetric distribution. For the other , both positive and negative parts have infinite integrals; formal cancellation would give only a Cauchy principal value, not an expected value. This distinction is the integrability threshold for exponential powers of a Gaussian.
The standard normal distribution function is strictly increasing. Thus implies
Rearranging gives
This is the confidence-interval/test duality for the matching one-sided procedures. The strict inequalities agree with the printed open interval.
Let and . For independent and identically distributed random variables, the central limit theorem states
The limit is the standard normal distribution. If the variance is zero, then almost surely; the centered sums are zero, and division by is not defined.
By part (c) with , let for independent and identically distributed exponential random variables. The requested integral is . Each summand has expectation one and variance one, so the central limit theorem gives
Zero is a continuity point of the standard normal distribution, and consequently