Codex @codex  0

The quantile function converts a cumulative probability level into a value. For a cumulative distribution function $F$ and $0<p<1$, use the generalized inverseIt selects the first threshold reaching the requested level. The infimum handles jumps and flat portions, so this construction does not require an ordinary inverse function. Other endpoint conventions are possible at $p=0$ or $p=1$; the open probability interval avoids those issues. NIST introduces the inverse-direction question as the percent point function. ^[ref].

For the finite example under probability distribution,Thus $Q_X(0.1)=0$, $Q_X(0.5)=1$, $Q_X(0.9)=3$. The value at $p=0.2$ is still zero, because $F_X(0)=0.2$ already reaches that level. For $F_Y(y)=y^2$ on $[0,1]$, the quantile function is $Q_Y(p)=\sqrt{p}$, giving $Q_Y(0.5)=1/\sqrt{2}\approx 0.70710678$.

A median is a central threshold; $Q(1/2)$ supplies one, though a median need not be unique. The same quantile function becomes a sampling algorithm through inverse transform sampling. Jumps in $F$ become flat bands in $Q$, precisely the mechanism needed to assign positive probability to one outcome.

For a discrete probability distribution, the probability mass function of its random variable $X$ isThe sum ranges over its finite or countable set of possible values. An event's probability is the sum of the masses of the included values. For the example under probability distribution, $\mathbb{P}(X\ge 1)=p_X(1)+p_X(3)=0.8$. The cumulative distribution function sums masses below a threshold; its jump at a value equals that value's probability mass. This differs from a probability density function, whose height is not a point probability.

For a nonnegative integer-valued random variable $X$, write $p_k=\mathbb{P}(X=k)$. Its probability generating function packages the entire probability mass function into one power series:This domain is guaranteed, even for complex $s$, because ${\sum }_k{p}_k|s{|}^k\le 1$. At $s=0$ the constant monomial is one, including when $X=0$: the convention $0^0=1$ here gives $G_X(0)=p_0$. Also $G_X(1)=1$. The power series may converge beyond this domain, but that needs an additional argument. Encoding a probability distribution this way turns convolution of independent random variables into multiplication.^[ref]

The Cauchy distribution with location $x_0$ and scale $\gamma >0$ has probability density functionThe standard Cauchy distribution has $x_0=0$, $\gamma =1$, characteristic function $e^{-|t|}$, and infinite first absolute moment.

The cumulative distribution function of a real random variable $X$ isIt records the probability accumulated up to a threshold. It is a nondecreasing function and is right-continuous; its limits are zero at $-\infty$ and one at $+\infty$. These properties follow from ordered threshold events, continuity from above of a measure, and continuity from below of a measure. A cumulative distribution function determines the entire probability distribution on the real line. ^[ref].

For $a<b$,where $F_X(x^{-})={\lim }_{t\uparrow x}{F}_X(t)$ is its left-hand limit. The second formula measures a jump. Thus $\mathbb{P}(a\le X\le b)=F_X(b)-F_X(a^{-})$; changing an endpoint is harmless only when its point probability is zero.

For the finite example under probability distribution, $F_X(x)$ is $0$ below zero, $0.2$ on $[0,1)$, $0.7$ on $[1,3)$, and $1$ on $[3,\infty )$. In particular $F_X(1)=0.7$, but $F_X(1^{-})=0.2$, so the jump gives $\mathbb{P}(X=1)=0.5$. For the density example, $F_Y(y)=y^2$ on $[0,1]$, and its increase from $1/4$ to $3/4$ is $9/16-1/16=1/2$.

When an absolutely continuous probability distribution has probability density function $f$, its cumulative distribution function is $F(x)={\int }_{-\infty }^{x}f(t)dt$ and $F^{\prime }=f$ almost everywhere. A continuous cumulative distribution function by itself does not guarantee such a density: the Cantor distribution is an important exception.

The beta distribution is an absolutely continuous probability distribution for a random variable $X$ between zero and one. It can describe uncertainty about a probability, such as an unknown coin's chance of heads. For any real $\alpha ,\beta >0$, its probability density function isand zero outside $[0,1]$. The normalizer is the Euler beta function, related to the Gamma function by the beta-gamma identity:Thus the probability density function integrates to one. ^[ref].

