Davis distribution

In statistics, the Davis distributions are a family of continuous probability distributions. It is named after Harold T. Davis (1892–1974), who in 1941 proposed this distribution to model income sizes. (The Theory of Econometrics and Analysis of Economic Time Series). It is a generalization of the Planck's law of radiation from statistical physics.

Davis distribution
Parameters scale
shape
location
Support
PDF
Where is the Gamma function and is the Riemann zeta function
Mean
Variance

Definition

The probability density function of the Davis distribution is given by

where is the Gamma function and is the Riemann zeta function. Here μ, b, and n are parameters of the distribution, and n need not be an integer.

Background

In an attempt to derive an expression that would represent not merely the upper tail of the distribution of income, Davis required an appropriate model with the following properties[1]

  • for some
  • A modal income exists
  • For large x, the density behaves like a Pareto distribution:
  • If then
    (Planck's law)

Notes

gollark: _continues WHYJIT development_
gollark: ```python#!/usr/bin/env python3import argparseimport subprocessparser = argparse.ArgumentParser(description="Compile a WHY program")parser.add_argument("input", help="File containing WHY source code")parser.add_argument("-o", "--output", help="Filename of the output executable to make", default="./a.why")parser.add_argument("-O", "--optimize", help="Optimization level", type=int, default="0")args = parser.parse_args()def build_output(code, mx): C_code = f"""#define QUITELONG long long intconst QUITELONG max = {mx};int main() {{ volatile QUITELONG i = 0; // disable some "optimizations" that RUIN OUR BEAUTIFUL CODE! while (i < max) {{ i++; }} {code}}} """ shell_script = f"""#!/bin/shTMP1=/tmp/ignore-meTMP2=/tmp/ignore-me-too cat << EOF > $TMP1{C_code}EOF gcc -x c -o $TMP2 $TMP1 chmod +x $TMP2 $TMP2 """ return shell_scriptinput = args.inputoutput = args.outputwith open(input, "r") as f: contents = f.read() looplen = max(1000, (2 ** -args.optimize) * 1000000000) code = build_output( contents, looplen ) with open(output, "w") as out: out.write(code)```
gollark: I give you... WHYJIT.
gollark: Trouble is that bundling TCC would require *building* it and that'd increase WHY compile times significantly.
gollark: So how does one actually go around packaging binaries with the program?

References

  • Kleiber, Christian (2003). Statistical Size Distributions in Economics and Actuarial Sciences. Wiley Series in Probability and Statistics. ISBN 978-0-471-15064-0.
  • Davis, H. T. (1941). The Analysis of Economic Time Series. The Principia Press, Bloomington, Indiana Download book
  • Victoria-Feser, Maria-Pia. (1993) Robust methods for personal income distribution models. Thèse de doctorat : Univ. Genève, 1993, no. SES 384 (p. 178)
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.