Bernoulli process explained

In probability and statistics, a Bernoulli process (named after Jacob Bernoulli) is a finite or infinite sequence of binary random variables, so it is a discrete-time stochastic process that takes only two values, canonically 0 and 1. The component Bernoulli variables Xi are identically distributed and independent. Prosaically, a Bernoulli process is a repeated coin flipping, possibly with an unfair coin (but with consistent unfairness). Every variable Xi in the sequence is associated with a Bernoulli trial or experiment. They all have the same Bernoulli distribution. Much of what can be said about the Bernoulli process can also be generalized to more than two outcomes (such as the process for a six-sided die); this generalization is known as the Bernoulli scheme.

The problem of determining the process, given only a limited sample of Bernoulli trials, may be called the problem of checking whether a coin is fair.

Definition

A Bernoulli process is a finite or infinite sequence of independent random variables X1X2X3, ..., such that

In other words, a Bernoulli process is a sequence of independent identically distributed Bernoulli trials.

Independence of the trials implies that the process is memoryless. Given that the probability p is known, past outcomes provide no information about future outcomes. (If p is unknown, however, the past informs about the future indirectly, through inferences about p.)

If the process is infinite, then from any point the future trials constitute a Bernoulli process identical to the whole process, the fresh-start property.

Interpretation

The two possible values of each Xi are often called "success" and "failure". Thus, when expressed as a number 0 or 1, the outcome may be called the number of successes on the ith "trial".

Two other common interpretations of the values are true or false and yes or no. Under any interpretation of the two values, the individual variables Xi may be called Bernoulli trials with parameter p.

In many applications time passes between trials, as the index i increases. In effect, the trials X1X2, ... Xi, ... happen at "points in time" 1, 2, ..., i, .... That passage of time and the associated notions of "past" and "future" are not necessary, however. Most generally, any Xi and Xj in the process are simply two from a set of random variables indexed by, the finite cases, or by, the infinite cases.

One experiment with only two possible outcomes, often referred to as "success" and "failure", usually encoded as 1 and 0, can be modeled as a Bernoulli distribution.[1] Several random variables and probability distributions beside the Bernoullis may be derived from the Bernoulli process:

The negative binomial variables may be interpreted as random waiting times.

Formal definition

The Bernoulli process can be formalized in the language of probability spaces as a random sequence of independent realisations of a random variable that can take values of heads or tails. The state space for an individual value is denoted by

2=\{H,T\}.

Borel algebra

Consider the countably infinite direct product of copies of

2=\{H,T\}

. It is common to examine either the one-sided set

\Omega=2N=\{H,T\}N

or the two-sided set

\Omega=2Z

. There is a natural topology on this space, called the product topology. The sets in this topology are finite sequences of coin flips, that is, finite-length strings of H and T (H stands for heads and T stands for tails), with the rest of (infinitely long) sequence taken as "don't care". These sets of finite sequences are referred to as cylinder sets in the product topology. The set of all such strings forms a sigma algebra, specifically, a Borel algebra. This algebra is then commonly written as

(\Omega,l{B})

where the elements of

l{B}

are the finite-length sequences of coin flips (the cylinder sets).

Bernoulli measure

If the chances of flipping heads or tails are given by the probabilities

\{p,1-p\}

, then one can define a natural measure on the product space, given by

P=\{p,1-p\}N

(or by

P=\{p,1-p\}Z

for the two-sided process). In another word, if a discrete random variable X has a Bernoulli distribution with parameter p, where 0 ≤ p ≤ 1, and its probability mass function is given by

pX(1)=P(X=1)=p

and

pX(0)=P(X=0)=1-p

.

We denote this distribution by Ber(p).

Given a cylinder set, that is, a specific sequence of coin flip results

[\omega1,\omega2,\omegan]

at times

1,2,,n

, the probability of observing this particular sequence is given by

P([\omega1,\omega2,,\omegan])=pk(1-p)n-k

where k is the number of times that H appears in the sequence, and nk is the number of times that T appears in the sequence. There are several different kinds of notations for the above; a common one is to write

P(X1=x1,X2=x2,,Xn=xn)=pk(1-p)n-k

where each

Xi

is a binary-valued random variable with

xi=[\omegai=H]

in Iverson bracket notation, meaning either

1

if

\omegai=H

or

0

if

\omegai=T

. This probability

P

is commonly called the Bernoulli measure.[2]

Note that the probability of any specific, infinitely long sequence of coin flips is exactly zero; this is because

\limn\toinftypn=0

, for any

0\lep<1

. A probability equal to 1 implies that any given infinite sequence has measure zero. Nevertheless, one can still say that some classes of infinite sequences of coin flips are far more likely than others, this is given by the asymptotic equipartition property.

To conclude the formal definition, a Bernoulli process is then given by the probability triple

(\Omega,l{B},P)

, as defined above.

Law of large numbers, binomial distribution and central limit theorem

See main article: article, Law of large numbers, Central limit theorem and Binomial distribution. Let us assume the canonical process with

H

represented by

1

and

T

represented by

0

. The law of large numbers states that the average of the sequence, i.e.,

\bar{X}n:=

1
n
n
\sum
i=1

Xi

, will approach the expected value almost certainly, that is, the events which do not satisfy this limit have zero probability. The expectation value of flipping heads, assumed to be represented by 1, is given by

p

. In fact, one has

E[Xi]=P([Xi=1])=p,

for any given random variable

Xi

