Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/GumbelDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/GumbelDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/GumbelDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,61 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>GumbelDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/GumbelDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <em class="jxr_javadoccomment"> * Test cases for {@link GumbelDistribution}.</em> +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">class</strong> <a name="GumbelDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/GumbelDistributionTest.html#GumbelDistributionTest">GumbelDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L24" href="#L24">24</a> @Override +<a class="jxr_linenumber" name="L25" href="#L25">25</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> location = (Double) parameters[0]; +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = (Double) parameters[1]; +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <strong class="jxr_keyword">return</strong> GumbelDistribution.of(location, scale); +<a class="jxr_linenumber" name="L29" href="#L29">29</a> } +<a class="jxr_linenumber" name="L30" href="#L30">30</a> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> @Override +<a class="jxr_linenumber" name="L32" href="#L32">32</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L34" href="#L34">34</a> {10.0, 0.0}, +<a class="jxr_linenumber" name="L35" href="#L35">35</a> {10.0, -0.1} +<a class="jxr_linenumber" name="L36" href="#L36">36</a> }; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> } +<a class="jxr_linenumber" name="L38" href="#L38">38</a> +<a class="jxr_linenumber" name="L39" href="#L39">39</a> @Override +<a class="jxr_linenumber" name="L40" href="#L40">40</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L41" href="#L41">41</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Location"</span>, <span class="jxr_string">"Scale"</span>}; +<a class="jxr_linenumber" name="L42" href="#L42">42</a> } +<a class="jxr_linenumber" name="L43" href="#L43">43</a> +<a class="jxr_linenumber" name="L44" href="#L44">44</a> @Override +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L47" href="#L47">47</a> } +<a class="jxr_linenumber" name="L48" href="#L48">48</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html>
Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/HypergeometricDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/HypergeometricDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/HypergeometricDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,260 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>HypergeometricDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/HypergeometricDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> +<a class="jxr_linenumber" name="L18" href="#L18">18</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L19" href="#L19">19</a> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> java.util.stream.Stream; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.apache.commons.numbers.core.Precision; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.apache.commons.rng.simple.RandomSource; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments; +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.MethodSource; +<a class="jxr_linenumber" name="L28" href="#L28">28</a> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <em class="jxr_javadoccomment"> * Test cases for {@link HypergeometricDistribution}.</em> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment"> * Extends {@link BaseDiscreteDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <strong class="jxr_keyword">class</strong> <a name="HypergeometricDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/HypergeometricDistributionTest.html#HypergeometricDistributionTest">HypergeometricDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseDiscreteDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseDiscreteDistributionTest.html#BaseDiscreteDistributionTest">BaseDiscreteDistributionTest</a> { +<a class="jxr_linenumber" name="L34" href="#L34">34</a> @Override +<a class="jxr_linenumber" name="L35" href="#L35">35</a> DiscreteDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> populationSize = (Integer) parameters[0]; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> numberOfSuccesses = (Integer) parameters[1]; +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> sampleSize = (Integer) parameters[2]; +<a class="jxr_linenumber" name="L39" href="#L39">39</a> <strong class="jxr_keyword">return</strong> HypergeometricDistribution.of(populationSize, numberOfSuccesses, sampleSize); +<a class="jxr_linenumber" name="L40" href="#L40">40</a> } +<a class="jxr_linenumber" name="L41" href="#L41">41</a> +<a class="jxr_linenumber" name="L42" href="#L42">42</a> @Override +<a class="jxr_linenumber" name="L43" href="#L43">43</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L44" href="#L44">44</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L45" href="#L45">45</a> {0, 3, 5}, +<a class="jxr_linenumber" name="L46" href="#L46">46</a> {-1, 3, 5}, +<a class="jxr_linenumber" name="L47" href="#L47">47</a> {5, -1, 5}, +<a class="jxr_linenumber" name="L48" href="#L48">48</a> {5, 3, -1}, +<a class="jxr_linenumber" name="L49" href="#L49">49</a> {5, 6, 5}, +<a class="jxr_linenumber" name="L50" href="#L50">50</a> {5, 3, 6}, +<a class="jxr_linenumber" name="L51" href="#L51">51</a> }; +<a class="jxr_linenumber" name="L52" href="#L52">52</a> } +<a class="jxr_linenumber" name="L53" href="#L53">53</a> +<a class="jxr_linenumber" name="L54" href="#L54">54</a> @Override +<a class="jxr_linenumber" name="L55" href="#L55">55</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L56" href="#L56">56</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"PopulationSize"</span>, <span class="jxr_string">"NumberOfSuccesses"</span>, <span class="jxr_string">"SampleSize"</span>}; +<a class="jxr_linenumber" name="L57" href="#L57">57</a> } +<a class="jxr_linenumber" name="L58" href="#L58">58</a> +<a class="jxr_linenumber" name="L59" href="#L59">59</a> @Override +<a class="jxr_linenumber" name="L60" href="#L60">60</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L61" href="#L61">61</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L62" href="#L62">62</a> } +<a class="jxr_linenumber" name="L63" href="#L63">63</a> +<a class="jxr_linenumber" name="L64" href="#L64">64</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L65" href="#L65">65</a> +<a class="jxr_linenumber" name="L66" href="#L66">66</a> @ParameterizedTest +<a class="jxr_linenumber" name="L67" href="#L67">67</a> @MethodSource +<a class="jxr_linenumber" name="L68" href="#L68">68</a> <strong class="jxr_keyword">void</strong> testAdditionalMoments(<strong class="jxr_keyword">int</strong> populationSize, +<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">int</strong> numberOfSuccesses, +<a class="jxr_linenumber" name="L70" href="#L70">70</a> <strong class="jxr_keyword">int</strong> sampleSize, +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">double</strong> mean, <strong class="jxr_keyword">double</strong> variance) { +<a class="jxr_linenumber" name="L72" href="#L72">72</a> <strong class="jxr_keyword">final</strong> HypergeometricDistribution dist = HypergeometricDistribution.of(populationSize, numberOfSuccesses, sampleSize); +<a class="jxr_linenumber" name="L73" href="#L73">73</a> testMoments(dist, mean, variance, DoubleTolerances.ulps(1)); +<a class="jxr_linenumber" name="L74" href="#L74">74</a> } +<a class="jxr_linenumber" name="L75" href="#L75">75</a> +<a class="jxr_linenumber" name="L76" href="#L76">76</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionalMoments() { +<a class="jxr_linenumber" name="L77" href="#L77">77</a> <strong class="jxr_keyword">return</strong> Stream.of( +<a class="jxr_linenumber" name="L78" href="#L78">78</a> Arguments.of(1500, 40, 100, 40d * 100d / 1500d, (100d * 40d * (1500d - 100d) * (1500d - 40d)) / ((1500d * 1500d * 1499d))), +<a class="jxr_linenumber" name="L79" href="#L79">79</a> Arguments.of(3000, 55, 200, 55d * 200d / 3000d, (200d * 55d * (3000d - 200d) * (3000d - 55d)) / ((3000d * 3000d * 2999d))) +<a class="jxr_linenumber" name="L80" href="#L80">80</a> ); +<a class="jxr_linenumber" name="L81" href="#L81">81</a> } +<a class="jxr_linenumber" name="L82" href="#L82">82</a> +<a class="jxr_linenumber" name="L83" href="#L83">83</a> @Test +<a class="jxr_linenumber" name="L84" href="#L84">84</a> <strong class="jxr_keyword">void</strong> testLargeValues() { +<a class="jxr_linenumber" name="L85" href="#L85">85</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> populationSize = 3456; +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> sampleSize = 789; +<a class="jxr_linenumber" name="L87" href="#L87">87</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> numberOfSucceses = 101; +<a class="jxr_linenumber" name="L88" href="#L88">88</a> <em class="jxr_comment">// data[i][3] contains P(x >= x).</em> +<a class="jxr_linenumber" name="L89" href="#L89">89</a> <em class="jxr_comment">// It is tested using survivalProbability(x - 1)</em> +<a class="jxr_linenumber" name="L90" href="#L90">90</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[][] data = { +<a class="jxr_linenumber" name="L91" href="#L91">91</a> {0.0, 2.75646034603961e-12, 2.75646034603961e-12, 1.0}, +<a class="jxr_linenumber" name="L92" href="#L92">92</a> {1.0, 8.55705370142386e-11, 8.83269973602783e-11, 0.999999999997244}, +<a class="jxr_linenumber" name="L93" href="#L93">93</a> {2.0, 1.31288129219665e-9, 1.40120828955693e-9, 0.999999999911673}, +<a class="jxr_linenumber" name="L94" href="#L94">94</a> {3.0, 1.32724172984193e-8, 1.46736255879763e-8, 0.999999998598792}, +<a class="jxr_linenumber" name="L95" href="#L95">95</a> {4.0, 9.94501711734089e-8, 1.14123796761385e-7, 0.999999985326375}, +<a class="jxr_linenumber" name="L96" href="#L96">96</a> {5.0, 5.89080768883643e-7, 7.03204565645028e-7, 0.999999885876203}, +<a class="jxr_linenumber" name="L97" href="#L97">97</a> {20.0, 0.0760051397707708, 0.27349758476299, 0.802507555007781}, +<a class="jxr_linenumber" name="L98" href="#L98">98</a> {21.0, 0.087144222047629, 0.360641806810619, 0.72650241523701}, +<a class="jxr_linenumber" name="L99" href="#L99">99</a> {22.0, 0.0940378846881819, 0.454679691498801, 0.639358193189381}, +<a class="jxr_linenumber" name="L100" href="#L100">100</a> {23.0, 0.0956897500614809, 0.550369441560282, 0.545320308501199}, +<a class="jxr_linenumber" name="L101" href="#L101">101</a> {24.0, 0.0919766921922999, 0.642346133752582, 0.449630558439718}, +<a class="jxr_linenumber" name="L102" href="#L102">102</a> {25.0, 0.083641637261095, 0.725987771013677, 0.357653866247418}, +<a class="jxr_linenumber" name="L103" href="#L103">103</a> {96.0, 5.93849188852098e-57, 1.0, 6.01900244560712e-57}, +<a class="jxr_linenumber" name="L104" href="#L104">104</a> {97.0, 7.96593036832547e-59, 1.0, 8.05105570861321e-59}, +<a class="jxr_linenumber" name="L105" href="#L105">105</a> {98.0, 8.44582921934367e-61, 1.0, 8.5125340287733e-61}, +<a class="jxr_linenumber" name="L106" href="#L106">106</a> {99.0, 6.63604297068222e-63, 1.0, 6.670480942963e-63}, +<a class="jxr_linenumber" name="L107" href="#L107">107</a> {100.0, 3.43501099007557e-65, 1.0, 3.4437972280786e-65}, +<a class="jxr_linenumber" name="L108" href="#L108">108</a> {101.0, 8.78623800302957e-68, 1.0, 8.78623800302957e-68}, +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <em class="jxr_comment">// Out of domain</em> +<a class="jxr_linenumber" name="L110" href="#L110">110</a> {sampleSize + 1, 0, 1.0, 0}, +<a class="jxr_linenumber" name="L111" href="#L111">111</a> }; +<a class="jxr_linenumber" name="L112" href="#L112">112</a> +<a class="jxr_linenumber" name="L113" href="#L113">113</a> testHypergeometricDistributionProbabilities(populationSize, sampleSize, numberOfSucceses, data); +<a class="jxr_linenumber" name="L114" href="#L114">114</a> } +<a class="jxr_linenumber" name="L115" href="#L115">115</a> +<a class="jxr_linenumber" name="L116" href="#L116">116</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> testHypergeometricDistributionProbabilities(<strong class="jxr_keyword">int</strong> populationSize, <strong class="jxr_keyword">int</strong> sampleSize, +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <strong class="jxr_keyword">int</strong> numberOfSuccesses, <strong class="jxr_keyword">double</strong>[][] data) { +<a class="jxr_linenumber" name="L118" href="#L118">118</a> <strong class="jxr_keyword">final</strong> HypergeometricDistribution dist = HypergeometricDistribution.of(populationSize, numberOfSuccesses, sampleSize); +<a class="jxr_linenumber" name="L119" href="#L119">119</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < data.length; ++i) { +<a class="jxr_linenumber" name="L120" href="#L120">120</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> x = (<strong class="jxr_keyword">int</strong>)data[i][0]; +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> pmf = data[i][1]; +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actualPmf = dist.probability(x); +<a class="jxr_linenumber" name="L123" href="#L123">123</a> TestUtils.assertRelativelyEquals(() -> <span class="jxr_string">"Expected equals for <"</span> + x + <span class="jxr_string">"> pmf"</span>, pmf, actualPmf, 1.0e-9); +<a class="jxr_linenumber" name="L124" href="#L124">124</a> +<a class="jxr_linenumber" name="L125" href="#L125">125</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> cdf = data[i][2]; +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actualCdf = dist.cumulativeProbability(x); +<a class="jxr_linenumber" name="L127" href="#L127">127</a> TestUtils.assertRelativelyEquals(() -> <span class="jxr_string">"Expected equals for <"</span> + x + <span class="jxr_string">"> cdf"</span>, cdf, actualCdf, 1.0e-9); +<a class="jxr_linenumber" name="L128" href="#L128">128</a> +<a class="jxr_linenumber" name="L129" href="#L129">129</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> cdf1 = data[i][3]; +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actualCdf1 = dist.survivalProbability(x - 1); +<a class="jxr_linenumber" name="L131" href="#L131">131</a> TestUtils.assertRelativelyEquals(() -> <span class="jxr_string">"Expected equals for <"</span> + x + <span class="jxr_string">"> cdf1"</span>, cdf1, actualCdf1, 1.0e-9); +<a class="jxr_linenumber" name="L132" href="#L132">132</a> } +<a class="jxr_linenumber" name="L133" href="#L133">133</a> } +<a class="jxr_linenumber" name="L134" href="#L134">134</a> +<a class="jxr_linenumber" name="L135" href="#L135">135</a> @Test +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">void</strong> testMoreLargeValues() { +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> populationSize = 26896; +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> sampleSize = 895; +<a class="jxr_linenumber" name="L139" href="#L139">139</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> numberOfSucceses = 55; +<a class="jxr_linenumber" name="L140" href="#L140">140</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[][] data = { +<a class="jxr_linenumber" name="L141" href="#L141">141</a> {0.0, 0.155168304750504, 0.155168304750504, 1.0}, +<a class="jxr_linenumber" name="L142" href="#L142">142</a> {1.0, 0.29437545000746, 0.449543754757964, 0.844831695249496}, +<a class="jxr_linenumber" name="L143" href="#L143">143</a> {2.0, 0.273841321577003, 0.723385076334967, 0.550456245242036}, +<a class="jxr_linenumber" name="L144" href="#L144">144</a> {3.0, 0.166488572570786, 0.889873648905753, 0.276614923665033}, +<a class="jxr_linenumber" name="L145" href="#L145">145</a> {4.0, 0.0743969744713231, 0.964270623377076, 0.110126351094247}, +<a class="jxr_linenumber" name="L146" href="#L146">146</a> {5.0, 0.0260542785784855, 0.990324901955562, 0.0357293766229237}, +<a class="jxr_linenumber" name="L147" href="#L147">147</a> {20.0, 3.57101101678792e-16, 1.0, 3.78252101622096e-16}, +<a class="jxr_linenumber" name="L148" href="#L148">148</a> {21.0, 2.00551638598312e-17, 1.0, 2.11509999433041e-17}, +<a class="jxr_linenumber" name="L149" href="#L149">149</a> {22.0, 1.04317070180562e-18, 1.0, 1.09583608347287e-18}, +<a class="jxr_linenumber" name="L150" href="#L150">150</a> {23.0, 5.03153504903308e-20, 1.0, 5.266538166725e-20}, +<a class="jxr_linenumber" name="L151" href="#L151">151</a> {24.0, 2.2525984149695e-21, 1.0, 2.35003117691919e-21}, +<a class="jxr_linenumber" name="L152" href="#L152">152</a> {25.0, 9.3677424515947e-23, 1.0, 9.74327619496943e-23}, +<a class="jxr_linenumber" name="L153" href="#L153">153</a> {50.0, 9.83633962945521e-69, 1.0, 9.8677629437617e-69}, +<a class="jxr_linenumber" name="L154" href="#L154">154</a> {51.0, 3.13448949497553e-71, 1.0, 3.14233143064882e-71}, +<a class="jxr_linenumber" name="L155" href="#L155">155</a> {52.0, 7.82755221928122e-74, 1.0, 7.84193567329055e-74}, +<a class="jxr_linenumber" name="L156" href="#L156">156</a> {53.0, 1.43662126065532e-76, 1.0, 1.43834540093295e-76}, +<a class="jxr_linenumber" name="L157" href="#L157">157</a> {54.0, 1.72312692517348e-79, 1.0, 1.7241402776278e-79}, +<a class="jxr_linenumber" name="L158" href="#L158">158</a> {55.0, 1.01335245432581e-82, 1.0, 1.01335245432581e-82}, +<a class="jxr_linenumber" name="L159" href="#L159">159</a> }; +<a class="jxr_linenumber" name="L160" href="#L160">160</a> testHypergeometricDistributionProbabilities(populationSize, sampleSize, numberOfSucceses, data); +<a class="jxr_linenumber" name="L161" href="#L161">161</a> } +<a class="jxr_linenumber" name="L162" href="#L162">162</a> +<a class="jxr_linenumber" name="L163" href="#L163">163</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L164" href="#L164">164</a> <em class="jxr_javadoccomment"> * Test Math-644 is ported from Commons Math 3 where the distribution had the function</em> +<a class="jxr_linenumber" name="L165" href="#L165">165</a> <em class="jxr_javadoccomment"> * upperCumulativeProbability(x) to compute P(X >= x). This has been replaced</em> +<a class="jxr_linenumber" name="L166" href="#L166">166</a> <em class="jxr_javadoccomment"> * in Commons Statistics with survivalProbability(x) which computes P(X > x). To</em> +<a class="jxr_linenumber" name="L167" href="#L167">167</a> <em class="jxr_javadoccomment"> * create the equivalent use survivalProbability(x - 1).</em> +<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L169" href="#L169">169</a> @Test +<a class="jxr_linenumber" name="L170" href="#L170">170</a> <strong class="jxr_keyword">void</strong> testMath644() { +<a class="jxr_linenumber" name="L171" href="#L171">171</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> N = 14761461; <em class="jxr_comment">// population</em> +<a class="jxr_linenumber" name="L172" href="#L172">172</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m = 1035; <em class="jxr_comment">// successes in population</em> +<a class="jxr_linenumber" name="L173" href="#L173">173</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = 1841; <em class="jxr_comment">// number of trials</em> +<a class="jxr_linenumber" name="L174" href="#L174">174</a> +<a class="jxr_linenumber" name="L175" href="#L175">175</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> k = 0; +<a class="jxr_linenumber" name="L176" href="#L176">176</a> <strong class="jxr_keyword">final</strong> HypergeometricDistribution dist = HypergeometricDistribution.of(N, m, n); +<a class="jxr_linenumber" name="L177" href="#L177">177</a> +<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_comment">// Compute upper cumulative probability using the survival probability</em> +<a class="jxr_linenumber" name="L179" href="#L179">179</a> Assertions.assertEquals(0, Precision.compareTo(1.0, dist.survivalProbability(k - 1), 1)); +<a class="jxr_linenumber" name="L180" href="#L180">180</a> Assertions.assertTrue(Precision.compareTo(dist.cumulativeProbability(k), 0.0, 1) > 0); +<a class="jxr_linenumber" name="L181" href="#L181">181</a> +<a class="jxr_linenumber" name="L182" href="#L182">182</a> <em class="jxr_comment">// another way to calculate the upper cumulative probability</em> +<a class="jxr_linenumber" name="L183" href="#L183">183</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> upper = 1.0 - dist.cumulativeProbability(k) + dist.probability(k); +<a class="jxr_linenumber" name="L184" href="#L184">184</a> Assertions.assertEquals(0, Precision.compareTo(1.0, upper, 1)); +<a class="jxr_linenumber" name="L185" href="#L185">185</a> } +<a class="jxr_linenumber" name="L186" href="#L186">186</a> +<a class="jxr_linenumber" name="L187" href="#L187">187</a> @Test +<a class="jxr_linenumber" name="L188" href="#L188">188</a> <strong class="jxr_keyword">void</strong> testZeroTrials() { +<a class="jxr_linenumber" name="L189" href="#L189">189</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = 11; <em class="jxr_comment">// population</em> +<a class="jxr_linenumber" name="L190" href="#L190">190</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m = 4; <em class="jxr_comment">// successes in population</em> +<a class="jxr_linenumber" name="L191" href="#L191">191</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> s = 0; <em class="jxr_comment">// number of trials</em> +<a class="jxr_linenumber" name="L192" href="#L192">192</a> +<a class="jxr_linenumber" name="L193" href="#L193">193</a> <strong class="jxr_keyword">final</strong> HypergeometricDistribution dist = HypergeometricDistribution.of(n, m, s); +<a class="jxr_linenumber" name="L194" href="#L194">194</a> +<a class="jxr_linenumber" name="L195" href="#L195">195</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 1; i <= n; i++) { +<a class="jxr_linenumber" name="L196" href="#L196">196</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p = dist.probability(i); +<a class="jxr_linenumber" name="L197" href="#L197">197</a> Assertions.assertEquals(0, p, () -> <span class="jxr_string">"p="</span> + p); +<a class="jxr_linenumber" name="L198" href="#L198">198</a> } +<a class="jxr_linenumber" name="L199" href="#L199">199</a> } +<a class="jxr_linenumber" name="L200" href="#L200">200</a> +<a class="jxr_linenumber" name="L201" href="#L201">201</a> @Test +<a class="jxr_linenumber" name="L202" href="#L202">202</a> <strong class="jxr_keyword">void</strong> testMath1356() { +<a class="jxr_linenumber" name="L203" href="#L203">203</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = 11; <em class="jxr_comment">// population</em> +<a class="jxr_linenumber" name="L204" href="#L204">204</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m = 11; <em class="jxr_comment">// successes in population</em> +<a class="jxr_linenumber" name="L205" href="#L205">205</a> +<a class="jxr_linenumber" name="L206" href="#L206">206</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> s = 0; s <= n; s++) { +<a class="jxr_linenumber" name="L207" href="#L207">207</a> <strong class="jxr_keyword">final</strong> HypergeometricDistribution dist = HypergeometricDistribution.of(n, m, s); +<a class="jxr_linenumber" name="L208" href="#L208">208</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p = dist.probability(s); +<a class="jxr_linenumber" name="L209" href="#L209">209</a> Assertions.assertEquals(1, p, () -> <span class="jxr_string">"p="</span> + p); +<a class="jxr_linenumber" name="L210" href="#L210">210</a> } +<a class="jxr_linenumber" name="L211" href="#L211">211</a> } +<a class="jxr_linenumber" name="L212" href="#L212">212</a> +<a class="jxr_linenumber" name="L213" href="#L213">213</a> @Test +<a class="jxr_linenumber" name="L214" href="#L214">214</a> <strong class="jxr_keyword">void</strong> testMath1021() { +<a class="jxr_linenumber" name="L215" href="#L215">215</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> N = 43130568; +<a class="jxr_linenumber" name="L216" href="#L216">216</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m = 42976365; +<a class="jxr_linenumber" name="L217" href="#L217">217</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = 50; +<a class="jxr_linenumber" name="L218" href="#L218">218</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler dist = +<a class="jxr_linenumber" name="L219" href="#L219">219</a> HypergeometricDistribution.of(N, m, n).createSampler(RandomSource.XO_SHI_RO_256_PP.create()); +<a class="jxr_linenumber" name="L220" href="#L220">220</a> +<a class="jxr_linenumber" name="L221" href="#L221">221</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < 100; i++) { +<a class="jxr_linenumber" name="L222" href="#L222">222</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> sample = dist.sample(); +<a class="jxr_linenumber" name="L223" href="#L223">223</a> Assertions.assertTrue(0 <= sample, () -> <span class="jxr_string">"sample="</span> + sample); +<a class="jxr_linenumber" name="L224" href="#L224">224</a> Assertions.assertTrue(sample <= n, () -> <span class="jxr_string">"sample="</span> + sample); +<a class="jxr_linenumber" name="L225" href="#L225">225</a> } +<a class="jxr_linenumber" name="L226" href="#L226">226</a> } +<a class="jxr_linenumber" name="L227" href="#L227">227</a> +<a class="jxr_linenumber" name="L228" href="#L228">228</a> @Test +<a class="jxr_linenumber" name="L229" href="#L229">229</a> <strong class="jxr_keyword">void</strong> testAdditionalCumulativeProbabilityHighPrecision() { +<a class="jxr_linenumber" name="L230" href="#L230">230</a> <em class="jxr_comment">// computed using R version 3.4.4</em> +<a class="jxr_linenumber" name="L231" href="#L231">231</a> testCumulativeProbabilityHighPrecision( +<a class="jxr_linenumber" name="L232" href="#L232">232</a> HypergeometricDistribution.of(500, 70, 300), +<a class="jxr_linenumber" name="L233" href="#L233">233</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {10, 8}, +<a class="jxr_linenumber" name="L234" href="#L234">234</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {2.4055720603264525e-17, 1.2848174992266236e-19}, +<a class="jxr_linenumber" name="L235" href="#L235">235</a> DoubleTolerances.relative(5e-14)); +<a class="jxr_linenumber" name="L236" href="#L236">236</a> } +<a class="jxr_linenumber" name="L237" href="#L237">237</a> +<a class="jxr_linenumber" name="L238" href="#L238">238</a> @Test +<a class="jxr_linenumber" name="L239" href="#L239">239</a> <strong class="jxr_keyword">void</strong> testAdditionalSurvivalProbabilityHighPrecision() { +<a class="jxr_linenumber" name="L240" href="#L240">240</a> <em class="jxr_comment">// computed using R version 3.4.4</em> +<a class="jxr_linenumber" name="L241" href="#L241">241</a> testSurvivalProbabilityHighPrecision( +<a class="jxr_linenumber" name="L242" href="#L242">242</a> HypergeometricDistribution.of(500, 70, 300), +<a class="jxr_linenumber" name="L243" href="#L243">243</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {68, 69}, +<a class="jxr_linenumber" name="L244" href="#L244">244</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {4.570379934029859e-16, 7.4187180434325268e-18}, +<a class="jxr_linenumber" name="L245" href="#L245">245</a> DoubleTolerances.relative(5e-14)); +<a class="jxr_linenumber" name="L246" href="#L246">246</a> } +<a class="jxr_linenumber" name="L247" href="#L247">247</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LaplaceDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LaplaceDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LaplaceDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,62 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>LaplaceDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/LaplaceDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <em class="jxr_javadoccomment"> * Test cases for {@link LaplaceDistribution}.</em> +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">class</strong> <a name="LaplaceDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/LaplaceDistributionTest.html#LaplaceDistributionTest">LaplaceDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L24" href="#L24">24</a> @Override +<a class="jxr_linenumber" name="L25" href="#L25">25</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> location = (Double) parameters[0]; +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = (Double) parameters[1]; +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <strong class="jxr_keyword">return</strong> LaplaceDistribution.of(location, scale); +<a class="jxr_linenumber" name="L29" href="#L29">29</a> } +<a class="jxr_linenumber" name="L30" href="#L30">30</a> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> +<a class="jxr_linenumber" name="L32" href="#L32">32</a> @Override +<a class="jxr_linenumber" name="L33" href="#L33">33</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L35" href="#L35">35</a> {0.0, 0.0}, +<a class="jxr_linenumber" name="L36" href="#L36">36</a> {0.0, -0.1} +<a class="jxr_linenumber" name="L37" href="#L37">37</a> }; +<a class="jxr_linenumber" name="L38" href="#L38">38</a> } +<a class="jxr_linenumber" name="L39" href="#L39">39</a> +<a class="jxr_linenumber" name="L40" href="#L40">40</a> @Override +<a class="jxr_linenumber" name="L41" href="#L41">41</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L42" href="#L42">42</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Location"</span>, <span class="jxr_string">"Scale"</span>}; +<a class="jxr_linenumber" name="L43" href="#L43">43</a> } +<a class="jxr_linenumber" name="L44" href="#L44">44</a> +<a class="jxr_linenumber" name="L45" href="#L45">45</a> @Override +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L47" href="#L47">47</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L48" href="#L48">48</a> } +<a class="jxr_linenumber" name="L49" href="#L49">49</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LevyDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LevyDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LevyDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,61 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>LevyDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/LevyDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <em class="jxr_javadoccomment"> * Test cases for {@link LevyDistribution}.</em> +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">class</strong> <a name="LevyDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/LevyDistributionTest.html#LevyDistributionTest">LevyDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L24" href="#L24">24</a> @Override +<a class="jxr_linenumber" name="L25" href="#L25">25</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> location = (Double) parameters[0]; +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = (Double) parameters[1]; +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <strong class="jxr_keyword">return</strong> LevyDistribution.of(location, scale); +<a class="jxr_linenumber" name="L29" href="#L29">29</a> } +<a class="jxr_linenumber" name="L30" href="#L30">30</a> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> @Override +<a class="jxr_linenumber" name="L32" href="#L32">32</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L34" href="#L34">34</a> {0.0, 0.0}, +<a class="jxr_linenumber" name="L35" href="#L35">35</a> {0.0, -0.1} +<a class="jxr_linenumber" name="L36" href="#L36">36</a> }; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> } +<a class="jxr_linenumber" name="L38" href="#L38">38</a> +<a class="jxr_linenumber" name="L39" href="#L39">39</a> @Override +<a class="jxr_linenumber" name="L40" href="#L40">40</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L41" href="#L41">41</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Location"</span>, <span class="jxr_string">"Scale"</span>}; +<a class="jxr_linenumber" name="L42" href="#L42">42</a> } +<a class="jxr_linenumber" name="L43" href="#L43">43</a> +<a class="jxr_linenumber" name="L44" href="#L44">44</a> @Override +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L47" href="#L47">47</a> } +<a class="jxr_linenumber" name="L48" href="#L48">48</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogNormalDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogNormalDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogNormalDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,114 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>LogNormalDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/LogNormalDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> +<a class="jxr_linenumber" name="L18" href="#L18">18</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L19" href="#L19">19</a> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <em class="jxr_javadoccomment"> * Test cases for {@link LogNormalDistribution}.</em> +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">class</strong> <a name="LogNormalDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/LogNormalDistributionTest.html#LogNormalDistributionTest">LogNormalDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L28" href="#L28">28</a> @Override +<a class="jxr_linenumber" name="L29" href="#L29">29</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> mu = (Double) parameters[0]; +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> sigma = (Double) parameters[1]; +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <strong class="jxr_keyword">return</strong> LogNormalDistribution.of(mu, sigma); +<a class="jxr_linenumber" name="L33" href="#L33">33</a> } +<a class="jxr_linenumber" name="L34" href="#L34">34</a> +<a class="jxr_linenumber" name="L35" href="#L35">35</a> @Override +<a class="jxr_linenumber" name="L36" href="#L36">36</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L38" href="#L38">38</a> {0.0, 0.0}, +<a class="jxr_linenumber" name="L39" href="#L39">39</a> {0.0, -0.1} +<a class="jxr_linenumber" name="L40" href="#L40">40</a> }; +<a class="jxr_linenumber" name="L41" href="#L41">41</a> } +<a class="jxr_linenumber" name="L42" href="#L42">42</a> +<a class="jxr_linenumber" name="L43" href="#L43">43</a> @Override +<a class="jxr_linenumber" name="L44" href="#L44">44</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Mu"</span>, <span class="jxr_string">"Sigma"</span>}; +<a class="jxr_linenumber" name="L46" href="#L46">46</a> } +<a class="jxr_linenumber" name="L47" href="#L47">47</a> +<a class="jxr_linenumber" name="L48" href="#L48">48</a> @Override +<a class="jxr_linenumber" name="L49" href="#L49">49</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L50" href="#L50">50</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L51" href="#L51">51</a> } +<a class="jxr_linenumber" name="L52" href="#L52">52</a> +<a class="jxr_linenumber" name="L53" href="#L53">53</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L54" href="#L54">54</a> +<a class="jxr_linenumber" name="L55" href="#L55">55</a> @Test +<a class="jxr_linenumber" name="L56" href="#L56">56</a> <strong class="jxr_keyword">void</strong> testCumulativeProbabilityExtremes() { +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <em class="jxr_comment">// Use a small shape parameter so that we can exceed 40 * shape</em> +<a class="jxr_linenumber" name="L58" href="#L58">58</a> testCumulativeProbability(LogNormalDistribution.of(1, 0.0001), +<a class="jxr_linenumber" name="L59" href="#L59">59</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {0.5, 10}, +<a class="jxr_linenumber" name="L60" href="#L60">60</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {0, 1.0}, +<a class="jxr_linenumber" name="L61" href="#L61">61</a> DoubleTolerances.equals()); +<a class="jxr_linenumber" name="L62" href="#L62">62</a> } +<a class="jxr_linenumber" name="L63" href="#L63">63</a> +<a class="jxr_linenumber" name="L64" href="#L64">64</a> @Test +<a class="jxr_linenumber" name="L65" href="#L65">65</a> <strong class="jxr_keyword">void</strong> testSurvivalProbabilityExtremes() { +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <em class="jxr_comment">// Use a small shape parameter so that we can exceed 40 * shape</em> +<a class="jxr_linenumber" name="L67" href="#L67">67</a> testSurvivalProbability(LogNormalDistribution.of(1, 0.0001), +<a class="jxr_linenumber" name="L68" href="#L68">68</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {0.5, 10}, +<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {1.0, 0.0}, +<a class="jxr_linenumber" name="L70" href="#L70">70</a> DoubleTolerances.equals()); +<a class="jxr_linenumber" name="L71" href="#L71">71</a> } +<a class="jxr_linenumber" name="L72" href="#L72">72</a> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <em class="jxr_javadoccomment"> * Check to make sure top-coding of extreme values works correctly.</em> +<a class="jxr_linenumber" name="L75" href="#L75">75</a> <em class="jxr_javadoccomment"> * Verifies fixes for JIRA MATH-167, MATH-414</em> +<a class="jxr_linenumber" name="L76" href="#L76">76</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L77" href="#L77">77</a> @Test +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <strong class="jxr_keyword">void</strong> testExtremeValues() { +<a class="jxr_linenumber" name="L79" href="#L79">79</a> <strong class="jxr_keyword">final</strong> LogNormalDistribution dist = LogNormalDistribution.of(0, 1); +<a class="jxr_linenumber" name="L80" href="#L80">80</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < 1e5; i++) { <em class="jxr_comment">// make sure no convergence exception</em> +<a class="jxr_linenumber" name="L81" href="#L81">81</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> upperTail = dist.cumulativeProbability(i); +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <strong class="jxr_keyword">if</strong> (i <= 72) { <em class="jxr_comment">// make sure not top-coded</em> +<a class="jxr_linenumber" name="L83" href="#L83">83</a> Assertions.assertTrue(upperTail < 1.0d); +<a class="jxr_linenumber" name="L84" href="#L84">84</a> } <strong class="jxr_keyword">else</strong> { <em class="jxr_comment">// make sure top coding not reversed</em> +<a class="jxr_linenumber" name="L85" href="#L85">85</a> Assertions.assertTrue(upperTail > 0.99999); +<a class="jxr_linenumber" name="L86" href="#L86">86</a> } +<a class="jxr_linenumber" name="L87" href="#L87">87</a> } +<a class="jxr_linenumber" name="L88" href="#L88">88</a> +<a class="jxr_linenumber" name="L89" href="#L89">89</a> Assertions.assertEquals(1, dist.cumulativeProbability(Double.MAX_VALUE)); +<a class="jxr_linenumber" name="L90" href="#L90">90</a> Assertions.assertEquals(0, dist.cumulativeProbability(-Double.MAX_VALUE)); +<a class="jxr_linenumber" name="L91" href="#L91">91</a> Assertions.assertEquals(1, dist.cumulativeProbability(Double.POSITIVE_INFINITY)); +<a class="jxr_linenumber" name="L92" href="#L92">92</a> Assertions.assertEquals(0, dist.cumulativeProbability(Double.NEGATIVE_INFINITY)); +<a class="jxr_linenumber" name="L93" href="#L93">93</a> } +<a class="jxr_linenumber" name="L94" href="#L94">94</a> +<a class="jxr_linenumber" name="L95" href="#L95">95</a> @Test +<a class="jxr_linenumber" name="L96" href="#L96">96</a> <strong class="jxr_keyword">void</strong> testTinyVariance() { +<a class="jxr_linenumber" name="L97" href="#L97">97</a> <strong class="jxr_keyword">final</strong> LogNormalDistribution dist = LogNormalDistribution.of(0, 1e-9); +<a class="jxr_linenumber" name="L98" href="#L98">98</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> t = dist.getVariance(); +<a class="jxr_linenumber" name="L99" href="#L99">99</a> Assertions.assertEquals(1e-18, t, 1e-20); +<a class="jxr_linenumber" name="L100" href="#L100">100</a> } +<a class="jxr_linenumber" name="L101" href="#L101">101</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogisticDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogisticDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/LogisticDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,108 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>LogisticDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/LogisticDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <em class="jxr_javadoccomment"> * Test cases for {@link LogisticDistribution}.</em> +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">class</strong> <a name="LogisticDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/LogisticDistributionTest.html#LogisticDistributionTest">LogisticDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L27" href="#L27">27</a> @Override +<a class="jxr_linenumber" name="L28" href="#L28">28</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> location = (Double) parameters[0]; +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = (Double) parameters[1]; +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <strong class="jxr_keyword">return</strong> LogisticDistribution.of(location, scale); +<a class="jxr_linenumber" name="L32" href="#L32">32</a> } +<a class="jxr_linenumber" name="L33" href="#L33">33</a> +<a class="jxr_linenumber" name="L34" href="#L34">34</a> @Override +<a class="jxr_linenumber" name="L35" href="#L35">35</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L37" href="#L37">37</a> {0.0, 0.0}, +<a class="jxr_linenumber" name="L38" href="#L38">38</a> {0.0, -0.1} +<a class="jxr_linenumber" name="L39" href="#L39">39</a> }; +<a class="jxr_linenumber" name="L40" href="#L40">40</a> } +<a class="jxr_linenumber" name="L41" href="#L41">41</a> +<a class="jxr_linenumber" name="L42" href="#L42">42</a> @Override +<a class="jxr_linenumber" name="L43" href="#L43">43</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L44" href="#L44">44</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Location"</span>, <span class="jxr_string">"Scale"</span>}; +<a class="jxr_linenumber" name="L45" href="#L45">45</a> } +<a class="jxr_linenumber" name="L46" href="#L46">46</a> +<a class="jxr_linenumber" name="L47" href="#L47">47</a> @Override +<a class="jxr_linenumber" name="L48" href="#L48">48</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L49" href="#L49">49</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L50" href="#L50">50</a> } +<a class="jxr_linenumber" name="L51" href="#L51">51</a> +<a class="jxr_linenumber" name="L52" href="#L52">52</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L53" href="#L53">53</a> +<a class="jxr_linenumber" name="L54" href="#L54">54</a> @Test +<a class="jxr_linenumber" name="L55" href="#L55">55</a> <strong class="jxr_keyword">void</strong> testExtremeDensity() { +<a class="jxr_linenumber" name="L56" href="#L56">56</a> <strong class="jxr_keyword">final</strong> LogisticDistribution dist = LogisticDistribution.of(0, 1.0); +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <em class="jxr_comment">// Direct density (with scale = 1):</em> +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <em class="jxr_comment">// exp(-x) / (1 + exp(-x))^2</em> +<a class="jxr_linenumber" name="L59" href="#L59">59</a> <em class="jxr_comment">// As x -> large negative then exp(-x) will overflow and a simple</em> +<a class="jxr_linenumber" name="L60" href="#L60">60</a> <em class="jxr_comment">// computation will be incorrect</em> +<a class="jxr_linenumber" name="L61" href="#L61">61</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x0 = -710; +<a class="jxr_linenumber" name="L62" href="#L62">62</a> Assertions.assertEquals(Double.POSITIVE_INFINITY, Math.exp(-x0)); +<a class="jxr_linenumber" name="L63" href="#L63">63</a> Assertions.assertEquals(Double.NaN, Math.exp(-x0) / Math.pow(1 + Math.exp(-x0), 2.0)); +<a class="jxr_linenumber" name="L64" href="#L64">64</a> +<a class="jxr_linenumber" name="L65" href="#L65">65</a> <em class="jxr_comment">// Computed using scipy.stats logistic.</em> +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <em class="jxr_comment">// These values of exp(x) will create overflow.</em> +<a class="jxr_linenumber" name="L67" href="#L67">67</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] x = {710, 720, 730, 740, 750}; +<a class="jxr_linenumber" name="L68" href="#L68">68</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] values = {4.47628622567513e-309, +<a class="jxr_linenumber" name="L69" href="#L69">69</a> 2.03223080241836e-313, 9.22631526816382e-318, +<a class="jxr_linenumber" name="L70" href="#L70">70</a> 4.19955798965060e-322, 0.00000000000000e+000}; +<a class="jxr_linenumber" name="L71" href="#L71">71</a> testDensity(dist, x, values, DoubleTolerances.equals()); +<a class="jxr_linenumber" name="L72" href="#L72">72</a> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <em class="jxr_comment">// Test symmetry.</em> +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xi : x) { +<a class="jxr_linenumber" name="L75" href="#L75">75</a> Assertions.assertEquals(dist.density(xi), dist.density(-xi), () -> Double.toString(xi)); +<a class="jxr_linenumber" name="L76" href="#L76">76</a> } +<a class="jxr_linenumber" name="L77" href="#L77">77</a> } +<a class="jxr_linenumber" name="L78" href="#L78">78</a> +<a class="jxr_linenumber" name="L79" href="#L79">79</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L80" href="#L80">80</a> <em class="jxr_javadoccomment"> * Test a value for log density when the density computation is zero.</em> +<a class="jxr_linenumber" name="L81" href="#L81">81</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L82" href="#L82">82</a> @Test +<a class="jxr_linenumber" name="L83" href="#L83">83</a> <strong class="jxr_keyword">void</strong> testExtremeLogDensity() { +<a class="jxr_linenumber" name="L84" href="#L84">84</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = 2.5; +<a class="jxr_linenumber" name="L85" href="#L85">85</a> <strong class="jxr_keyword">final</strong> LogisticDistribution dist = LogisticDistribution.of(0, scale); +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <em class="jxr_comment">// Direct density (with scale = s):</em> +<a class="jxr_linenumber" name="L87" href="#L87">87</a> <em class="jxr_comment">// exp(-x / s) / (1 + exp(-x / s))^2</em> +<a class="jxr_linenumber" name="L88" href="#L88">88</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = 1e160; +<a class="jxr_linenumber" name="L89" href="#L89">89</a> Assertions.assertEquals(0.0, dist.density(x)); +<a class="jxr_linenumber" name="L90" href="#L90">90</a> <em class="jxr_comment">// Log computation</em> +<a class="jxr_linenumber" name="L91" href="#L91">91</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = -x / scale - 2 * Math.log1p(Math.exp(-x / scale)); +<a class="jxr_linenumber" name="L92" href="#L92">92</a> Assertions.assertNotEquals(Double.NEGATIVE_INFINITY, expected, <span class="jxr_string">"Density is zero but log density should not be -infinity"</span>); +<a class="jxr_linenumber" name="L93" href="#L93">93</a> Assertions.assertEquals(expected, dist.logDensity(x), Math.abs(expected) * 1e-15); +<a class="jxr_linenumber" name="L94" href="#L94">94</a> } +<a class="jxr_linenumber" name="L95" href="#L95">95</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html>
