3 <meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
4 <title>Pareto Distribution</title>
5 <link rel="stylesheet" href="../../../math.css" type="text/css">
6 <meta name="generator" content="DocBook XSL Stylesheets V1.79.1">
7 <link rel="home" href="../../../index.html" title="Math Toolkit 2.11.0">
8 <link rel="up" href="../dists.html" title="Distributions">
9 <link rel="prev" href="normal_dist.html" title="Normal (Gaussian) Distribution">
10 <link rel="next" href="poisson_dist.html" title="Poisson Distribution">
12 <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
13 <table cellpadding="2" width="100%"><tr>
14 <td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../../boost.png"></td>
15 <td align="center"><a href="../../../../../../../index.html">Home</a></td>
16 <td align="center"><a href="../../../../../../../libs/libraries.htm">Libraries</a></td>
17 <td align="center"><a href="http://www.boost.org/users/people.html">People</a></td>
18 <td align="center"><a href="http://www.boost.org/users/faq.html">FAQ</a></td>
19 <td align="center"><a href="../../../../../../../more/index.htm">More</a></td>
22 <div class="spirit-nav">
23 <a accesskey="p" href="normal_dist.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../dists.html"><img src="../../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../../index.html"><img src="../../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="poisson_dist.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
26 <div class="titlepage"><div><div><h4 class="title">
27 <a name="math_toolkit.dist_ref.dists.pareto"></a><a class="link" href="pareto.html" title="Pareto Distribution">Pareto Distribution</a>
28 </h4></div></div></div>
29 <pre class="programlisting"><span class="preprocessor">#include</span> <span class="special"><</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">math</span><span class="special">/</span><span class="identifier">distributions</span><span class="special">/</span><span class="identifier">pareto</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">></span></pre>
30 <pre class="programlisting"><span class="keyword">namespace</span> <span class="identifier">boost</span><span class="special">{</span> <span class="keyword">namespace</span> <span class="identifier">math</span><span class="special">{</span>
32 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">RealType</span> <span class="special">=</span> <span class="keyword">double</span><span class="special">,</span>
33 <span class="keyword">class</span> <a class="link" href="../../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a> <span class="special">=</span> <a class="link" href="../../pol_ref/pol_ref_ref.html" title="Policy Class Reference">policies::policy<></a> <span class="special">></span>
34 <span class="keyword">class</span> <span class="identifier">pareto_distribution</span><span class="special">;</span>
36 <span class="keyword">typedef</span> <span class="identifier">pareto_distribution</span><span class="special"><></span> <span class="identifier">pareto</span><span class="special">;</span>
38 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">RealType</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
39 <span class="keyword">class</span> <span class="identifier">pareto_distribution</span>
40 <span class="special">{</span>
41 <span class="keyword">public</span><span class="special">:</span>
42 <span class="keyword">typedef</span> <span class="identifier">RealType</span> <span class="identifier">value_type</span><span class="special">;</span>
43 <span class="comment">// Constructor:</span>
44 <span class="identifier">pareto_distribution</span><span class="special">(</span><span class="identifier">RealType</span> <span class="identifier">scale</span> <span class="special">=</span> <span class="number">1</span><span class="special">,</span> <span class="identifier">RealType</span> <span class="identifier">shape</span> <span class="special">=</span> <span class="number">1</span><span class="special">)</span>
45 <span class="comment">// Accessors:</span>
46 <span class="identifier">RealType</span> <span class="identifier">scale</span><span class="special">()</span><span class="keyword">const</span><span class="special">;</span>
47 <span class="identifier">RealType</span> <span class="identifier">shape</span><span class="special">()</span><span class="keyword">const</span><span class="special">;</span>
48 <span class="special">};</span>
50 <span class="special">}}</span> <span class="comment">// namespaces</span>
53 The <a href="http://en.wikipedia.org/wiki/pareto_distribution" target="_top">Pareto
54 distribution</a> is a continuous distribution with the <a href="http://en.wikipedia.org/wiki/Probability_density_function" target="_top">probability
55 density function (pdf)</a>:
57 <div class="blockquote"><blockquote class="blockquote"><p>
58 <span class="serif_italic">f(x; α, β) = αβ<sup>α</sup> / x<sup>α+ 1</sup></span>
59 </p></blockquote></div>
61 For shape parameter α > 0, and scale parameter β > 0. If x < β, the
65 The <a href="http://mathworld.wolfram.com/ParetoDistribution.html" target="_top">Pareto
66 distribution</a> often describes the larger compared to the smaller.
67 A classic example is that 80% of the wealth is owned by 20% of the population.
70 The following graph illustrates how the PDF varies with the scale parameter
73 <div class="blockquote"><blockquote class="blockquote"><p>
74 <span class="inlinemediaobject"><img src="../../../../graphs/pareto_pdf1.svg" align="middle"></span>
76 </p></blockquote></div>
78 And this graph illustrates how the PDF varies with the shape parameter
81 <div class="blockquote"><blockquote class="blockquote"><p>
82 <span class="inlinemediaobject"><img src="../../../../graphs/pareto_pdf2.svg" align="middle"></span>
84 </p></blockquote></div>
86 <a name="math_toolkit.dist_ref.dists.pareto.h0"></a>
87 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.related_distributions"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.related_distributions">Related
91 <a name="math_toolkit.dist_ref.dists.pareto.h1"></a>
92 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.member_functions"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.member_functions">Member Functions</a>
94 <pre class="programlisting"><span class="identifier">pareto_distribution</span><span class="special">(</span><span class="identifier">RealType</span> <span class="identifier">scale</span> <span class="special">=</span> <span class="number">1</span><span class="special">,</span> <span class="identifier">RealType</span> <span class="identifier">shape</span> <span class="special">=</span> <span class="number">1</span><span class="special">);</span>
97 Constructs a <a href="http://en.wikipedia.org/wiki/pareto_distribution" target="_top">pareto
98 distribution</a> with shape <span class="emphasis"><em>shape</em></span> and scale <span class="emphasis"><em>scale</em></span>.
101 Requires that the <span class="emphasis"><em>shape</em></span> and <span class="emphasis"><em>scale</em></span>
102 parameters are both greater than zero, otherwise calls <a class="link" href="../../error_handling.html#math_toolkit.error_handling.domain_error">domain_error</a>.
104 <pre class="programlisting"><span class="identifier">RealType</span> <span class="identifier">scale</span><span class="special">()</span><span class="keyword">const</span><span class="special">;</span>
107 Returns the <span class="emphasis"><em>scale</em></span> parameter of this distribution.
109 <pre class="programlisting"><span class="identifier">RealType</span> <span class="identifier">shape</span><span class="special">()</span><span class="keyword">const</span><span class="special">;</span>
112 Returns the <span class="emphasis"><em>shape</em></span> parameter of this distribution.
115 <a name="math_toolkit.dist_ref.dists.pareto.h2"></a>
116 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.non_member_accessors"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.non_member_accessors">Non-member
120 All the <a class="link" href="../nmp.html" title="Non-Member Properties">usual non-member accessor
121 functions</a> that are generic to all distributions are supported:
122 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.cdf">Cumulative Distribution Function</a>,
123 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.pdf">Probability Density Function</a>,
124 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.quantile">Quantile</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.hazard">Hazard Function</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.chf">Cumulative Hazard Function</a>,
125 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.mean">mean</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.median">median</a>,
126 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.mode">mode</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.variance">variance</a>,
127 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.sd">standard deviation</a>,
128 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.skewness">skewness</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.kurtosis">kurtosis</a>, <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.kurtosis_excess">kurtosis_excess</a>,
129 <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.range">range</a> and <a class="link" href="../nmp.html#math_toolkit.dist_ref.nmp.support">support</a>.
132 The supported domain of the random variable is [scale, ∞].
135 <a name="math_toolkit.dist_ref.dists.pareto.h3"></a>
136 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.accuracy"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.accuracy">Accuracy</a>
139 The Pareto distribution is implemented in terms of the standard library
140 <code class="computeroutput"><span class="identifier">exp</span></code> functions plus <a class="link" href="../../powers/expm1.html" title="expm1">expm1</a> and so should have very
141 small errors, usually only a few epsilon.
144 If probability is near to unity (or the complement of a probability near
145 zero) see also <a class="link" href="../../stat_tut/overview/complements.html#why_complements">why complements?</a>.
148 <a name="math_toolkit.dist_ref.dists.pareto.h4"></a>
149 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.implementation"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.implementation">Implementation</a>
152 In the following table α is the shape parameter of the distribution, and
153 β is its scale parameter, <span class="emphasis"><em>x</em></span> is the random variate,
154 <span class="emphasis"><em>p</em></span> is the probability and its complement <span class="emphasis"><em>q
157 <div class="informaltable"><table class="table">
183 Using the relation: pdf p = αβ<sup>α</sup>/x<sup>α +1</sup>
195 Using the relation: cdf p = 1 - (β / x)<sup>α</sup>
207 Using the relation: q = 1 - p = -(β / x)<sup>α</sup>
219 Using the relation: x = β / (1 - p)<sup>1/α</sup>
226 quantile from the complement
231 Using the relation: x = β / (q)<sup>1/α</sup>
243 αβ / (β - 1)
255 βα<sup>2</sup> / (β - 1)<sup>2</sup> (β - 2)
279 Refer to <a href="http://mathworld.wolfram.com/ParetoDistribution.html" target="_top">Weisstein,
280 Eric W. "Pareto Distribution." From MathWorld--A Wolfram
293 Refer to <a href="http://mathworld.wolfram.com/ParetoDistribution.html" target="_top">Weisstein,
294 Eric W. "Pareto Distribution." From MathWorld--A Wolfram
307 Refer to <a href="http://mathworld.wolfram.com/ParetoDistribution.html" target="_top">Weisstein,
308 Eric W. "pareto Distribution." From MathWorld--A Wolfram
316 <a name="math_toolkit.dist_ref.dists.pareto.h5"></a>
317 <span class="phrase"><a name="math_toolkit.dist_ref.dists.pareto.references"></a></span><a class="link" href="pareto.html#math_toolkit.dist_ref.dists.pareto.references">References</a>
319 <div class="itemizedlist"><ul class="itemizedlist" style="list-style-type: disc; ">
320 <li class="listitem">
321 <a href="http://en.wikipedia.org/wiki/pareto_distribution" target="_top">Pareto
324 <li class="listitem">
325 <a href="http://mathworld.wolfram.com/paretoDistribution.html" target="_top">Weisstein,
326 Eric W. "Pareto Distribution." From MathWorld--A Wolfram
329 <li class="listitem">
330 Handbook of Statistical Distributions with Applications, K Krishnamoorthy,
331 ISBN 1-58488-635-8, Chapter 23, pp 257 - 267. (Note the meaning of
332 a and b is reversed in Wolfram and Krishnamoorthy).
336 <table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
337 <td align="left"></td>
338 <td align="right"><div class="copyright-footer">Copyright © 2006-2019 Nikhar
339 Agrawal, Anton Bikineev, Paul A. Bristow, Marco Guazzone, Christopher Kormanyos,
340 Hubert Holin, Bruno Lalande, John Maddock, Jeremy Murphy, Matthew Pulver, Johan
341 Råde, Gautam Sewani, Benjamin Sobotta, Nicholas Thompson, Thijs van den Berg,
342 Daryle Walker and Xiaogang Zhang<p>
343 Distributed under the Boost Software License, Version 1.0. (See accompanying
344 file LICENSE_1_0.txt or copy at <a href="http://www.boost.org/LICENSE_1_0.txt" target="_top">http://www.boost.org/LICENSE_1_0.txt</a>)
349 <div class="spirit-nav">
350 <a accesskey="p" href="normal_dist.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../dists.html"><img src="../../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../../index.html"><img src="../../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="poisson_dist.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>