Boost C++ Libraries: Ticket #2651: Documentation errors, potential application problems https://svn.boost.org/trac10/ticket/2651 <p> Some formulae in the documentation of pareto_distribution are wrong, e.g. the true quantile is beta/(1-p)<sup>(1/alpha). </sup></p> <p> Further I have a remark on the implementation. According to the documentation page, computation errors are low except p is near to 1. Unfortunately, this is just the region of interest in extreme value applications, and the Pareto distribution is widely used there. </p> <p> I have not tested the implementation yet, but there may be some problems with extreme value applications. </p> en-us Boost C++ Libraries /htdocs/site/boost.png https://svn.boost.org/trac10/ticket/2651 Trac 1.4.3 Steven Watanabe Wed, 14 Jan 2009 03:06:16 GMT component changed; owner set https://svn.boost.org/trac10/ticket/2651#comment:1 https://svn.boost.org/trac10/ticket/2651#comment:1 <ul> <li><strong>owner</strong> set to <span class="trac-author">John Maddock</span> </li> <li><strong>component</strong> <span class="trac-field-old">None</span> → <span class="trac-field-new">math</span> </li> </ul> Ticket John Maddock Wed, 14 Jan 2009 10:29:16 GMT status changed; cc set https://svn.boost.org/trac10/ticket/2651#comment:2 https://svn.boost.org/trac10/ticket/2651#comment:2 <ul> <li><strong>cc</strong> <span class="trac-author">pbristow@…</span> added </li> <li><strong>status</strong> <span class="trac-field-old">new</span> → <span class="trac-field-new">assigned</span> </li> </ul> <p> I'm CC'ing Paul Bristow on this one as I believe he wrote this one. </p> <p> The quantile in the docs looks like a simple type (alpha and beta wrong way around). </p> <p> Re: errors for p near 1, I'm not sure I understand this remark - Paul can you clarify what the issue is if any? </p> <p> John. </p> Ticket John Maddock Fri, 23 Jan 2009 11:46:50 GMT <link>https://svn.boost.org/trac10/ticket/2651#comment:3 </link> <guid isPermaLink="false">https://svn.boost.org/trac10/ticket/2651#comment:3</guid> <description> <p> Just to clarify this is an issue with getting the docs wording correct, not with the implementation. </p> <ul><li>Asuming the p values passed to the quantile function are exact binary values then there is no issue with accuracy, but... </li><li>As with all distributions if you pass a p value close to 1, and that value is inexact, then what we're really need in order to calculate the quantile is 1-p which will be mostly garbage digits for p ~ 1. </li></ul><p> To put this another way: </p> <p> quantile(my_pareto, 0.99) </p> <p> is potentially less accurate than </p> <p> quantile(complement(my_pareto, 0.01)) </p> <p> since 1 - 0.99 != 0.01 when using fp-arithmetic. </p> <p> This assumes that the 0.01 value (for example) is either a constant, or computed by some manner other than subtracting 0.99 from 1. </p> <p> Not sure at present, whether we should reword this, or just remove the accuracy caveat from the docs, given that this is an issue with all distributions, and all implementations thereof not just ours. </p> <p> HTH, John. </p> </description> <category>Ticket</category> </item> <item> <dc:creator>Paul A. Bristow</dc:creator> <pubDate>Wed, 28 Jan 2009 16:57:14 GMT</pubDate> <title>status changed; resolution set https://svn.boost.org/trac10/ticket/2651#comment:4 https://svn.boost.org/trac10/ticket/2651#comment:4 <ul> <li><strong>status</strong> <span class="trac-field-old">assigned</span> → <span class="trac-field-new">closed</span> </li> <li><strong>resolution</strong> → <span class="trac-field-new">fixed</span> </li> </ul> <p> Pareto errors corrected and major revisions to make need for complements more obvious, in TOC, and indexed. </p> Ticket