Tutorials › Stochastic Differential Equations in Pharmacokinetics › Defining Standard Brownian Motion

Brownian Motion and the Wiener Process · Tutorial 14 of 100

Defining Standard Brownian Motion

Learn the defining assumptions of standard Brownian motion, derive its basic distributions, and distinguish a particular random process from the Wiener law shared by all standard Brownian motions.

Intermediate 10 min read

What You'll Learn

  • State the initial-value, increment, and path-continuity conditions defining standard Brownian motion.
  • Calculate increment distributions and probabilities from elapsed time.
  • Derive the mean, variance, and covariance of Brownian values.
  • Explain how a Brownian process differs from its Wiener law on path space.
  • Identify data patterns that would be inconsistent with the defining assumptions.

Why Define the Noise Precisely?

In a pharmacokinetic model, a stochastic driver represents variation in a state between observation times, not merely scatter in the measured concentrations. To specify a Brownian driver, it is not enough to say that it is “random” or “normally distributed.” Its starting value, the distribution and dependence of its increments, and the regularity of its paths are all part of the definition. These conditions determine which probability calculations are justified.

The term standard Brownian motion describes a stochastic process with a particular law and path behavior. A specific process is defined on a probability space and has particular sample paths; its law describes the distribution of those paths. Keeping these two meanings separate is useful when comparing models, simulations, or mathematical constructions.

Definition of Standard Brownian Motion

Definition: A real-valued process \(W=(W_t)_{t\geq 0}\), defined on a probability space, is a standard Brownian motion if:
  • \(W_0=0\) almost surely;
  • for every finite collection of disjoint time intervals, the increments over those intervals are mutually independent;
  • for every \(0\leq s<t\), \(W_t-W_s\sim N(0,t-s)\); and
  • with probability one, the function \(t\mapsto W_t\) is continuous on \([0,\infty)\).

The variance parameter in \(N(0,t-s)\) is the elapsed time, not its square. The independent-increment condition means that increments on disjoint intervals are independent random variables, even if the intervals are not adjacent. Equivalently, each increment after time \(s\) is independent of the process history up to \(s\), represented by \(\sigma(W_u:0\leq u\leq s)\).

The continuity condition is almost sure: there may be a probability-zero set of outcomes on which the path is not continuous. The definition does not assert that every possible outcome has a continuous path. Nor does it say that a path is deterministic or differentiable. Continuity concerns the realized function of time; the Gaussian and independence conditions concern probability distributions.

Worked Example: Distribution of an Increment

Let \(W\) be standard Brownian motion. Find the distribution of \(W_2-W_{0.5}\), and calculate the probability that this increment exceeds 1.

The interval length is \(2-0.5=1.5\). By the definition, \[ W_2-W_{0.5}\sim N(0,1.5). \] To standardize the threshold, divide by the standard deviation \(\sqrt{1.5}\): \[ P(W_2-W_{0.5}>1) = P\left(Z>\frac{1}{\sqrt{1.5}}\right) = P(Z>0.8165), \] where \(Z\sim N(0,1)\). Using the standard normal distribution function, \(\Phi(0.8165)\approx 0.7929\), so the probability is approximately \(1-0.7929=0.2071\).

The result depends on the time difference, 1.5, rather than on the absolute location of the interval. An increment over an interval of the same length has the same distribution.

Consequences for Means and Covariances

The definition gives the distribution of increments directly. It also determines the first and second moments of the process at individual times. The following result is useful when checking calculations involving Brownian values.

Theorem: If \(W\) is standard Brownian motion, then for all \(s,t\geq 0\), \[ E[W_t]=0,\qquad \operatorname{Var}(W_t)=t,\qquad \operatorname{Cov}(W_s,W_t)=\min(s,t). \]

