On the sum of power laws

Many real-life data sets have power laws or Zipfian distributions. An integer-valued random variable X follows a power law with parameter a if P(X = k) is proportional to ka. Panos asked what the sum of two power laws was. He cites Wilke at al. who claim that the sum of two power laws X and Y with parameters a and b is a power law with parameter min(ab).

I relate this problem to the sum of exponentials. Any engineer knows that if a>b, then eat + ebt will be approximately eat for t sufficiently large. Hence, the sum of power law distributions X and Y is a power law distribution with parameter min(ab) if you are only interested in large values of k in P(X + Y = k).

However, the sum of two power laws is not a power law. Egghe showed in The distribution of N-grams that even if the words follow a power law, the n-grams won’t!

Published by

Daniel Lemire

A computer science professor at the University of Quebec (TELUQ).

3 thoughts on “On the sum of power laws”

  1. It seems that power law is really the same as Pareto distribution . This paper gives some closed formulas for distributions of sums of Pareto, which are themselves not Pareto

    If two power laws have different parameters, as you go to infinity, odds of encountering the one with higher a becomes vs. one with lower a goes to 0, so I also expect that for large values, heavier tail distribution will dominate

    BTW, I also wondered about distribution of bigrams when unigrams are power-law distributed, David Cantrell in sci.math gave an approximate formula for the cdf involving Lambert’s W function
    http://groups.google.com/group/sci.math/browse_thread/thread/8de7cee65f65ff70/810470b85f36523b?lnk=st&q=group%3Asci.math#810470b85f36523b

Leave a Reply

Your email address will not be published. The comment form expects plain text. If you need to format your text, you can use HTML elements such strong, blockquote, cite, code and em. For formatting code as HTML automatically, I recommend tohtml.com.

You may subscribe to this blog by email.