Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/UniformDiscreteDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/UniformDiscreteDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/UniformDiscreteDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,280 @@ +<!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>UniformDiscreteDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/UniformDiscreteDistributionTest.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.math3.util.MathArrays; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments; +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.CsvSource; +<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 UniformDiscreteDistribution}.</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="UniformDiscreteDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/UniformDiscreteDistributionTest.html#UniformDiscreteDistributionTest">UniformDiscreteDistributionTest</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> lower = (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> upper = (Integer) parameters[1]; +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">return</strong> UniformDiscreteDistribution.of(lower, upper); +<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> 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> <em class="jxr_comment">// MATH-1141</em> +<a class="jxr_linenumber" name="L46" href="#L46">46</a> {1, 0}, +<a class="jxr_linenumber" name="L47" href="#L47">47</a> {3, 2}, +<a class="jxr_linenumber" name="L48" href="#L48">48</a> }; +<a class="jxr_linenumber" name="L49" href="#L49">49</a> } +<a class="jxr_linenumber" name="L50" href="#L50">50</a> +<a class="jxr_linenumber" name="L51" href="#L51">51</a> @Override +<a class="jxr_linenumber" name="L52" href="#L52">52</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L53" href="#L53">53</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"SupportLowerBound"</span>, <span class="jxr_string">"SupportUpperBound"</span>}; +<a class="jxr_linenumber" name="L54" href="#L54">54</a> } +<a class="jxr_linenumber" name="L55" href="#L55">55</a> +<a class="jxr_linenumber" name="L56" href="#L56">56</a> @Override +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <em class="jxr_comment">// Tolerance is 4.440892098500626E-16</em> +<a class="jxr_linenumber" name="L59" href="#L59">59</a> <strong class="jxr_keyword">return</strong> 2 * RELATIVE_EPS; +<a class="jxr_linenumber" name="L60" href="#L60">60</a> } +<a class="jxr_linenumber" name="L61" href="#L61">61</a> +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L63" href="#L63">63</a> +<a class="jxr_linenumber" name="L64" href="#L64">64</a> @ParameterizedTest +<a class="jxr_linenumber" name="L65" href="#L65">65</a> @MethodSource +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <strong class="jxr_keyword">void</strong> testAdditionalMoments(<strong class="jxr_keyword">int</strong> lower, <strong class="jxr_keyword">int</strong> upper, <strong class="jxr_keyword">double</strong> mean, <strong class="jxr_keyword">double</strong> variance) { +<a class="jxr_linenumber" name="L67" href="#L67">67</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L68" href="#L68">68</a> testMoments(dist, mean, variance, DoubleTolerances.equals()); +<a class="jxr_linenumber" name="L69" href="#L69">69</a> } +<a class="jxr_linenumber" name="L70" href="#L70">70</a> +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionalMoments() { +<a class="jxr_linenumber" name="L72" href="#L72">72</a> <strong class="jxr_keyword">return</strong> Stream.of( +<a class="jxr_linenumber" name="L73" href="#L73">73</a> Arguments.of(0, 5, 2.5, 35 / 12.0), +<a class="jxr_linenumber" name="L74" href="#L74">74</a> Arguments.of(0, 1, 0.5, 3 / 12.0) +<a class="jxr_linenumber" name="L75" href="#L75">75</a> ); +<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> <em class="jxr_comment">// MATH-1396</em> +<a class="jxr_linenumber" name="L79" href="#L79">79</a> @Test +<a class="jxr_linenumber" name="L80" href="#L80">80</a> <strong class="jxr_keyword">void</strong> testLargeRangeSubtractionOverflow() { +<a class="jxr_linenumber" name="L81" href="#L81">81</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> hi = Integer.MAX_VALUE / 2 + 10; +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> lower = -hi; +<a class="jxr_linenumber" name="L83" href="#L83">83</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> upper = hi - 1; +<a class="jxr_linenumber" name="L84" href="#L84">84</a> +<a class="jxr_linenumber" name="L85" href="#L85">85</a> <em class="jxr_comment">// range = upper - lower + 1 would overflow</em> +<a class="jxr_linenumber" name="L86" href="#L86">86</a> Assertions.assertTrue(upper - lower < 0); +<a class="jxr_linenumber" name="L87" href="#L87">87</a> +<a class="jxr_linenumber" name="L88" href="#L88">88</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L89" href="#L89">89</a> +<a class="jxr_linenumber" name="L90" href="#L90">90</a> Assertions.assertEquals(0.5 / hi, dist.probability(123456)); +<a class="jxr_linenumber" name="L91" href="#L91">91</a> Assertions.assertEquals(0.5, dist.cumulativeProbability(-1)); +<a class="jxr_linenumber" name="L92" href="#L92">92</a> +<a class="jxr_linenumber" name="L93" href="#L93">93</a> Assertions.assertEquals((Math.pow(2d * hi, 2) - 1) / 12, dist.getVariance()); +<a class="jxr_linenumber" name="L94" href="#L94">94</a> } +<a class="jxr_linenumber" name="L95" href="#L95">95</a> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> <em class="jxr_comment">// MATH-1396</em> +<a class="jxr_linenumber" name="L97" href="#L97">97</a> @Test +<a class="jxr_linenumber" name="L98" href="#L98">98</a> <strong class="jxr_keyword">void</strong> testLargeRangeAdditionOverflow() { +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> hi = Integer.MAX_VALUE / 2 + 10; +<a class="jxr_linenumber" name="L100" href="#L100">100</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> lower = hi - 1; +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> upper = hi + 1; +<a class="jxr_linenumber" name="L102" href="#L102">102</a> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <em class="jxr_comment">// mean = (lower + upper) / 2 would overflow</em> +<a class="jxr_linenumber" name="L104" href="#L104">104</a> Assertions.assertTrue(lower + upper < 0); +<a class="jxr_linenumber" name="L105" href="#L105">105</a> +<a class="jxr_linenumber" name="L106" href="#L106">106</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L107" href="#L107">107</a> +<a class="jxr_linenumber" name="L108" href="#L108">108</a> Assertions.assertEquals(1d / 3d, dist.probability(hi)); +<a class="jxr_linenumber" name="L109" href="#L109">109</a> Assertions.assertEquals(2d / 3d, dist.cumulativeProbability(hi)); +<a class="jxr_linenumber" name="L110" href="#L110">110</a> +<a class="jxr_linenumber" name="L111" href="#L111">111</a> Assertions.assertEquals(hi, dist.getMean()); +<a class="jxr_linenumber" name="L112" href="#L112">112</a> } +<a class="jxr_linenumber" name="L113" href="#L113">113</a> +<a class="jxr_linenumber" name="L114" href="#L114">114</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L115" href="#L115">115</a> <em class="jxr_javadoccomment"> * Test the inverse CDF returns the correct x from the CDF result.</em> +<a class="jxr_linenumber" name="L116" href="#L116">116</a> <em class="jxr_javadoccomment"> * Test cases created to generate rounding errors on the inversion.</em> +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L118" href="#L118">118</a> @ParameterizedTest +<a class="jxr_linenumber" name="L119" href="#L119">119</a> @CsvSource(value = { +<a class="jxr_linenumber" name="L120" href="#L120">120</a> <em class="jxr_comment">// Extreme bounds</em> +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <span class="jxr_string">"-2147483648, -2147483648"</span>, +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <span class="jxr_string">"-2147483648, -2147483647"</span>, +<a class="jxr_linenumber" name="L123" href="#L123">123</a> <span class="jxr_string">"-2147483648, -2147483646"</span>, +<a class="jxr_linenumber" name="L124" href="#L124">124</a> <span class="jxr_string">"-2147483648, -2147483638"</span>, +<a class="jxr_linenumber" name="L125" href="#L125">125</a> <span class="jxr_string">"2147483647, 2147483647"</span>, +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <span class="jxr_string">"2147483646, 2147483647"</span>, +<a class="jxr_linenumber" name="L127" href="#L127">127</a> <span class="jxr_string">"2147483645, 2147483647"</span>, +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <span class="jxr_string">"2147483637, 2147483647"</span>, +<a class="jxr_linenumber" name="L129" href="#L129">129</a> <em class="jxr_comment">// icdf(cdf(x)) requires rounding up</em> +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <span class="jxr_string">"3, 40"</span>, +<a class="jxr_linenumber" name="L131" href="#L131">131</a> <span class="jxr_string">"71, 201"</span>, +<a class="jxr_linenumber" name="L132" href="#L132">132</a> <span class="jxr_string">"223, 267"</span>, +<a class="jxr_linenumber" name="L133" href="#L133">133</a> <span class="jxr_string">"45, 125"</span>, +<a class="jxr_linenumber" name="L134" href="#L134">134</a> <span class="jxr_string">"53, 81"</span>, +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_comment">// icdf(cdf(x)) requires rounding down</em> +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <span class="jxr_string">"48, 247"</span>, +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <span class="jxr_string">"141, 222"</span>, +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <span class="jxr_string">"106, 223"</span>, +<a class="jxr_linenumber" name="L139" href="#L139">139</a> <span class="jxr_string">"156, 201"</span>, +<a class="jxr_linenumber" name="L140" href="#L140">140</a> <span class="jxr_string">"86, 265"</span>, +<a class="jxr_linenumber" name="L141" href="#L141">141</a> }) +<a class="jxr_linenumber" name="L142" href="#L142">142</a> <strong class="jxr_keyword">void</strong> testInverseCDF(<strong class="jxr_keyword">int</strong> lower, <strong class="jxr_keyword">int</strong> upper) { +<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L144" href="#L144">144</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] x = MathArrays.sequence(upper - lower, lower, 1); +<a class="jxr_linenumber" name="L145" href="#L145">145</a> testCumulativeProbabilityInverseMapping(dist, x); +<a class="jxr_linenumber" name="L146" href="#L146">146</a> } +<a class="jxr_linenumber" name="L147" href="#L147">147</a> +<a class="jxr_linenumber" name="L148" href="#L148">148</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L149" href="#L149">149</a> <em class="jxr_javadoccomment"> * Test the inverse SF returns the correct x from the SF result.</em> +<a class="jxr_linenumber" name="L150" href="#L150">150</a> <em class="jxr_javadoccomment"> * Test cases created to generate rounding errors on the inversion.</em> +<a class="jxr_linenumber" name="L151" href="#L151">151</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L152" href="#L152">152</a> @ParameterizedTest +<a class="jxr_linenumber" name="L153" href="#L153">153</a> @CsvSource(value = { +<a class="jxr_linenumber" name="L154" href="#L154">154</a> <em class="jxr_comment">// Extreme bounds</em> +<a class="jxr_linenumber" name="L155" href="#L155">155</a> <span class="jxr_string">"-2147483648, -2147483648"</span>, +<a class="jxr_linenumber" name="L156" href="#L156">156</a> <span class="jxr_string">"-2147483648, -2147483647"</span>, +<a class="jxr_linenumber" name="L157" href="#L157">157</a> <span class="jxr_string">"-2147483648, -2147483646"</span>, +<a class="jxr_linenumber" name="L158" href="#L158">158</a> <span class="jxr_string">"-2147483648, -2147483638"</span>, +<a class="jxr_linenumber" name="L159" href="#L159">159</a> <span class="jxr_string">"2147483647, 2147483647"</span>, +<a class="jxr_linenumber" name="L160" href="#L160">160</a> <span class="jxr_string">"2147483646, 2147483647"</span>, +<a class="jxr_linenumber" name="L161" href="#L161">161</a> <span class="jxr_string">"2147483645, 2147483647"</span>, +<a class="jxr_linenumber" name="L162" href="#L162">162</a> <span class="jxr_string">"2147483637, 2147483647"</span>, +<a class="jxr_linenumber" name="L163" href="#L163">163</a> <em class="jxr_comment">// isf(sf(x)) requires rounding up</em> +<a class="jxr_linenumber" name="L164" href="#L164">164</a> <span class="jxr_string">"52, 91"</span>, +<a class="jxr_linenumber" name="L165" href="#L165">165</a> <span class="jxr_string">"81, 106"</span>, +<a class="jxr_linenumber" name="L166" href="#L166">166</a> <span class="jxr_string">"79, 268"</span>, +<a class="jxr_linenumber" name="L167" href="#L167">167</a> <span class="jxr_string">"54, 249"</span>, +<a class="jxr_linenumber" name="L168" href="#L168">168</a> <span class="jxr_string">"189, 267"</span>, +<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_comment">// isf(sf(x)) requires rounding down</em> +<a class="jxr_linenumber" name="L170" href="#L170">170</a> <span class="jxr_string">"105, 279"</span>, +<a class="jxr_linenumber" name="L171" href="#L171">171</a> <span class="jxr_string">"42, 261"</span>, +<a class="jxr_linenumber" name="L172" href="#L172">172</a> <span class="jxr_string">"37, 133"</span>, +<a class="jxr_linenumber" name="L173" href="#L173">173</a> <span class="jxr_string">"59, 214"</span>, +<a class="jxr_linenumber" name="L174" href="#L174">174</a> <span class="jxr_string">"33, 118"</span>, +<a class="jxr_linenumber" name="L175" href="#L175">175</a> }) +<a class="jxr_linenumber" name="L176" href="#L176">176</a> <strong class="jxr_keyword">void</strong> testInverseSF(<strong class="jxr_keyword">int</strong> lower, <strong class="jxr_keyword">int</strong> upper) { +<a class="jxr_linenumber" name="L177" href="#L177">177</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L178" href="#L178">178</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] x = MathArrays.sequence(upper - lower, lower, 1); +<a class="jxr_linenumber" name="L179" href="#L179">179</a> testSurvivalProbabilityInverseMapping(dist, x); +<a class="jxr_linenumber" name="L180" href="#L180">180</a> } +<a class="jxr_linenumber" name="L181" href="#L181">181</a> +<a class="jxr_linenumber" name="L182" href="#L182">182</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L183" href="#L183">183</a> <em class="jxr_javadoccomment"> * Test the probability in a range uses the exact computation of</em> +<a class="jxr_linenumber" name="L184" href="#L184">184</a> <em class="jxr_javadoccomment"> * {@code (x1 - x0) / (upper - lower + 1)} assuming x0 and x1 are within [lower, upper].</em> +<a class="jxr_linenumber" name="L185" href="#L185">185</a> <em class="jxr_javadoccomment"> * This test will fail if the distribution uses the default implementation in</em> +<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment"> * {@link AbstractDiscreteDistribution}.</em> +<a class="jxr_linenumber" name="L187" href="#L187">187</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L188" href="#L188">188</a> @ParameterizedTest +<a class="jxr_linenumber" name="L189" href="#L189">189</a> @CsvSource(value = { +<a class="jxr_linenumber" name="L190" href="#L190">190</a> <em class="jxr_comment">// Extreme bounds</em> +<a class="jxr_linenumber" name="L191" href="#L191">191</a> <span class="jxr_string">"-2147483648, -2147483648"</span>, +<a class="jxr_linenumber" name="L192" href="#L192">192</a> <span class="jxr_string">"-2147483648, -2147483647"</span>, +<a class="jxr_linenumber" name="L193" href="#L193">193</a> <span class="jxr_string">"-2147483648, -2147483646"</span>, +<a class="jxr_linenumber" name="L194" href="#L194">194</a> <span class="jxr_string">"-2147483648, -2147483638"</span>, +<a class="jxr_linenumber" name="L195" href="#L195">195</a> <span class="jxr_string">"2147483647, 2147483647"</span>, +<a class="jxr_linenumber" name="L196" href="#L196">196</a> <span class="jxr_string">"2147483646, 2147483647"</span>, +<a class="jxr_linenumber" name="L197" href="#L197">197</a> <span class="jxr_string">"2147483645, 2147483647"</span>, +<a class="jxr_linenumber" name="L198" href="#L198">198</a> <span class="jxr_string">"2147483637, 2147483647"</span>, +<a class="jxr_linenumber" name="L199" href="#L199">199</a> <em class="jxr_comment">// Range is a prime number</em> +<a class="jxr_linenumber" name="L200" href="#L200">200</a> <span class="jxr_string">"-10, 2"</span>, <em class="jxr_comment">// 13</em> +<a class="jxr_linenumber" name="L201" href="#L201">201</a> <span class="jxr_string">"10, 16"</span>, <em class="jxr_comment">// 7</em> +<a class="jxr_linenumber" name="L202" href="#L202">202</a> <span class="jxr_string">"-20, -10"</span>, <em class="jxr_comment">// 11</em> +<a class="jxr_linenumber" name="L203" href="#L203">203</a> <em class="jxr_comment">// Range is even</em> +<a class="jxr_linenumber" name="L204" href="#L204">204</a> <span class="jxr_string">"-10, 3"</span>, <em class="jxr_comment">// 14</em> +<a class="jxr_linenumber" name="L205" href="#L205">205</a> <span class="jxr_string">"10, 17"</span>, <em class="jxr_comment">// 8</em> +<a class="jxr_linenumber" name="L206" href="#L206">206</a> <span class="jxr_string">"-20, -9"</span>, <em class="jxr_comment">// 12</em> +<a class="jxr_linenumber" name="L207" href="#L207">207</a> <em class="jxr_comment">// Large range</em> +<a class="jxr_linenumber" name="L208" href="#L208">208</a> <span class="jxr_string">"-2147483648, 2147483647"</span>, +<a class="jxr_linenumber" name="L209" href="#L209">209</a> <span class="jxr_string">"-2147483648, 1263781682"</span>, +<a class="jxr_linenumber" name="L210" href="#L210">210</a> <span class="jxr_string">"-2147483648, 1781682"</span>, +<a class="jxr_linenumber" name="L211" href="#L211">211</a> <span class="jxr_string">"-2147483648, -231781682"</span>, +<a class="jxr_linenumber" name="L212" href="#L212">212</a> <span class="jxr_string">"-1324234584, 2147483647"</span>, +<a class="jxr_linenumber" name="L213" href="#L213">213</a> <span class="jxr_string">"-324234584, 2147483647"</span>, +<a class="jxr_linenumber" name="L214" href="#L214">214</a> <span class="jxr_string">"6234584, 2147483647"</span>, +<a class="jxr_linenumber" name="L215" href="#L215">215</a> <span class="jxr_string">"-1256362376, 125637"</span>, +<a class="jxr_linenumber" name="L216" href="#L216">216</a> <span class="jxr_string">"-62378468, 1325657374"</span>, +<a class="jxr_linenumber" name="L217" href="#L217">217</a> }) +<a class="jxr_linenumber" name="L218" href="#L218">218</a> <strong class="jxr_keyword">void</strong> testProbabilityRange(<strong class="jxr_keyword">int</strong> lower, <strong class="jxr_keyword">int</strong> upper) { +<a class="jxr_linenumber" name="L219" href="#L219">219</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(lower, upper); +<a class="jxr_linenumber" name="L220" href="#L220">220</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> r = (<strong class="jxr_keyword">double</strong>) upper - lower + 1; +<a class="jxr_linenumber" name="L221" href="#L221">221</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> stride = r < 20 ? 1 : (<strong class="jxr_keyword">long</strong>) (r / 20); +<a class="jxr_linenumber" name="L222" href="#L222">222</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">long</strong> x0 = lower; x0 <= upper; x0 += stride) { +<a class="jxr_linenumber" name="L223" href="#L223">223</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">long</strong> x1 = x0; x1 <= upper; x1 += stride) { +<a class="jxr_linenumber" name="L224" href="#L224">224</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p = (x1 - x0) / r; +<a class="jxr_linenumber" name="L225" href="#L225">225</a> Assertions.assertEquals(p, dist.probability((<strong class="jxr_keyword">int</strong>) x0, (<strong class="jxr_keyword">int</strong>) x1)); +<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> } +<a class="jxr_linenumber" name="L229" href="#L229">229</a> +<a class="jxr_linenumber" name="L230" href="#L230">230</a> @Test +<a class="jxr_linenumber" name="L231" href="#L231">231</a> <strong class="jxr_keyword">void</strong> testProbabilityRangeEdgeCases() { +<a class="jxr_linenumber" name="L232" href="#L232">232</a> <strong class="jxr_keyword">final</strong> UniformDiscreteDistribution dist = UniformDiscreteDistribution.of(3, 5); +<a class="jxr_linenumber" name="L233" href="#L233">233</a> +<a class="jxr_linenumber" name="L234" href="#L234">234</a> Assertions.assertThrows(DistributionException.<strong class="jxr_keyword">class</strong>, () -> dist.probability(4, 3)); +<a class="jxr_linenumber" name="L235" href="#L235">235</a> +<a class="jxr_linenumber" name="L236" href="#L236">236</a> <em class="jxr_comment">// x0 >= upper</em> +<a class="jxr_linenumber" name="L237" href="#L237">237</a> Assertions.assertEquals(0, dist.probability(5, 6)); +<a class="jxr_linenumber" name="L238" href="#L238">238</a> Assertions.assertEquals(0, dist.probability(15, 16)); +<a class="jxr_linenumber" name="L239" href="#L239">239</a> <em class="jxr_comment">// x1 < lower</em> +<a class="jxr_linenumber" name="L240" href="#L240">240</a> Assertions.assertEquals(0, dist.probability(-3, 1)); +<a class="jxr_linenumber" name="L241" href="#L241">241</a> +<a class="jxr_linenumber" name="L242" href="#L242">242</a> <em class="jxr_comment">// x0 == x1</em> +<a class="jxr_linenumber" name="L243" href="#L243">243</a> Assertions.assertEquals(0, dist.probability(3, 3)); +<a class="jxr_linenumber" name="L244" href="#L244">244</a> Assertions.assertEquals(0, dist.probability(4, 4)); +<a class="jxr_linenumber" name="L245" href="#L245">245</a> Assertions.assertEquals(0, dist.probability(5, 5)); +<a class="jxr_linenumber" name="L246" href="#L246">246</a> Assertions.assertEquals(0, dist.probability(6, 6)); +<a class="jxr_linenumber" name="L247" href="#L247">247</a> +<a class="jxr_linenumber" name="L248" href="#L248">248</a> <em class="jxr_comment">// x0+1 == x1</em> +<a class="jxr_linenumber" name="L249" href="#L249">249</a> Assertions.assertEquals(1.0 / 3, dist.probability(3, 4)); +<a class="jxr_linenumber" name="L250" href="#L250">250</a> Assertions.assertEquals(1.0 / 3, dist.probability(4, 5)); +<a class="jxr_linenumber" name="L251" href="#L251">251</a> +<a class="jxr_linenumber" name="L252" href="#L252">252</a> <em class="jxr_comment">// x1 > upper</em> +<a class="jxr_linenumber" name="L253" href="#L253">253</a> Assertions.assertEquals(1, dist.probability(2, 6)); +<a class="jxr_linenumber" name="L254" href="#L254">254</a> Assertions.assertEquals(2.0 / 3, dist.probability(3, 6)); +<a class="jxr_linenumber" name="L255" href="#L255">255</a> Assertions.assertEquals(1.0 / 3, dist.probability(4, 6)); +<a class="jxr_linenumber" name="L256" href="#L256">256</a> Assertions.assertEquals(0, dist.probability(5, 6)); +<a class="jxr_linenumber" name="L257" href="#L257">257</a> +<a class="jxr_linenumber" name="L258" href="#L258">258</a> <em class="jxr_comment">// x0 < lower</em> +<a class="jxr_linenumber" name="L259" href="#L259">259</a> Assertions.assertEquals(0, dist.probability(-2, 2)); +<a class="jxr_linenumber" name="L260" href="#L260">260</a> Assertions.assertEquals(1.0 / 3, dist.probability(-2, 3)); +<a class="jxr_linenumber" name="L261" href="#L261">261</a> Assertions.assertEquals(2.0 / 3, dist.probability(-2, 4)); +<a class="jxr_linenumber" name="L262" href="#L262">262</a> Assertions.assertEquals(1.0, dist.probability(-2, 5)); +<a class="jxr_linenumber" name="L263" href="#L263">263</a> +<a class="jxr_linenumber" name="L264" href="#L264">264</a> <em class="jxr_comment">// x1 > upper && x0 < lower</em> +<a class="jxr_linenumber" name="L265" href="#L265">265</a> Assertions.assertEquals(1, dist.probability(-2, 6)); +<a class="jxr_linenumber" name="L266" href="#L266">266</a> } +<a class="jxr_linenumber" name="L267" href="#L267">267</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/UserGuideTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/UserGuideTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/UserGuideTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,156 @@ +<!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>UserGuideTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/UserGuideTest.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.IntStream; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.apache.commons.rng.UniformRandomProvider; +<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> +<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> <em class="jxr_javadoccomment"> * Test code used in the distributions section of the user guide.</em> +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <strong class="jxr_keyword">class</strong> <a name="UserGuideTest" href="../../../../../org/apache/commons/statistics/distribution/UserGuideTest.html#UserGuideTest">UserGuideTest</a> { +<a class="jxr_linenumber" name="L30" href="#L30">30</a> @Test +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <strong class="jxr_keyword">void</strong> testCDF() { +<a class="jxr_linenumber" name="L32" href="#L32">32</a> TDistribution t = TDistribution.of(29); +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <strong class="jxr_keyword">double</strong> lowerTail = t.cumulativeProbability(-2.656); <em class="jxr_comment">// P(T(29) &lt;= -2.656)</em> +<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">double</strong> upperTail = t.survivalProbability(2.75); <em class="jxr_comment">// P(T(29) &gt; 2.75)</em> +<a class="jxr_linenumber" name="L35" href="#L35">35</a> +<a class="jxr_linenumber" name="L36" href="#L36">36</a> Assertions.assertTrue(lowerTail > upperTail, +<a class="jxr_linenumber" name="L37" href="#L37">37</a> () -> String.format(<span class="jxr_string">"Since 2.75 > |-2.656|, expected %s > %s"</span>, lowerTail, upperTail)); +<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> @Test +<a class="jxr_linenumber" name="L41" href="#L41">41</a> <strong class="jxr_keyword">void</strong> testProbability() { +<a class="jxr_linenumber" name="L42" href="#L42">42</a> PoissonDistribution pd = PoissonDistribution.of(1.23); +<a class="jxr_linenumber" name="L43" href="#L43">43</a> <strong class="jxr_keyword">double</strong> p1 = pd.probability(5); +<a class="jxr_linenumber" name="L44" href="#L44">44</a> <strong class="jxr_keyword">double</strong> p2 = pd.probability(5, 5); +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">double</strong> p3 = pd.probability(4, 5); +<a class="jxr_linenumber" name="L46" href="#L46">46</a> +<a class="jxr_linenumber" name="L47" href="#L47">47</a> Assertions.assertEquals(0, p2); +<a class="jxr_linenumber" name="L48" href="#L48">48</a> Assertions.assertEquals(p1, p3); +<a class="jxr_linenumber" name="L49" href="#L49">49</a> } +<a class="jxr_linenumber" name="L50" href="#L50">50</a> +<a class="jxr_linenumber" name="L51" href="#L51">51</a> @Test +<a class="jxr_linenumber" name="L52" href="#L52">52</a> <strong class="jxr_keyword">void</strong> testInverseCDF() { +<a class="jxr_linenumber" name="L53" href="#L53">53</a> NormalDistribution n = NormalDistribution.of(0, 1); +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <strong class="jxr_keyword">double</strong> x1 = n.inverseCumulativeProbability(1e-300); +<a class="jxr_linenumber" name="L55" href="#L55">55</a> <strong class="jxr_keyword">double</strong> x2 = n.inverseSurvivalProbability(1e-300); +<a class="jxr_linenumber" name="L56" href="#L56">56</a> +<a class="jxr_linenumber" name="L57" href="#L57">57</a> Assertions.assertEquals(x1, -x2); +<a class="jxr_linenumber" name="L58" href="#L58">58</a> Assertions.assertEquals(-37.0471, x1, 1e-3); +<a class="jxr_linenumber" name="L59" href="#L59">59</a> } +<a class="jxr_linenumber" name="L60" href="#L60">60</a> +<a class="jxr_linenumber" name="L61" href="#L61">61</a> @Test +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <strong class="jxr_keyword">void</strong> testProperties() { +<a class="jxr_linenumber" name="L63" href="#L63">63</a> ChiSquaredDistribution chi2 = ChiSquaredDistribution.of(42); +<a class="jxr_linenumber" name="L64" href="#L64">64</a> <strong class="jxr_keyword">double</strong> df = chi2.getDegreesOfFreedom(); <em class="jxr_comment">// 42</em> +<a class="jxr_linenumber" name="L65" href="#L65">65</a> <strong class="jxr_keyword">double</strong> mean = chi2.getMean(); <em class="jxr_comment">// 42</em> +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <strong class="jxr_keyword">double</strong> <strong class="jxr_keyword">var</strong> = chi2.getVariance(); <em class="jxr_comment">// 84</em> +<a class="jxr_linenumber" name="L67" href="#L67">67</a> +<a class="jxr_linenumber" name="L68" href="#L68">68</a> CauchyDistribution cauchy = CauchyDistribution.of(1.23, 4.56); +<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">double</strong> location = cauchy.getLocation(); <em class="jxr_comment">// 1.23</em> +<a class="jxr_linenumber" name="L70" href="#L70">70</a> <strong class="jxr_keyword">double</strong> scale = cauchy.getScale(); <em class="jxr_comment">// 4.56</em> +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">double</strong> undefined1 = cauchy.getMean(); <em class="jxr_comment">// NaN</em> +<a class="jxr_linenumber" name="L72" href="#L72">72</a> <strong class="jxr_keyword">double</strong> undefined2 = cauchy.getVariance(); <em class="jxr_comment">// NaN</em> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> +<a class="jxr_linenumber" name="L74" href="#L74">74</a> Assertions.assertEquals(42, df); +<a class="jxr_linenumber" name="L75" href="#L75">75</a> Assertions.assertEquals(42, mean); +<a class="jxr_linenumber" name="L76" href="#L76">76</a> Assertions.assertEquals(84, <strong class="jxr_keyword">var</strong>); +<a class="jxr_linenumber" name="L77" href="#L77">77</a> Assertions.assertEquals(1.23, location); +<a class="jxr_linenumber" name="L78" href="#L78">78</a> Assertions.assertEquals(4.56, scale); +<a class="jxr_linenumber" name="L79" href="#L79">79</a> Assertions.assertEquals(Double.NaN, undefined1); +<a class="jxr_linenumber" name="L80" href="#L80">80</a> Assertions.assertEquals(Double.NaN, undefined2); +<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> testDomain() { +<a class="jxr_linenumber" name="L85" href="#L85">85</a> BinomialDistribution b = BinomialDistribution.of(13, 0.15); +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">int</strong> lower = b.getSupportLowerBound(); <em class="jxr_comment">// 0</em> +<a class="jxr_linenumber" name="L87" href="#L87">87</a> <strong class="jxr_keyword">int</strong> upper = b.getSupportUpperBound(); <em class="jxr_comment">// 13</em> +<a class="jxr_linenumber" name="L88" href="#L88">88</a> +<a class="jxr_linenumber" name="L89" href="#L89">89</a> Assertions.assertEquals(0, lower); +<a class="jxr_linenumber" name="L90" href="#L90">90</a> Assertions.assertEquals(13, upper); +<a class="jxr_linenumber" name="L91" href="#L91">91</a> } +<a class="jxr_linenumber" name="L92" href="#L92">92</a> +<a class="jxr_linenumber" name="L93" href="#L93">93</a> @Test +<a class="jxr_linenumber" name="L94" href="#L94">94</a> <strong class="jxr_keyword">void</strong> testSampling() { +<a class="jxr_linenumber" name="L95" href="#L95">95</a> <em class="jxr_comment">// From Commons RNG Simple</em> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> UniformRandomProvider rng = RandomSource.KISS.create(123L); +<a class="jxr_linenumber" name="L97" href="#L97">97</a> +<a class="jxr_linenumber" name="L98" href="#L98">98</a> NormalDistribution n = NormalDistribution.of(0, 1); +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <strong class="jxr_keyword">double</strong> x = n.createSampler(rng).sample(); +<a class="jxr_linenumber" name="L100" href="#L100">100</a> +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <em class="jxr_comment">// Generate a number of samples</em> +<a class="jxr_linenumber" name="L102" href="#L102">102</a> GeometricDistribution g = GeometricDistribution.of(0.75); +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <strong class="jxr_keyword">int</strong>[] k = IntStream.generate(g.createSampler(rng)::sample).limit(100).toArray(); +<a class="jxr_linenumber" name="L104" href="#L104">104</a> +<a class="jxr_linenumber" name="L105" href="#L105">105</a> Assertions.assertTrue(-5 < x && x < 5, () -> Double.toString(x)); +<a class="jxr_linenumber" name="L106" href="#L106">106</a> Assertions.assertEquals(100, k.length); +<a class="jxr_linenumber" name="L107" href="#L107">107</a> } +<a class="jxr_linenumber" name="L108" href="#L108">108</a> +<a class="jxr_linenumber" name="L109" href="#L109">109</a> @Test +<a class="jxr_linenumber" name="L110" href="#L110">110</a> <strong class="jxr_keyword">void</strong> testComplement() { +<a class="jxr_linenumber" name="L111" href="#L111">111</a> ChiSquaredDistribution chi2 = ChiSquaredDistribution.of(42); +<a class="jxr_linenumber" name="L112" href="#L112">112</a> <strong class="jxr_keyword">double</strong> q1 = 1 - chi2.cumulativeProbability(168); +<a class="jxr_linenumber" name="L113" href="#L113">113</a> <strong class="jxr_keyword">double</strong> q2 = chi2.survivalProbability(168); +<a class="jxr_linenumber" name="L114" href="#L114">114</a> +<a class="jxr_linenumber" name="L115" href="#L115">115</a> Assertions.assertEquals(0, q1); +<a class="jxr_linenumber" name="L116" href="#L116">116</a> Assertions.assertNotEquals(0, q2); +<a class="jxr_linenumber" name="L117" href="#L117">117</a> +<a class="jxr_linenumber" name="L118" href="#L118">118</a> <em class="jxr_comment">// For the table</em> +<a class="jxr_linenumber" name="L119" href="#L119">119</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> eps = Math.pow(2, -53); +<a class="jxr_linenumber" name="L120" href="#L120">120</a> Assertions.assertEquals(1.110223e-16, eps, 1e-3); +<a class="jxr_linenumber" name="L121" href="#L121">121</a> Assertions.assertEquals(eps, 1 - chi2.cumulativeProbability(166)); +<a class="jxr_linenumber" name="L122" href="#L122">122</a> Assertions.assertEquals(eps, 1 - chi2.cumulativeProbability(167)); +<a class="jxr_linenumber" name="L123" href="#L123">123</a> Assertions.assertEquals(0, 1 - chi2.cumulativeProbability(168)); +<a class="jxr_linenumber" name="L124" href="#L124">124</a> Assertions.assertEquals(0, 1 - chi2.cumulativeProbability(200)); +<a class="jxr_linenumber" name="L125" href="#L125">125</a> Assertions.assertEquals(1.16583e-16, chi2.survivalProbability(166), 1e-3); +<a class="jxr_linenumber" name="L126" href="#L126">126</a> Assertions.assertEquals(7.95907e-17, chi2.survivalProbability(167), 1e-3); +<a class="jxr_linenumber" name="L127" href="#L127">127</a> Assertions.assertEquals(5.42987e-17, chi2.survivalProbability(168), 1e-3); +<a class="jxr_linenumber" name="L128" href="#L128">128</a> Assertions.assertEquals(1.19056e-22, chi2.survivalProbability(200), 1e-3); +<a class="jxr_linenumber" name="L129" href="#L129">129</a> } +<a class="jxr_linenumber" name="L130" href="#L130">130</a> +<a class="jxr_linenumber" name="L131" href="#L131">131</a> @Test +<a class="jxr_linenumber" name="L132" href="#L132">132</a> <strong class="jxr_keyword">void</strong> testInverseComplement() { +<a class="jxr_linenumber" name="L133" href="#L133">133</a> ChiSquaredDistribution chi2 = ChiSquaredDistribution.of(42); +<a class="jxr_linenumber" name="L134" href="#L134">134</a> <strong class="jxr_keyword">double</strong> q = 5.43e-17; +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_comment">// Incorrect: p = 1 - q == 1.0 !!!</em> +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">double</strong> x1 = chi2.inverseCumulativeProbability(1 - q); +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_comment">// Correct: invert q</em> +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <strong class="jxr_keyword">double</strong> x2 = chi2.inverseSurvivalProbability(q); +<a class="jxr_linenumber" name="L139" href="#L139">139</a> +<a class="jxr_linenumber" name="L140" href="#L140">140</a> Assertions.assertEquals(Double.POSITIVE_INFINITY, x1); +<a class="jxr_linenumber" name="L141" href="#L141">141</a> Assertions.assertEquals(168.0, x2, 0.1); +<a class="jxr_linenumber" name="L142" href="#L142">142</a> } +<a class="jxr_linenumber" name="L143" href="#L143">143</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/WeibullDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/WeibullDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/WeibullDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,120 @@ +<!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>WeibullDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/WeibullDistributionTest.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.gamma.LogGamma; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments; +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.CsvSource; +<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 WeibullDistribution}.</em> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. 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="WeibullDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/WeibullDistributionTest.html#WeibullDistributionTest">WeibullDistributionTest</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="L34" href="#L34">34</a> @Override +<a class="jxr_linenumber" name="L35" href="#L35">35</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> shape = (Double) parameters[0]; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> scale = (Double) parameters[1]; +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">return</strong> WeibullDistribution.of(shape, scale); +<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> 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.0, 2.0}, +<a class="jxr_linenumber" name="L46" href="#L46">46</a> {-0.1, 2.0}, +<a class="jxr_linenumber" name="L47" href="#L47">47</a> {1.0, 0.0}, +<a class="jxr_linenumber" name="L48" href="#L48">48</a> {1.0, -0.1}, +<a class="jxr_linenumber" name="L49" href="#L49">49</a> }; +<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> @Override +<a class="jxr_linenumber" name="L53" href="#L53">53</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"Shape"</span>, <span class="jxr_string">"Scale"</span>}; +<a class="jxr_linenumber" name="L55" href="#L55">55</a> } +<a class="jxr_linenumber" name="L56" href="#L56">56</a> +<a class="jxr_linenumber" name="L57" href="#L57">57</a> @Override +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L59" href="#L59">59</a> <strong class="jxr_keyword">return</strong> 1e-14; +<a class="jxr_linenumber" name="L60" href="#L60">60</a> } +<a class="jxr_linenumber" name="L61" href="#L61">61</a> +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L63" href="#L63">63</a> +<a class="jxr_linenumber" name="L64" href="#L64">64</a> @ParameterizedTest +<a class="jxr_linenumber" name="L65" href="#L65">65</a> @MethodSource +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <strong class="jxr_keyword">void</strong> testAdditionalMoments(<strong class="jxr_keyword">double</strong> shape, <strong class="jxr_keyword">double</strong> scale, <strong class="jxr_keyword">double</strong> mean, <strong class="jxr_keyword">double</strong> variance) { +<a class="jxr_linenumber" name="L67" href="#L67">67</a> <strong class="jxr_keyword">final</strong> WeibullDistribution dist = WeibullDistribution.of(shape, scale); +<a class="jxr_linenumber" name="L68" href="#L68">68</a> testMoments(dist, mean, variance, createRelTolerance(1e-15)); +<a class="jxr_linenumber" name="L69" href="#L69">69</a> } +<a class="jxr_linenumber" name="L70" href="#L70">70</a> +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionalMoments() { +<a class="jxr_linenumber" name="L72" href="#L72">72</a> <em class="jxr_comment">// In R: 3.5*gamma(1+(1/2.5)) (or empirically: mean(rweibull(10000, 2.5, 3.5)))</em> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <strong class="jxr_keyword">double</strong> mu1 = 3.5 * Math.exp(LogGamma.value(1 + (1 / 2.5))); +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <strong class="jxr_keyword">double</strong> mu2 = 2.222 * Math.exp(LogGamma.value(1 + (1 / 10.4))); +<a class="jxr_linenumber" name="L75" href="#L75">75</a> <strong class="jxr_keyword">return</strong> Stream.of( +<a class="jxr_linenumber" name="L76" href="#L76">76</a> Arguments.of(2.5, 3.5, mu1, +<a class="jxr_linenumber" name="L77" href="#L77">77</a> (3.5 * 3.5) * +<a class="jxr_linenumber" name="L78" href="#L78">78</a> Math.exp(LogGamma.value(1 + (2 / 2.5))) - +<a class="jxr_linenumber" name="L79" href="#L79">79</a> (mu1 * mu1)), +<a class="jxr_linenumber" name="L80" href="#L80">80</a> Arguments.of(10.4, 2.222, mu2, +<a class="jxr_linenumber" name="L81" href="#L81">81</a> (2.222 * 2.222) * +<a class="jxr_linenumber" name="L82" href="#L82">82</a> Math.exp(LogGamma.value(1 + (2 / 10.4))) - +<a class="jxr_linenumber" name="L83" href="#L83">83</a> (mu2 * mu2)) +<a class="jxr_linenumber" name="L84" href="#L84">84</a> ); +<a class="jxr_linenumber" name="L85" href="#L85">85</a> } +<a class="jxr_linenumber" name="L86" href="#L86">86</a> +<a class="jxr_linenumber" name="L87" href="#L87">87</a> @ParameterizedTest +<a class="jxr_linenumber" name="L88" href="#L88">88</a> @CsvSource({ +<a class="jxr_linenumber" name="L89" href="#L89">89</a> <span class="jxr_string">"1, 2"</span>, +<a class="jxr_linenumber" name="L90" href="#L90">90</a> <span class="jxr_string">"0.1, 2.34"</span>, +<a class="jxr_linenumber" name="L91" href="#L91">91</a> }) +<a class="jxr_linenumber" name="L92" href="#L92">92</a> <strong class="jxr_keyword">void</strong> testAdditionalParameterAccessors(<strong class="jxr_keyword">double</strong> shape, <strong class="jxr_keyword">double</strong> scale) { +<a class="jxr_linenumber" name="L93" href="#L93">93</a> <strong class="jxr_keyword">final</strong> WeibullDistribution dist = WeibullDistribution.of(shape, scale); +<a class="jxr_linenumber" name="L94" href="#L94">94</a> Assertions.assertEquals(shape, dist.getShape()); +<a class="jxr_linenumber" name="L95" href="#L95">95</a> Assertions.assertEquals(scale, dist.getScale()); +<a class="jxr_linenumber" name="L96" href="#L96">96</a> } +<a class="jxr_linenumber" name="L97" href="#L97">97</a> +<a class="jxr_linenumber" name="L98" href="#L98">98</a> @Test +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <strong class="jxr_keyword">void</strong> testInverseCumulativeProbabilitySmallPAccuracy() { +<a class="jxr_linenumber" name="L100" href="#L100">100</a> <strong class="jxr_keyword">final</strong> WeibullDistribution dist = WeibullDistribution.of(2, 3); +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> t = dist.inverseCumulativeProbability(1e-17); +<a class="jxr_linenumber" name="L102" href="#L102">102</a> <em class="jxr_comment">// Analytically, answer is solution to 1e-17 = 1-exp(-(x/3)^2)</em> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <em class="jxr_comment">// x = sqrt(-9*log(1-1e-17))</em> +<a class="jxr_linenumber" name="L104" href="#L104">104</a> <em class="jxr_comment">// If we're not careful, answer will be 0. Answer below is computed with care in Octave:</em> +<a class="jxr_linenumber" name="L105" href="#L105">105</a> Assertions.assertEquals(9.48683298050514e-9, t, 1e-17); +<a class="jxr_linenumber" name="L106" href="#L106">106</a> } +<a class="jxr_linenumber" name="L107" href="#L107">107</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/ZipfDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ZipfDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ZipfDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,161 @@ +<!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>ZipfDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/ZipfDistributionTest.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.math3.util.MathArrays; +<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 ZipfDistribution}.</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="ZipfDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/ZipfDistributionTest.html#ZipfDistributionTest">ZipfDistributionTest</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> n = (Integer) parameters[0]; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> e = (Double) parameters[1]; +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">return</strong> ZipfDistribution.of(n, e); +<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> 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, 1.0}, +<a class="jxr_linenumber" name="L46" href="#L46">46</a> {-1, 1.0}, +<a class="jxr_linenumber" name="L47" href="#L47">47</a> {1, -0.1}, +<a class="jxr_linenumber" name="L48" href="#L48">48</a> }; +<a class="jxr_linenumber" name="L49" href="#L49">49</a> } +<a class="jxr_linenumber" name="L50" href="#L50">50</a> +<a class="jxr_linenumber" name="L51" href="#L51">51</a> @Override +<a class="jxr_linenumber" name="L52" href="#L52">52</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L53" href="#L53">53</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"NumberOfElements"</span>, <span class="jxr_string">"Exponent"</span>}; +<a class="jxr_linenumber" name="L54" href="#L54">54</a> } +<a class="jxr_linenumber" name="L55" href="#L55">55</a> +<a class="jxr_linenumber" name="L56" href="#L56">56</a> @Override +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <strong class="jxr_keyword">return</strong> 1e-14; +<a class="jxr_linenumber" name="L59" href="#L59">59</a> } +<a class="jxr_linenumber" name="L60" href="#L60">60</a> +<a class="jxr_linenumber" name="L61" href="#L61">61</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L62" href="#L62">62</a> +<a class="jxr_linenumber" name="L63" href="#L63">63</a> @ParameterizedTest +<a class="jxr_linenumber" name="L64" href="#L64">64</a> @MethodSource +<a class="jxr_linenumber" name="L65" href="#L65">65</a> <strong class="jxr_keyword">void</strong> testAdditionlSurvivalProbabilityHighPrecision(<strong class="jxr_keyword">int</strong> n, <strong class="jxr_keyword">double</strong> e, <strong class="jxr_keyword">int</strong>[] x, <strong class="jxr_keyword">double</strong>[] expected, <a name="DoubleTolerance" href="../../../../../org/apache/commons/statistics/distribution/DoubleTolerance.html#DoubleTolerance">DoubleTolerance</a> tol) { +<a class="jxr_linenumber" name="L66" href="#L66">66</a> testSurvivalProbabilityHighPrecision(ZipfDistribution.of(n, e), x, expected, tol); +<a class="jxr_linenumber" name="L67" href="#L67">67</a> } +<a class="jxr_linenumber" name="L68" href="#L68">68</a> +<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionlSurvivalProbabilityHighPrecision() { +<a class="jxr_linenumber" name="L70" href="#L70">70</a> <em class="jxr_comment">// computed using scipy.stats (1.7.1) zipfian</em> +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">return</strong> Stream.of( +<a class="jxr_linenumber" name="L72" href="#L72">72</a> Arguments.of(60, 10, +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {57, 59}, +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {2.3189337454689757e-18, 1.6521739576668957e-18}, +<a class="jxr_linenumber" name="L75" href="#L75">75</a> DoubleTolerances.absolute(1e-25)), +<a class="jxr_linenumber" name="L76" href="#L76">76</a> Arguments.of(60, 50.5, +<a class="jxr_linenumber" name="L77" href="#L77">77</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {57, 59}, +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {8.8488396450491320e-90, 1.5972093932264611e-90}, +<a class="jxr_linenumber" name="L79" href="#L79">79</a> DoubleTolerances.absolute(1e-95)), +<a class="jxr_linenumber" name="L80" href="#L80">80</a> Arguments.of(60, 100.5, +<a class="jxr_linenumber" name="L81" href="#L81">81</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {57, 59}, +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {5.9632998443758656e-178, 1.9760564023408183e-179}, +<a class="jxr_linenumber" name="L83" href="#L83">83</a> DoubleTolerances.absolute(1e-185)) +<a class="jxr_linenumber" name="L84" href="#L84">84</a> ); +<a class="jxr_linenumber" name="L85" href="#L85">85</a> } +<a class="jxr_linenumber" name="L86" href="#L86">86</a> +<a class="jxr_linenumber" name="L87" href="#L87">87</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L88" href="#L88">88</a> <em class="jxr_javadoccomment"> * Test the high precision survival probability computation when the exponent creates</em> +<a class="jxr_linenumber" name="L89" href="#L89">89</a> <em class="jxr_javadoccomment"> * an overflow in the intermediate. The result should not be infinite or NaN and</em> +<a class="jxr_linenumber" name="L90" href="#L90">90</a> <em class="jxr_javadoccomment"> * it should be a complement to the CDF value.</em> +<a class="jxr_linenumber" name="L91" href="#L91">91</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L92" href="#L92">92</a> @Test +<a class="jxr_linenumber" name="L93" href="#L93">93</a> <strong class="jxr_keyword">void</strong> testAdditionalSurvivalAndCumulativeProbabilityComplement() { +<a class="jxr_linenumber" name="L94" href="#L94">94</a> <em class="jxr_comment">// Requires (x+1)^a to overflow</em> +<a class="jxr_linenumber" name="L95" href="#L95">95</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> n = 60; +<a class="jxr_linenumber" name="L96" href="#L96">96</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = 200.5; +<a class="jxr_linenumber" name="L97" href="#L97">97</a> Assertions.assertEquals(Double.POSITIVE_INFINITY, Math.pow(n, a)); +<a class="jxr_linenumber" name="L98" href="#L98">98</a> ZipfDistribution dist = ZipfDistribution.of(n, a); +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] points = MathArrays.natural(n); +<a class="jxr_linenumber" name="L100" href="#L100">100</a> testSurvivalAndCumulativeProbabilityComplement(dist, points, createTolerance()); +<a class="jxr_linenumber" name="L101" href="#L101">101</a> } +<a class="jxr_linenumber" name="L102" href="#L102">102</a> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L104" href="#L104">104</a> <em class="jxr_javadoccomment"> * Test sampling for various number of points and exponents.</em> +<a class="jxr_linenumber" name="L105" href="#L105">105</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L106" href="#L106">106</a> @Test +<a class="jxr_linenumber" name="L107" href="#L107">107</a> <strong class="jxr_keyword">void</strong> testSamplingExtended() { +<a class="jxr_linenumber" name="L108" href="#L108">108</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> sampleSize = 1000; +<a class="jxr_linenumber" name="L109" href="#L109">109</a> +<a class="jxr_linenumber" name="L110" href="#L110">110</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] numPointsValues = { +<a class="jxr_linenumber" name="L111" href="#L111">111</a> 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40, 45, 50, 60, 70, 80, 90, 100 +<a class="jxr_linenumber" name="L112" href="#L112">112</a> }; +<a class="jxr_linenumber" name="L113" href="#L113">113</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] exponentValues = { +<a class="jxr_linenumber" name="L114" href="#L114">114</a> 1e-10, 1e-9, 1e-8, 1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1, 2e-1, 5e-1, +<a class="jxr_linenumber" name="L115" href="#L115">115</a> 1. - 1e-9, 1.0, 1. + 1e-9, 1.1, 1.2, 1.3, 1.5, 1.6, 1.7, 1.8, 2.0, +<a class="jxr_linenumber" name="L116" href="#L116">116</a> 2.5, 3.0, 4., 5., 6., 7., 8., 9., 10., 20., 30., 100., 150. +<a class="jxr_linenumber" name="L117" href="#L117">117</a> }; +<a class="jxr_linenumber" name="L118" href="#L118">118</a> +<a class="jxr_linenumber" name="L119" href="#L119">119</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> numPoints : numPointsValues) { +<a class="jxr_linenumber" name="L120" href="#L120">120</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> exponent : exponentValues) { +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <strong class="jxr_keyword">double</strong> weightSum = 0.; +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] weights = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[numPoints]; +<a class="jxr_linenumber" name="L123" href="#L123">123</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = numPoints; i >= 1; i -= 1) { +<a class="jxr_linenumber" name="L124" href="#L124">124</a> weights[i - 1] = Math.pow(i, -exponent); +<a class="jxr_linenumber" name="L125" href="#L125">125</a> weightSum += weights[i - 1]; +<a class="jxr_linenumber" name="L126" href="#L126">126</a> } +<a class="jxr_linenumber" name="L127" href="#L127">127</a> +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <em class="jxr_comment">// Use fixed seed, the test is expected to fail for more than 50% of all</em> +<a class="jxr_linenumber" name="L129" href="#L129">129</a> <em class="jxr_comment">// seeds because each test case can fail with probability 0.001, the chance</em> +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_comment">// that all test cases do not fail is 0.999^(32*22) = 0.49442874426</em> +<a class="jxr_linenumber" name="L131" href="#L131">131</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler distribution = +<a class="jxr_linenumber" name="L132" href="#L132">132</a> ZipfDistribution.of(numPoints, exponent).createSampler( +<a class="jxr_linenumber" name="L133" href="#L133">133</a> RandomSource.XO_SHI_RO_256_PP.create(1)); +<a class="jxr_linenumber" name="L134" href="#L134">134</a> +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] expectedCounts = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[numPoints]; +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong>[] observedCounts = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">long</strong>[numPoints]; +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < numPoints; i++) { +<a class="jxr_linenumber" name="L138" href="#L138">138</a> expectedCounts[i] = sampleSize * (weights[i] / weightSum); +<a class="jxr_linenumber" name="L139" href="#L139">139</a> } +<a class="jxr_linenumber" name="L140" href="#L140">140</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] sample = TestUtils.sample(sampleSize, distribution); +<a class="jxr_linenumber" name="L141" href="#L141">141</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> s : sample) { +<a class="jxr_linenumber" name="L142" href="#L142">142</a> observedCounts[s - 1]++; +<a class="jxr_linenumber" name="L143" href="#L143">143</a> } +<a class="jxr_linenumber" name="L144" href="#L144">144</a> TestUtils.assertChiSquareAccept(expectedCounts, observedCounts, 0.001); +<a class="jxr_linenumber" name="L145" href="#L145">145</a> } +<a class="jxr_linenumber" name="L146" href="#L146">146</a> } +<a class="jxr_linenumber" name="L147" href="#L147">147</a> } +<a class="jxr_linenumber" name="L148" href="#L148">148</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>
