Power Law Size Distributions Overview Introduction Principles of - - PowerPoint PPT Presentation

power law size distributions
SMART_READER_LITE
LIVE PREVIEW

Power Law Size Distributions Overview Introduction Principles of - - PowerPoint PPT Presentation

Power Law Size Distributions Power Law Size Distributions Overview Introduction Principles of Complex Systems Examples Zipfs law Course 300, Fall, 2008 Wild vs. Mild CCDFs References Prof. Peter Dodds Department of Mathematics &


slide-1
SLIDE 1

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 1/33

Power Law Size Distributions

Principles of Complex Systems Course 300, Fall, 2008

  • Prof. Peter Dodds

Department of Mathematics & Statistics University of Vermont

Licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License.

slide-2
SLIDE 2

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 2/33

Outline

Overview Introduction Examples Zipf’s law Wild vs. Mild CCDFs References

slide-3
SLIDE 3

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 3/33

The Don

Extreme deviations in test cricket

100 10 20 30 90 40 50 60 70 80

Don Bradman’s batting average = 166% next best.

slide-4
SLIDE 4

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 5/33

Size distributions

The sizes of many systems’ elements appear to obey an inverse power-law size distribution: P(size = x) ∼ c x−γ where xmin < x < xmax and γ > 1

◮ Typically, 2 < γ < 3. ◮ xmin = lower cutoff ◮ xmax = upper cutoff

slide-5
SLIDE 5

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 6/33

Size distributions

◮ Usually, only the tail of the distribution obeys a power

law: P(x) ∼ c x−γ as x → ∞.

◮ Still use term ‘power law distribution’

slide-6
SLIDE 6

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 7/33

Size distributions

Many systems have discrete sizes k:

◮ Word frequency ◮ Node degree (as we have seen): # hyperlinks, etc. ◮ number of citations for articles, court decisions, etc.

P(k) ∼ c k−γ where kmin ≤ k ≤ kmax

slide-7
SLIDE 7

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 8/33

Size distributions

Power law size distributions are sometimes called Pareto distributions after Italian scholar Vilfredo Pareto.

◮ Pareto noted wealth in Italy was distributed unevenly

(80–20 rule).

◮ Term used especially by economists

slide-8
SLIDE 8

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 9/33

Size distributions

◮ Negative linear relationship in log-log space:

log P(x) = log c − γ log x

slide-9
SLIDE 9

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 11/33

Size distributions

Examples:

◮ Earthquake magnitude (Gutenberg Richter law):

P(M) ∝ M−3

◮ Number of war deaths: P(d) ∝ d−1.8 ◮ Sizes of forest fires ◮ Sizes of cities: P(n) ∝ n−2.1 ◮ Number of links to and from websites

slide-10
SLIDE 10

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 12/33

Size distributions

Examples:

◮ Number of citations to papers: P(k) ∝ k−3. ◮ Individual wealth (maybe): P(W) ∝ W −2. ◮ Distributions of tree trunk diameters: P(d) ∝ d−2. ◮ The gravitational force at a random point in the

universe: P(F) ∝ F −5/2.

◮ Diameter of moon craters: P(d) ∝ d−3. ◮ Word frequency: e.g., P(k) ∝ k−2.2 (variable)

(Note: Exponents range in error; see M.E.J. Newman arxiv.org/cond-mat/0412004v3 (⊞))

slide-11
SLIDE 11

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 13/33

Size distributions

Power-law distributions are..

◮ often called ‘heavy-tailed’ ◮ or said to have ‘fat tails’

Important!:

◮ Inverse power laws aren’t the only ones:

◮ lognormals, stretched exponentials, ...

slide-12
SLIDE 12

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 15/33

Zipfian rank-frequency plots

George Kingsley Zipf:

◮ noted various rank distributions

followed power laws, often with exponent -1 (word frequency, city sizes...) “Human Behaviour and the Principle of Least-Effort” [2] Addison-Wesley,

Cambridge MA, 1949.

◮ We’ll study Zipf’s law in depth...

slide-13
SLIDE 13

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 16/33

Zipfian rank-frequency plots

Zipf’s way:

◮ si = the size of the ith ranked object. ◮ i = 1 corresponds to the largest size. ◮ s1 could be the frequency of occurrence of the most

common word in a text.

◮ Zipf’s observation:

si ∝ i−α

slide-14
SLIDE 14

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 18/33

Power law distributions

Gaussians versus power-law distributions:

◮ Example: Height versus wealth. ◮ Mild versus Wild (Mandelbrot) ◮ Mediocristan versus Extremistan

(See “The Black Swan” by Nassim Taleb [1])

slide-15
SLIDE 15

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 19/33