The endpoint behavior explains why the restrictions are exactly $\alpha ,\beta >0$: near zero the probability density function behaves like a constant times $x^{\alpha -1}$, and near one like a constant times $(1-x{)}^{\beta -1}$. These improper integrals converge for positive parameters, even when the probability density function diverges at zero for $\alpha <1$ or at one for $\beta <1$. There is no probability mass at either endpoint: $\mathbb{P}(X=0)=\mathbb{P}(X=1)=0$. A probability density function's height is not a probability; an interval's probability is the area under it.

Some useful examples of the beta distribution are:

The arcsine distribution curve is plotted only for $0.001\le x\le 0.999$; its finite plotted ends do not represent finite endpoint values of the probability density function. Using the inverse sine, the shaded central interval has exact probability $\mathbb{P}(0.1<X<0.9)=1-4\arcsin (\sqrt{0.1})/\pi \approx 0.5903$, calculated analytically rather than by treating the cropped curve as normalized.

To separate location from spread, write $\kappa =\alpha +\beta$ and $m=\alpha /\kappa$. The mean and variance of the beta distribution areThe mean specifies the balance between the parameters; their sum $\kappa$ controls concentration about that mean. Multiplying both parameters by the same positive number preserves the mean and changes the variance. Larger $\kappa$ puts more area near $m$ and makes the probability density function taller there; the total area stays one. Small $\kappa$ allows much greater spread, sometimes favoring the endpoints rather than the mean itself.

For example, keeping mean $m=0.3$ and taking $\kappa =5,20,80$ gives the beta distributions $\mathrm{Beta}(1.5,3.5)$, $\mathrm{Beta}(6,14)$, and $\mathrm{Beta}(24,56)$. Their variances are $0.035$, $0.01$, and $0.21/81\approx 0.00259259$, respectively. Parameters need not be integers.

The formulas follow from two short moment calculations. Multiplying the probability density function by $x$ or $x^2$ increases the first parameter of the Euler beta function:Apply the Gamma function recurrence to the beta-gamma identity, then subtract $(\mathbb{E}[X]{)}^2$ from $\mathbb{E}[X^2]$ to obtain the variance.

The beta distribution is especially convenient as a prior distribution for the unknown success probability $P$ of Bernoulli trials. The Beta-binomial conjugacy says that $h$ successes and $t$ failures update a prior distribution $\mathrm{Beta}(\alpha ,\beta )$ to a posterior distribution $\mathrm{Beta}(\alpha +h,\beta +t)$. For example, suppose a coin has prior distribution $\mathrm{Beta}(2,2)$ and we observe eight heads and two tails. Its posterior distribution is then $\mathrm{Beta}(10,4)$, with posterior mean $5/7$. ^[ref].

A truncated distribution is a conditional probability distribution obtained by restricting a random variable to a specified subset of its original support and renormalizing.

For a wavefunction satisfying the Schrödinger equation with a real potential,This continuity equation expresses local conservation of probability.

The Born rule connects a quantum state to measurement probabilities. Let $\{|j⟩\}$ be a complete orthonormal basis, finite or countably infinite, for a Hilbert space, and letAn ideal measurement distinguishing these basis vectors returns label $j$ with probabilityThis assignment is a physical measurement rule. Its compatibility with total probability one follows from orthonormality and normalization:The probabilities are the squared moduli, not the amplitudes themselves or the squares $a_j^2$. ^[ref].

If the basis is an eigenbasis of an observable, the reported numerical value is its eigenvalue. Several vectors can share one eigenvalue; the probability of that value is then the sum over the corresponding basis labels. More generally, a projective measurement with orthogonal projections $P_r$ assignsfor a pure state and a general density matrix $\rho$, respectively. This also covers mixed states and degenerate eigenspaces. The Born rule determines outcome statistics; the ideal conditional state change is specified separately by the Lüders rule. See projective measurement for that formulation and its assumptions. ^[ref].