Proof. Since \(W_t-W_0\sim N(0,t)\) and \(W_0=0\) almost surely, \(W_t\sim N(0,t)\). Therefore \(E[W_t]=0\) and \(\operatorname{Var}(W_t)=t\). To obtain the covariance, first suppose \(s\leq t\) and write \[ W_t=W_s+(W_t-W_s). \] The increment \(W_t-W_s\) is independent of \(W_s\), has mean zero, and has finite variance. Thus \[ \operatorname{Cov}(W_s,W_t) =\operatorname{Cov}(W_s,W_s)+\operatorname{Cov}(W_s,W_t-W_s) =s+0=s. \] If \(t<s\), the same argument with the times reversed gives covariance \(t\). In either case the covariance is \(\min(s,t)\). \(\square\)

In particular, Brownian values at two different times are generally not independent: their covariance is positive when both times are positive. It is increments on disjoint intervals that are independent. For example, \(W_s\) and \(W_t-W_s\) are independent for \(0\leq s<t\), while \(W_s\) and \(W_t\) are not independent when \(s>0\).

Worked Example: Covariance and an Independent Increment

For standard Brownian motion, calculate \(\operatorname{Cov}(W_2,W_5)\), \(\operatorname{Var}(W_5-W_2)\), and the correlation between \(W_2\) and \(W_5\).

The covariance theorem gives \[ \operatorname{Cov}(W_2,W_5)=\min(2,5)=2. \] The increment \(W_5-W_2\) has distribution \(N(0,5-2)=N(0,3)\), so its variance is 3. The correlation is \[ \operatorname{Corr}(W_2,W_5) =\frac{\operatorname{Cov}(W_2,W_5)} {\sqrt{\operatorname{Var}(W_2)\operatorname{Var}(W_5)}} =\frac{2}{\sqrt{2\cdot 5}} =\sqrt{0.4}\approx 0.6325. \] Although \(W_2\) and \(W_5\) are correlated, \(W_2\) and \(W_5-W_2\) are independent by the independent-increment condition. Their covariance is zero, as can also be checked by subtracting: \[ \operatorname{Cov}(W_2,W_5-W_2) =\operatorname{Cov}(W_2,W_5)-\operatorname{Var}(W_2) =2-2=0. \]

A Process and Its Wiener Law

A process \(W\) is a collection of random variables on a probability space. At a particular outcome \(\omega\), it produces a path \(t\mapsto W_t(\omega)\). The law of the process describes the probabilities assigned to possible paths; it does not identify which path occurred on a particular outcome.

Definition: Let \(C_0[0,\infty)\) be the space of continuous real-valued functions \(x\) on \([0,\infty)\) satisfying \(x(0)=0\), equipped with its Borel sets for the topology of uniform convergence on bounded time intervals. The Wiener measure is the probability law on this path space induced by a standard Brownian motion. Thus it describes the distribution of a Brownian path, not a particular realized path.

The defining conditions determine a unique Wiener measure. The key reason is that the increment distributions and independence determine every finite-dimensional distribution, while continuous paths are determined by their values at rational times. Here is the precise argument.

Theorem: Any two standard Brownian motions induce the same probability law on \(C_0[0,\infty)\).

Proof. Choose any finite set of times \(0\leq t_1<\cdots<t_m\). The vector of Brownian values \((W_{t_1},\ldots,W_{t_m})\) is obtained by adding the increments over the successive intervals \([0,t_1], [t_1,t_2],\ldots,[t_{m-1},t_m]\). Those increments are independent, and each has a normal distribution with variance equal to its interval length. Consequently, the joint distribution of the vector is determined by the definition and is the same for every standard Brownian motion.

It remains to pass from finite-dimensional distributions to measures on path space. Evaluations \(x\mapsto x(q)\) at nonnegative rational times \(q\) generate the Borel sigma-field of \(C_0[0,\infty)\). To see why continuity matters here, for any continuous path the supremum of its distance from another continuous path over a bounded interval equals the supremum over rational times in that interval. Thus open sets for uniform convergence on bounded intervals can be described using countably many rational-time evaluations. The two path laws agree on every finite collection of these evaluations by the preceding argument. The corresponding cylinder sets form a generating class closed under finite intersections, so uniqueness of probability measures implies that the laws agree on the Borel sigma-field. \(\square\)

