Skip to main content

Distributions

Lots to do here

This is incomplete and missing a most important thing: Distribution ‘Arithmetic’ (e.g. what is the sum of two Uniform distributions?)

Think of these as a class in Python or some other OOP language. If you identify that an Random Variable ‘fits’ one of these, you can use precalculated/derived formulae for things like the Expected Value E(X)E(X) and Variance Var(X)Var(X) (for baby Statisticians like you and I.)

Don't Overthink it, Sunshine

A ‘distribution’ doesn’t have to be anything fancy and can just be a small list that tells people what the probabilities are for each value (“realization”) of your Random Variable. E.g. your shiny Random Variable XX can only take on values {1,2,3}\{1, 2, 3\} and you assign each value (“realization”) probabilities of {29,49,39}\{\frac{2}{9}, \frac{4}{9}, \frac{3}{9}\}.

P(X=1)=29P(X = 1) = \frac{2}{9}, P(X=2)=49P(X = 2) = \frac{4}{9}, P(X=3)=39P(X = 3) = \frac{3}{9}.

Bam, done, Distribution.

This applies to joint distributions as well! Let’s say XX can take on {1,2}\{1,2\} and YY can take on {a,b}\{a, b\}. You’d do a nice table of P(X=1,Y=a)P(X = 1, Y = a), P(X=2,Y=a)P(X = 2, Y = a), P(X=1,Y=b)P(X = 1, Y = b), and P(X=2,Y=b)P(X = 2, Y = b). Note that I’m not saying anything about independence or anything. Just a table of values.

All your probability values, joint or not, must add up to 11.

There are other rules of course. But that’s all really. A function is nice because you won’t have to labor over this list for a lot of values (and if you’re laboring over R\Reals, you’ll be laboring for a long time indeed.)

Important: This applies to those Contingency Tables as well! Think about it!

The Normal Distribution ♥️​

TODO: Write more about this and the CLT.

Here’s a lovely derivation that’s worth watching.

Discrete Distributions​

Bernoulli​

Simplest one. If an experiment with probability pp succeeds, you get a 1, else 0.

X∼Ber(p)X \sim Ber(p)
PMFE[X]E[X]Var(X)Var(X)
P(X=1)=pP(X = 1) = p, P(X=0)=1−pP(X = 0) = 1-p ppp(1−p)p(1-p)

Examples: Coin flip results in a heads. Your friend liked Nights in Rodanthe.

Binomial​

Extends Bernoulli to nn independent trials of your experiment with probability of success pp

X∼Bin(n,p)X \sim Bin(n,p)
PMFE[X]E[X]Var(X)Var(X)
P(X=k)=(nk)pk(1−p)n−kP(X = k) = {n \choose k} p^k(1-p)^{n-k}npnpnp(1−p)np(1-p)

Examples: Number of heads in 12 coin flips. Bits being corrupted (think Hamming Codes).

Poisson​

X∼Poi(λ)X \sim Poi(\lambda)

Approximates Binomial to cases when probablity pp is very small and the number of trials nn is very large. It’s pretty important and used in queueing theory.

  • Think the number of things that happen over a fixed interval of time at a constant average rate.
  • Think of events that are rare, independent (of course), and relatively uniformly distributed in time or space: mutations in a stretch of DNA, calls into a call center, insurance claims filed throughout a day, patients arriving at a clinic, rare diseases reported per month, meteors observed in the night sky, bankruptcies observed in a year…

Geometric​

X∼Geo(p)X \sim Geo(p)

Number of independent trials until first success (which has a probability pp).

Negative Binomial​

X∼NegBin(p)X \sim NegBin(p)

This is the number of independent trials you’ll need until a fixed number of successes rr with probability pp. Kinda extends the idea of the Geometric.