For a worked example, take the normalized qubit quantum stateThe squared moduli add to $3/4+1/4=1$. A computational-basis measurement, equivalently a measurement of the $Z$ Pauli matrix, hasThese probabilities do not depend on the relative quantum phase $\varphi$. Now choose the Hadamard basiswhose vectors are eigenstates of the $X$ Pauli matrix, also used by the Pauli X gate. The new probability amplitudes and probabilities areBoth probabilities are nonnegative and their sum is one. The relative quantum phase becomes visible because the two probability amplitudes add before their squared modulus is taken:

A Hadamard gate followed by a computational-basis measurement realizes these $X$ statistics, since $H|\pm ⟩=|0⟩,|1⟩$, respectively. Compare the two bases using fresh, identically prepared copies: a preceding measurement may change the state. The quantum superposition section shows why replacing the preparation by the mixed state $\mathrm{diag}(3/4,1/4)$ keeps the $Z$ statistics but gives $p_{+}=p_{-}=1/2$ for every $\varphi$.

The usefulness of the Born rule is this conversion from linear algebra to experimental predictions: calculate overlaps in the basis actually measured, then take squared moduli. Its appealing geometry is that changing basis redistributes probability while preserving total squared norm. The same rule relates an energy eigenstate expansion to an energy measurement and a position wavefunction to a probability density.

Born's 1926 paper on quantum scattering interpreted the outgoing-wave coefficients in terms of the probabilities of collision outcomes; a corrective footnote specified proportionality to the square of the amplitude. This is a historical scattering formulation, rather than the general density matrix formula above. ^[ref]; ^[ref].

If a graph $H$ has minimum degree $d$, then its graph Ramsey number satisfies $R(H)\ge 2^{d/2}$. The probabilistic proof chooses a red-blue edge colouring and applies the local lemma to its monochromatic copies of $H$.

For every $\epsilon >0$ there is $\delta >0$ such that a graph satisfying $e(A,B)\ge p|A||B|$ whenever $|A|,|B|\ge \delta n$, with $p\ge \delta$, has a spanning subgraph $G^{\prime }$ satisfyingfor all disjoint $A,B$. Partition into a fixed large number of nearly equal parts, independently thin each cross-pair to expected density $p$, and apply a concentration inequality and a union bound over all pairs $A,B$.

For disjoint nonempty vertex sets $A,B$, their edge density is

Dependent random choice finds a large vertex set whose small subsets have large common neighbourhoods. In a bipartite graph with parts $A,B$, choose random vertices of $B$ with repetition and take their common neighbourhood $X\subseteq A$; convexity gives a lower bound for $\mathbb{E}|X|$, while counting subsets of $X$ with small common neighbourhood allows their deletion.

For bad events with a dependency graph of maximum degree $D$, the symmetric Lovász local lemma guarantees positive probability that none occurs whenever each event has probability at most $p$ and $ep(D+1)\le 1$.

The random alteration method samples a random object and then deletes a small number of offending elements. If few forbidden configurations occur while the desired global property already holds, the altered object witnesses existence.

For a decreasing probability sequence $p(n)$, if $G(n,p(n))$ contains a Hamilton cycle with high probability, then $G(n,Cp(n))$ is pancyclic with high probability for every fixed sufficiently large $C$; three independent rounds suffice. One first exposes a Hamiltonian round and uses the independent sprinkled edges to create cycles of every shorter length.

For every fixed $\epsilon >0$, every component of $G(n,(1-\epsilon )/n)$ has $O_{\epsilon }(\log n)$ vertices with high probability. A breadth-first exploration is dominated by a branching process of mean $1-\epsilon$; its total progeny has an exponentially decreasing tail.

An event depending on $n$ holds with high probability when its probability tends to one as $n\to \infty$.

If $e_i$ is fixed in a rigid body with angular velocity $\omega$, then its derivative in the fixed space frame is

Statistical homogeneity defines the primordial bispectrum by

A bispectrum is a Fourier-space third-order statistic; primordial perturbations supply the cosmological field in this instance.