Worked Example: Same Law, Different Process Values

Suppose \(W\) is standard Brownian motion and define \(\widetilde W_t=-W_t\) for every \(t\geq 0\). Verify that \(\widetilde W\) is standard Brownian motion, and explain how it can differ from \(W\).

First, \(\widetilde W_0=-W_0=0\) almost surely, and its paths are continuous whenever the paths of \(W\) are continuous. For \(s<t\), \[ \widetilde W_t-\widetilde W_s=-(W_t-W_s)\sim N(0,t-s), \] because changing the sign of a zero-mean normal random variable does not change its distribution. Negating each increment also preserves independence, so \(\widetilde W\) satisfies all the defining conditions.

Therefore \(W\) and \(\widetilde W\) have the same Wiener law. But they are not equal at time 1 except on the event \(W_1=0\). Since \(W_1\sim N(0,1)\), \(P(W_1=0)=0\). Thus their values at time 1 differ with probability one, even though the two processes have the same law. Equality of laws means equality of distributions, not equality of sample paths.

Model Checks, Alternatives, and Limitations

When Brownian motion is used as a model assumption, data can be checked against its implications, but finite observations cannot establish the full definition. For observed increments over intervals of known lengths, useful diagnostics include:

  • Mean and variance: increments should be centered, with variance proportional to elapsed time.
  • Distribution: standardized increments should be broadly compatible with a normal distribution, allowing for sampling uncertainty.
  • Dependence: increments over non-overlapping intervals should not show systematic dependence.
  • Path behavior: obvious jumps are inconsistent with the almost-sure continuity assumption, although sparse measurements cannot verify continuity between sampling times.

These checks are especially important when observations are sparse: a smooth-looking sequence of measured concentrations does not demonstrate that the latent process has continuous paths, and irregular measurements do not by themselves identify which source of variability produced them. A model’s observation process and its latent process must be distinguished when interpreting such patterns.

Standard Brownian motion is a convenient model for continuous random variation with independent, time-homogeneous Gaussian increments. It is not a universal description of randomness. A deterministic smooth input has no random increments; a jump process allows discontinuities; and a mean-reverting Gaussian process generally has dependent increments. These alternatives may be more appropriate when the scientific mechanism or observed behavior contradicts Brownian assumptions. Choosing among them is a modeling decision, not something resolved by calling all unexplained variability “noise.”

A common mistake is to infer independence of \(W_s\) and \(W_t\) from the fact that Brownian motion has independent increments. The definition applies to increments on disjoint intervals, not to the process values themselves. For \(0<s<t\), the values share the earlier part of the path, and their covariance is \(s\), not zero. The correction is to express the quantities of interest as increments first, then check whether their time intervals are disjoint.

Key takeaway: Standard Brownian motion starts at zero, has independent Gaussian increments with variance equal to elapsed time, and has continuous paths almost surely. Its Wiener law specifies the distribution of complete paths; it is not one particular path or process realization.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. What is the distribution of \(W_4-W_1\) for standard Brownian motion?
  2. Are \(W_1\) and \(W_4-W_1\) independent? What about \(W_1\) and \(W_4\)?
  3. State the covariance of \(W_3\) and \(W_7\).
  4. What does it mean for two Brownian motions to have the same Wiener law? Does it require their paths to match on the same probability space?

Solutions. (1) The interval length is \(4-1=3\), so \(W_4-W_1\sim N(0,3)\). (2) \(W_1\) and \(W_4-W_1\) are independent because they concern disjoint intervals, \([0,1]\) and \([1,4]\). In contrast, \(\operatorname{Cov}(W_1,W_4)=\min(1,4)=1\), so \(W_1\) and \(W_4\) are not independent. (3) \(\operatorname{Cov}(W_3,W_7)=\min(3,7)=3\). (4) The two processes assign the same probabilities to Borel sets of continuous paths. Their realized paths need not match; for example, a Brownian motion and its negative have the same law but differ at time 1 with probability one.

Clinical Biostats