out of the infinite sequence of Bernoulli trials that compose the Bernoulli process.

N(k,n)={n\choosek}=

n!
k!(n-k)!

If the probability of flipping heads is given by p, then the total probability of seeing a string of length n with k heads is

P([Sn=k])={n\choosek}pk(1-p)n-k,

where

Sn=\sum

n
i=1

Xi

.The probability measure thus defined is known as the Binomial distribution.

As we can see from the above formula that, if n=1, the Binomial distribution will turn into a Bernoulli distribution. So we can know that the Bernoulli distribution is exactly a special case of Binomial distribution when n equals to 1.

Of particular interest is the question of the value of

Sn

for a sufficiently long sequences of coin flips, that is, for the limit

n\toinfty

. In this case, one may make use of Stirling's approximation to the factorial, and write

n!=\sqrt{2\pin}nne-n\left(1+

l{O}\left(1
n
\right)\right)

Inserting this into the expression for P(k,n), one obtains the Normal distribution; this is the content of the central limit theorem, and this is the simplest example thereof.

The combination of the law of large numbers, together with the central limit theorem, leads to an interesting and perhaps surprising result: the asymptotic equipartition property. Put informally, one notes that, yes, over many coin flips, one will observe H exactly p fraction of the time, and that this corresponds exactly with the peak of the Gaussian. The asymptotic equipartition property essentially states that this peak is infinitely sharp, with infinite fall-off on either side. That is, given the set of all possible infinitely long strings of H and T occurring in the Bernoulli process, this set is partitioned into two: those strings that occur with probability 1, and those that occur with probability 0. This partitioning is known as the Kolmogorov 0-1 law.

The size of this set is interesting, also, and can be explicitly determined: the logarithm of it is exactly the entropy of the Bernoulli process. Once again, consider the set of all strings of length n. The size of this set is

2n

. Of these, only a certain subset are likely; the size of this set is

2nH

for

H\le1

. By using Stirling's approximation, putting it into the expression for P(k,n), solving for the location and width of the peak, and finally taking

n\toinfty

one finds that

H=-plog2p-(1-p)log2(1-p)

This value is the Bernoulli entropy of a Bernoulli process. Here, H stands for entropy; not to be confused with the same symbol H standing for heads.

John von Neumann posed a question about the Bernoulli process regarding the possibility of a given process being isomorphic to another, in the sense of the isomorphism of dynamical systems. The question long defied analysis, but was finally and completely answered with the Ornstein isomorphism theorem. This breakthrough resulted in the understanding that the Bernoulli process is unique and universal; in a certain sense, it is the single most random process possible; nothing is 'more' random than the Bernoulli process (although one must be careful with this informal statement; certainly, systems that are mixing are, in a certain sense, "stronger" than the Bernoulli process, which is merely ergodic but not mixing. However, such processes do not consist of independent random variables: indeed, many purely deterministic, non-random systems can be mixing).

Dynamical systems

The Bernoulli process can also be understood to be a dynamical system, as an example of an ergodic system and specifically, a measure-preserving dynamical system, in one of several different ways. One way is as a shift space, and the other is as an odometer. These are reviewed below.

Bernoulli shift

See main article: article, Bernoulli scheme and Dyadic transformation. One way to create a dynamical system out of the Bernoulli process is as a shift space. There is a natural translation symmetry on the product space

\Omega=2N

given by the shift operator

T(X0,X1,X2,)=(X1,X2,)

The Bernoulli measure, defined above, is translation-invariant; that is, given any cylinder set

\sigma\inl{B}

, one has

P(T-1(\sigma))=P(\sigma)

and thus the Bernoulli measure is a Haar measure; it is an invariant measure on the product space.

Instead of the probability measure

f\circT-1

defined by

\left(f\circT-1\right)(\sigma)=f(T-1(\sigma))

is again some function

l{B}\toR.

Thus, the map

T

induces another map

l{L}T

on the space of all functions

l{B}\toR.

That is, given some

f:l{B}\toR

, one defines

l{L}Tf=f\circT-1

The map

l{L}T

is a linear operator, as (obviously) one has

l{L}T(f+g)=l{L}T(f)+l{L}T(g)

and

l{L}T(af)=al{L}T(f)

for functions

f,g

and constant

a

. This linear operator is called the transfer operator or the Ruelle–Frobenius–Perron operator. This operator has a spectrum, that is, a collection of eigenfunctions and corresponding eigenvalues. The largest eigenvalue is the Frobenius–Perron eigenvalue, and in this case, it is 1. The associated eigenvector is the invariant measure: in this case, it is the Bernoulli measure. That is,

l{L}T(P)=P.

If one restricts

l{L}T

to act on polynomials, then the eigenfunctions are (curiously) the Bernoulli polynomials![3] [4] This coincidence of naming was presumably not known to Bernoulli.

The 2x mod 1 map

Notes and References

  1. Book: A modern introduction to probability and statistics. 9781852338961. 45–46. Dekking. F. M.. Kraaikamp. C.. Lopuhaä. H. P.. Meester. L. E.. 2005.
  2. Book: Klenke, Achim . Probability Theory . 2006 . Springer-Verlag . 978-1-84800-047-6.
  3. Pierre Gaspard, "r-adic one-dimensional maps and the Euler summation formula", Journal of Physics A, 25 (letter) L483-L485 (1992).
  4. Dean J. Driebe, Fully Chaotic Maps and Broken Time Symmetry, (1999) Kluwer Academic Publishers, Dordrecht Netherlands