Turkeys...

From “The Black Swan” [1]

slide-16
SLIDE 16

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 20/33

Taleb’s table [1]

Mediocristan/Extremistan

◮ Most typical member is mediocre/Most typical is either

giant or tiny

◮ Winners get a small segment/Winner take almost all

effects

◮ When you observe for a while, you know what’s going on/

It takes a very long time to figure out what’s going on

◮ Prediction is easy/Prediction is hard ◮ History crawls/History makes jumps ◮ Tyranny of the collective/Tyranny of the accidental

slide-17
SLIDE 17

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 22/33

Complementary Cumulative Distribution Function:

CCDF:

P≥(x) = P(x′ ≥ x) = 1 − P(x′ < x)

= ∞

x′=x

P(x′)dx′

∝ ∞

x′=x

(x′)−γdx′

= 1 −γ + 1(x′)−γ+1

x′=x ◮

∝ x−γ+1

slide-18
SLIDE 18

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 23/33

Complementary Cumulative Distribution Function:

CCDF:

P≥(x) ∝ x−γ+1

◮ Use when tail of P follows a power law. ◮ Increases exponent by one. ◮ Useful in cleaning up data.

slide-19
SLIDE 19

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 24/33

Complementary Cumulative Distribution Function:

◮ Discrete variables:

P≥(k) = P(k′ ≥ k) =

  • k′=k

P(k) ∝ k−γ+1

◮ Use integrals to approximate sums.

slide-20
SLIDE 20

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 25/33

Size distributions

Brown Corpus (1,015,945 words):

CCDF:

−2.5 −2 −1.5 −1 −0.5 0.5 1 0.5 1 1.5 2 2.5 3 3.5

n N> n

Zipf:

0.5 1 1.5 2 2.5 3 3.5 −2.5 −2 −1.5 −1 −0.5 0.5 1

rank i ni

◮ The, of, and, to, a, ... = ‘objects’ ◮ ‘Size’ = word frequency ◮ Beep: CCDF and Zipf plots are related...

slide-21
SLIDE 21

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 26/33

Size distributions

Observe:

◮ NP≥(x) = the number of objects with size at least x

where N = total number of objects.

◮ If an object has size xi, then NP≥(xi) is its rank i. ◮ So

xi ∝ i−α = (NP≥(xi))−α ∝ x(−γ+1)(−α)

i

Since P≥(x) ∼ x−γ+1, α = 1 γ − 1 A rank distribution exponent of α = 1 corresponds to a size distribution exponent γ = 2.

slide-22
SLIDE 22

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 27/33

Details on the lack of scale:

Let’s find the mean:

x = xmax

x=xmin

xP(x)dx = c xmax

x=xmin

xx−γdx = c 2 − γ

  • x2−γ

max − x2−γ min

  • .
slide-23
SLIDE 23

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 28/33

The mean:

x ∼ c 2 − γ

  • x2−γ

max − x2−γ min

  • .

◮ Mean blows up with upper cutoff if γ < 2. ◮ Mean depends on lower cutoff if γ > 2. ◮ γ < 2: Typical sample is large. ◮ γ > 2: Typical sample is small.

slide-24
SLIDE 24

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 29/33

And in general...

Moments:

◮ All moments depend only on cutoffs. ◮ No internal scale dominates (even matters). ◮ Compare to a Gaussian, exponential, etc.

slide-25
SLIDE 25

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 30/33

Moments

For many real size distributions:

2 < γ < 3

◮ mean is finite (depends on lower cutoff) ◮ σ2 = variance is ‘infinite’ (depends on upper cutoff) ◮ Width of distribution is ‘infinite’

slide-26
SLIDE 26

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 31/33

Moments

Standard deviation is a mathematical convenience!:

◮ Variance is nice analytically... ◮ Another measure of distribution width:

Mean average deviation (MAD) = |x − x|

◮ MAD is unpleasant analytically...

slide-27
SLIDE 27

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 32/33

How sample sizes grow...

Given P(x) ∼ cx−γ:

◮ We can show that after n samples, we expect the

largest sample to be x1 n1/(γ−1)

◮ Sampling from a ‘mild’ distribution gives a much

slower growth with n.

◮ e.g., for P(x) = λe−λx, we find

x1 1 λ ln n.

slide-28
SLIDE 28

Power Law Size Distributions Overview

Introduction Examples Zipf’s law Wild vs. Mild CCDFs

References Frame 33/33

References I

  • N. N. Taleb.

The Black Swan. Random House, New York, 2007.

  • G. K. Zipf.

Human Behaviour and the Principle of Least-Effort. Addison-Wesley, Cambridge, MA, 1949.