17 type(display_type) :: disp
18 integer(IK) :: nsam1, nsam2
19 integer(IK) :: itry, ntry
= 10
23 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
24 call disp%show(
"! Compute the KS probability of sample originating from a Uniform distribution in range `[0, 1)`.")
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
30 integer(IK),
allocatable :: iweight1(:)
31 real(TKG),
allocatable :: rweight1(:)
32 real(TKG),
allocatable :: sample1(:)
33 real(TKG) :: probKS, quanKS, statKS
35 call disp%show(
"nsam1 = getUnifRand(1, 10)")
37 call disp%show(
"sample1 = getUnifRand(0., 1., nsam1)")
41 call disp%show(
"statKS = getDisKolm(sample1) ! assuming unweighted samples.")
45 call disp%show(
"call setProbKS(probKS, quanKS, statKS, nsam1) ! assuming unweighted samples.")
46 call setProbKS(probKS, quanKS, statKS, nsam1)
50 call disp%show(
"iweight1 = getUnifRand(1, 9, nsam1)")
54 call disp%show(
"rweight1 = iweight1")
56 call disp%show(
"statKS = getDisKolm(sample1, iweight1, sum(iweight1))")
57 statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
60 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(iweight1))")
61 call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
64 call disp%show(
"statKS = getDisKolm(sample1, rweight1, sum(rweight1))")
65 statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
68 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(rweight1), sum(rweight1**2))")
69 call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
77 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
78 call disp%show(
"! Compute the KS probability of a sample against a Normal distribution.")
79 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
84 integer(IK),
allocatable :: iweight1(:)
85 real(TKG),
allocatable :: rweight1(:)
86 real(TKG),
allocatable :: sample1(:)
87 real(TKG) :: probKS, quanKS, statKS
89 call disp%show(
"nsam1 = getUnifRand(5, 10)")
91 call disp%show(
"sample1 = getNormRand(mean = getFilled(0., nsam1))")
95 call disp%show(
"statKS = getDisKolm(sample1, getNormCDF_RKS) ! assuming unweighted samples.")
99 call disp%show(
"call setProbKS(probKS, quanKS, statKS, nsam1) ! assuming unweighted samples.")
100 call setProbKS(probKS, quanKS, statKS, nsam1)
104 call disp%show(
"iweight1 = getUnifRand(1, 9, nsam1)")
108 call disp%show(
"rweight1 = iweight1")
110 call disp%show(
"statKS = getDisKolm(sample1, iweight1, sum(iweight1), getNormCDF_RKS)")
111 statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
114 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(iweight1))")
115 call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
118 call disp%show(
"statKS = getDisKolm(sample1, rweight1, sum(rweight1), getNormCDF_RKS)")
119 statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
122 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(rweight1), sum(rweight1**2))")
123 call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
131 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
132 call disp%show(
"! Compute the two sample KS probability.")
133 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
138 integer(IK),
allocatable :: iweight1(:), iweight2(:)
139 real(TKG),
allocatable :: rweight1(:), rweight2(:)
140 real(TKG),
allocatable :: sample1(:), sample2(:)
141 real(TKG) :: probKS, quanKS, statKS
143 call disp%show(
"nsam1 = getUnifRand(1, 10); nsam2 = getUnifRand(1, 10)")
145 call disp%show(
"sample1 = getUnifRand(0., 1., nsam1)")
149 call disp%show(
"sample2 = getUnifRand(0., 1., nsam2)")
153 call disp%show(
"statKS = getDisKolm(sample1, sample2) ! assuming unweighted samples.")
157 call disp%show(
"call setProbKS(probKS, quanKS, statKS, nsam1, nsam2) ! assuming unweighted samples.")
158 call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
162 call disp%show(
"iweight1 = getUnifRand(1, 9, nsam1)")
166 call disp%show(
"rweight1 = iweight1")
168 call disp%show(
"statKS = getDisKolm(sample1, iweight1, sum(iweight1), sample2)")
169 statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
172 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(iweight1), nsam2)")
173 call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
176 call disp%show(
"statKS = getDisKolm(sample1, rweight1, sum(rweight1), sample2)")
177 statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
180 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(rweight1), nsam2, sum(rweight1**2))")
181 call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
185 call disp%show(
"iweight2 = getUnifRand(1, 9, nsam2)")
189 call disp%show(
"rweight2 = iweight2")
191 call disp%show(
"statKS = getDisKolm(sample1, iweight1, sum(iweight1), sample2, iweight2, sum(iweight2))")
192 statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
195 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(iweight1), sum(iweight2))")
196 call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
199 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(rweight1), sum(iweight2), sum(rweight1**2))")
200 call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
203 call disp%show(
"statKS = getDisKolm(sample1, rweight1, sum(rweight1), sample2, rweight2, sum(rweight2))")
204 statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
207 call disp%show(
"call setProbKS(probKS, quanKS, statKS, sum(rweight1), sum(rweight2), sum(rweight1**2), sum(rweight2**2))")
208 call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
217 function getNormCDF_RKS(x)
result(cdf)
220 real(RKG),
intent(in) :: x
225 function getUnifCDF_RKS(x)
result(cdf)
228 real(RKG),
intent(in) :: x
Generate and return an array of the specified rank and shape of arbitrary intrinsic type and kind wit...
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate and return the sorted elements of the input scalar string or contiguous vector in ascending ...
Generate an equally-weighted (verbose or flattened) array of the input weighted array of rank 1 or 2.
Generate and return the Cumulative Distribution Function (CDF) of the univariate Normal distribution.
Generate and return a scalar or array of arbitrary rank of random values from the univariate Normal d...
Generate and return the Cumulative Distribution Function (CDF) of a univariate Standard Uniform distr...
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Generate and return the Kolmogorov distance of a sample1 of size nsam1 from another sample sample2 of...
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
This module contains procedures and generic interfaces for convenient allocation and filling of array...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for various sorting tasks.
This module contains procedures and generic interfaces for flattening (duplicating the elements of) a...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for computing the Kolmogorov statistical distance.
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
Generate and return an object of type display_type.
6nsam1
= getUnifRand(
1,
10)
7sample1
= getUnifRand(
0.,
1., nsam1)
13call setProbKS(probKS, quanKS, statKS, nsam1)
15+0.361705840,
+0.923046768
17iweight1
= getUnifRand(
1,
9, nsam1)
21statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
24call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
26+0.970500708E-2,
+1.63221681
27statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
30call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
32+0.361705840,
+0.923046768
34nsam1
= getUnifRand(
1,
10)
35sample1
= getUnifRand(
0.,
1., nsam1)
37+0.841557503,
+0.644738078,
+0.872535110E-1,
+0.593524575,
+0.801190555,
+0.363895416,
+0.654039562
41call setProbKS(probKS, quanKS, statKS, nsam1)
43+0.444119096,
+0.864124179
45iweight1
= getUnifRand(
1,
9, nsam1)
47+6,
+5,
+3,
+7,
+8,
+9,
+8
49statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
52call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
54+0.501275063E-4,
+2.30148983
55statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
58call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
60+0.393775344,
+0.899244010
62nsam1
= getUnifRand(
1,
10)
63sample1
= getUnifRand(
0.,
1., nsam1)
65+0.158734679,
+0.778290927
69call setProbKS(probKS, quanKS, statKS, nsam1)
71+0.922698200,
+0.550118089
73iweight1
= getUnifRand(
1,
9, nsam1)
77statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
80call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
82+0.234961510E-3,
+2.12711716
83statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
86call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
88+0.375344157,
+0.912768483
90nsam1
= getUnifRand(
1,
10)
91sample1
= getUnifRand(
0.,
1., nsam1)
93+0.490797877,
+0.610904098E-1,
+0.178990126,
+0.192978442
97call setProbKS(probKS, quanKS, statKS, nsam1)
99+0.106184959,
+1.21152186
101iweight1
= getUnifRand(
1,
9, nsam1)
105statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
108call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
110+0.113993287E-1,
+1.60738134
111statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
114call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
116+0.185917497,
+1.08967745
118nsam1
= getUnifRand(
1,
10)
119sample1
= getUnifRand(
0.,
1., nsam1)
121+0.760805011E-1,
+0.738783360,
+0.575072885,
+0.429375708,
+0.180011272,
+0.869740188
125call setProbKS(probKS, quanKS, statKS, nsam1)
127+0.997105777,
+0.400844783
129iweight1
= getUnifRand(
1,
9, nsam1)
131+7,
+1,
+4,
+4,
+4,
+2
133statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
136call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
138+0.167069435E-1,
+1.54678333
139statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
142call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
144+0.624479890,
+0.751596630
146nsam1
= getUnifRand(
1,
10)
147sample1
= getUnifRand(
0.,
1., nsam1)
149+0.555587590,
+0.185561538,
+0.268361509,
+0.799393177,
+0.512172997,
+0.694734454
153call setProbKS(probKS, quanKS, statKS, nsam1)
155+0.946110487,
+0.524465859
157iweight1
= getUnifRand(
1,
9, nsam1)
159+7,
+9,
+4,
+9,
+8,
+9
161statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
164call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
166+0.128818154E-1,
+1.58825183
167statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
170call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
172+0.882701755,
+0.585643947
174nsam1
= getUnifRand(
1,
10)
175sample1
= getUnifRand(
0.,
1., nsam1)
177+0.213127136E-1,
+0.811957777,
+0.474444687,
+0.213783145,
+0.462083697
181call setProbKS(probKS, quanKS, statKS, nsam1)
183+0.571957409,
+0.783045590
185iweight1
= getUnifRand(
1,
9, nsam1)
189statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
192call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
194+0.248998404E-2,
+1.82874501
195statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
198call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
200+0.675803661,
+0.721147537
202nsam1
= getUnifRand(
1,
10)
203sample1
= getUnifRand(
0.,
1., nsam1)
205+0.250843346,
+0.946807146,
+0.533075809,
+0.442408502,
+0.467391491,
+0.493385077,
+0.194902897,
+0.662707448,
+0.126336217E-1,
+0.849234760
209call setProbKS(probKS, quanKS, statKS, nsam1)
211+0.919054627,
+0.553697944
213iweight1
= getUnifRand(
1,
9, nsam1)
215+8,
+4,
+4,
+7,
+8,
+2,
+3,
+9,
+2,
+3
217statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
220call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
222+0.350856781E-1,
+1.42181289
223statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
226call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
228+0.881862402,
+0.586319327
230nsam1
= getUnifRand(
1,
10)
231sample1
= getUnifRand(
0.,
1., nsam1)
233+0.374048710,
+0.764225721,
+0.837047338,
+0.413477421E-2,
+0.868876219,
+0.771119058,
+0.726211667E-1
237call setProbKS(probKS, quanKS, statKS, nsam1)
239+0.337042451,
+0.942291379
241iweight1
= getUnifRand(
1,
9, nsam1)
243+9,
+7,
+8,
+5,
+7,
+3,
+7
245statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
248call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
250+0.231504440E-3,
+2.12886477
251statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
254call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
256+0.489709258,
+0.834148705
258nsam1
= getUnifRand(
1,
10)
259sample1
= getUnifRand(
0.,
1., nsam1)
261+0.447513461,
+0.247067630,
+0.908288717,
+0.374045551,
+0.271606803,
+0.462080121,
+0.146921694,
+0.882143974E-1,
+0.740250707
265call setProbKS(probKS, quanKS, statKS, nsam1)
267+0.273713827,
+0.996552289
269iweight1
= getUnifRand(
1,
9, nsam1)
271+8,
+5,
+7,
+6,
+6,
+5,
+3,
+8,
+1
273statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
276call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
278+0.125169754E-5,
+2.67342758
279statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
282call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
284+0.174543738,
+1.10409713
291nsam1
= getUnifRand(
5,
10)
292sample1
= getNormRand(mean
= getFilled(
0., nsam1))
294+0.136388883,
+0.605567694,
+0.404625863,
+0.756903112,
+0.121487409,
-1.22627771
298call setProbKS(probKS, quanKS, statKS, nsam1)
300+0.272295237,
+0.997865140
302iweight1
= getUnifRand(
1,
9, nsam1)
304+5,
+4,
+4,
+6,
+1,
+7
306statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
309call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
311+0.171067119E-1,
+1.54295683
312statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
315call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
317+0.708747506,
+0.701493502
319nsam1
= getUnifRand(
5,
10)
320sample1
= getNormRand(mean
= getFilled(
0., nsam1))
322-1.04057145,
+1.82433641,
+2.53895116,
-0.816127479,
+0.863076866,
-1.39803541,
+0.696354628,
+1.75025964,
-1.15852964,
+1.16431034
326call setProbKS(probKS, quanKS, statKS, nsam1)
328+0.121229529,
+1.18384838
330iweight1
= getUnifRand(
1,
9, nsam1)
332+4,
+8,
+7,
+7,
+1,
+4,
+7,
+4,
+1,
+7
334statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
337call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
339+0.00000000,
+3.14854980
340statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
343call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
345+0.646111965E-1,
+1.31005371
347nsam1
= getUnifRand(
5,
10)
348sample1
= getNormRand(mean
= getFilled(
0., nsam1))
350+1.51843762,
+0.729435623,
+1.86606550,
-0.757727101E-1,
+0.320981413,
-0.609201491,
+0.147198111,
-0.837240040,
+0.360612899
354call setProbKS(probKS, quanKS, statKS, nsam1)
356+0.574480355,
+0.781520367
358iweight1
= getUnifRand(
1,
9, nsam1)
360+6,
+9,
+9,
+9,
+2,
+1,
+5,
+9,
+4
362statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
365call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
367+0.229477882E-3,
+2.12989712
368statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
371call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
373+0.530221641,
+0.808596313
375nsam1
= getUnifRand(
5,
10)
376sample1
= getNormRand(mean
= getFilled(
0., nsam1))
378+0.692703247,
-0.669361472,
-1.39046586,
+0.490267962,
-1.26310551,
-0.811202288,
-0.160265043,
+1.52044880
382call setProbKS(probKS, quanKS, statKS, nsam1)
384+0.640726447,
+0.741952777
386iweight1
= getUnifRand(
1,
9, nsam1)
388+1,
+7,
+8,
+6,
+2,
+8,
+3,
+2
390statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
393call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
395+0.178813934E-5,
+2.63790798
396statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
399call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
401+0.175895274,
+1.10234559
403nsam1
= getUnifRand(
5,
10)
404sample1
= getNormRand(mean
= getFilled(
0., nsam1))
406+0.780188739,
-0.114481777,
-1.20702732,
+0.602363721E-1,
-0.651725233,
+0.177392483,
-1.56050837,
+0.787280440
410call setProbKS(probKS, quanKS, statKS, nsam1)
412+0.801329732,
+0.643942893
414iweight1
= getUnifRand(
1,
9, nsam1)
416+3,
+9,
+9,
+8,
+3,
+1,
+3,
+6
418statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
421call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
423+0.140721798E-1,
+1.57427776
424statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
427call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
429+0.828294694,
+0.625866771
431nsam1
= getUnifRand(
5,
10)
432sample1
= getNormRand(mean
= getFilled(
0., nsam1))
434+1.65811586,
-0.600132048,
+1.07393289,
+1.59756124,
-0.420666814
438call setProbKS(probKS, quanKS, statKS, nsam1)
440+0.175397635,
+1.10298932
442iweight1
= getUnifRand(
1,
9, nsam1)
446statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
449call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
451+0.454783440E-4,
+2.31206369
452statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
455call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
457+0.361312151,
+0.923347116
459nsam1
= getUnifRand(
5,
10)
460sample1
= getNormRand(mean
= getFilled(
0., nsam1))
462-0.586329043,
+0.984112263,
+1.10989738,
+0.124827862,
+2.29397345,
+0.238590166,
+0.576355994,
+0.894027948,
-0.572956860,
-0.729726970
466call setProbKS(probKS, quanKS, statKS, nsam1)
468+0.499060452,
+0.828171313
470iweight1
= getUnifRand(
1,
9, nsam1)
472+2,
+6,
+1,
+6,
+4,
+9,
+8,
+9,
+2,
+3
474statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
477call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
479+0.596046448E-7,
+2.95233846
480statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
483call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
485+0.117876351,
+1.18976080
487nsam1
= getUnifRand(
5,
10)
488sample1
= getNormRand(mean
= getFilled(
0., nsam1))
490-0.442256600,
+0.839802086,
+0.884793103,
+1.87960017,
-0.398809940
494call setProbKS(probKS, quanKS, statKS, nsam1)
496+0.314313829,
+0.960878611
498iweight1
= getUnifRand(
1,
9, nsam1)
502statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
505call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
507+0.00000000,
+3.19033813
508statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
511call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
513+0.563831329E-1,
+1.33579755
515nsam1
= getUnifRand(
5,
10)
516sample1
= getNormRand(mean
= getFilled(
0., nsam1))
518+0.490279287,
+0.226869941,
+0.119143456,
+0.767090559,
+2.50107622,
-1.58514166,
+0.948038340,
-1.78756726,
-1.73981082
522call setProbKS(probKS, quanKS, statKS, nsam1)
524+0.429648578,
+0.873974442
526iweight1
= getUnifRand(
1,
9, nsam1)
528+5,
+8,
+8,
+9,
+8,
+9,
+4,
+7,
+3
530statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
533call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
535+0.544786453E-3,
+2.02587104
536statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
539call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
541+0.591543794,
+0.771249831
543nsam1
= getUnifRand(
5,
10)
544sample1
= getNormRand(mean
= getFilled(
0., nsam1))
546+1.36534190,
-2.05721831,
+0.706652641,
+2.85573173,
+0.250651240,
-0.123919308,
-0.197943762,
+0.653735101,
+0.397777975
550call setProbKS(probKS, quanKS, statKS, nsam1)
552+0.292132497,
+0.979934871
554iweight1
= getUnifRand(
1,
9, nsam1)
556+9,
+9,
+9,
+2,
+2,
+9,
+3,
+4,
+6
558statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
561call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
563+0.188153982E-2,
+1.86665726
564statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
567call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
569+0.688537002,
+0.713573337
576nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
577sample1
= getUnifRand(
0.,
1., nsam1)
579+0.725660682,
+0.386335015,
+0.353589177,
+0.889835060,
+0.550283134
580sample2
= getUnifRand(
0.,
1., nsam2)
586call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
588+0.362166107,
+0.922695935
590iweight1
= getUnifRand(
1,
9, nsam1)
594statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
597call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
599+0.669028938,
+0.725169897
600statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
603call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
605+0.744714260,
+0.679696739
607iweight2
= getUnifRand(
1,
9, nsam2)
611statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
614call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
616+0.669028938,
+0.725169897
617call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
619+0.744714260,
+0.679696739
620statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
623call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
625+0.744714260,
+0.679696739
627nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
628sample1
= getUnifRand(
0.,
1., nsam1)
630+0.994117320,
+0.510379612,
+0.981798708,
+0.766516864,
+0.243390024,
+0.728488863
631sample2
= getUnifRand(
0.,
1., nsam2)
633+0.436963081,
+0.424052179,
+0.859232068,
+0.966595292,
+0.780610621,
+0.529457927,
+0.986880779,
+0.849796593,
+0.582064927,
+0.174764156
637call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
639+0.968186021,
+0.493102282
641iweight1
= getUnifRand(
1,
9, nsam1)
643+9,
+4,
+8,
+7,
+4,
+9
645statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
648call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
650+0.337330997,
+0.942061067
651statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
654call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
656+0.795125067,
+0.647986054
658iweight2
= getUnifRand(
1,
9, nsam2)
660+4,
+9,
+6,
+6,
+4,
+3,
+9,
+5,
+7,
+3
662statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
665call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
667+0.528317690E-1,
+1.34791911
668call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
670+0.797113299,
+0.646694601
671statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
674call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
676+0.930143535,
+0.542505622
678nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
679sample1
= getUnifRand(
0.,
1., nsam1)
681+0.185329139,
+0.661870956,
+0.229286671
682sample2
= getUnifRand(
0.,
1., nsam2)
684+0.334512770,
+0.547339618,
+0.729911864,
+0.550948262,
+0.979734600,
+0.762786150
688call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
690+0.198367178,
+1.07466364
692iweight1
= getUnifRand(
1,
9, nsam1)
696statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
699call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
701+0.161435723,
+1.12166584
702statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
705call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
707+0.643276572,
+0.740440488
709iweight2
= getUnifRand(
1,
9, nsam2)
711+8,
+8,
+1,
+2,
+8,
+3
713statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
716call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
718+0.591423512E-1,
+1.32682514
719call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
721+0.784895837,
+0.654573083
722statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
725call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
727+0.901241183,
+0.570089817
729nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
730sample1
= getUnifRand(
0.,
1., nsam1)
732+0.220915675E-1,
+0.256340861,
+0.824155390
733sample2
= getUnifRand(
0.,
1., nsam2)
735+0.707827032,
+0.735501528,
+0.594927311
739call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
741+0.319724321,
+0.956373096
743iweight1
= getUnifRand(
1,
9, nsam1)
747statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
750call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
752+0.711643100E-1,
+1.29148483
753statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
756call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
758+0.247217238,
+1.02194273
760iweight2
= getUnifRand(
1,
9, nsam2)
764statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
767call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
769+0.745058060E-5,
+2.49964929
770call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
772+0.704075098E-1,
+1.29355311
773statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
776call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
778+0.253495038,
+1.01575279
780nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
781sample1
= getUnifRand(
0.,
1., nsam1)
783+0.471772552,
+0.533313811,
+0.979318619,
+0.232793093
784sample2
= getUnifRand(
0.,
1., nsam2)
786+0.412247002,
+0.340840816,
+0.661848009,
+0.631533146,
+0.472228229,
+0.638137341,
+0.822668970
790call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
792+0.897244513,
+0.573555171
794iweight1
= getUnifRand(
1,
9, nsam1)
798statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
801call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
803+0.706809402,
+0.702655911
804statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
807call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
809+0.968598187,
+0.492401004
811iweight2
= getUnifRand(
1,
9, nsam2)
813+4,
+1,
+3,
+9,
+6,
+9,
+9
815statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
818call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
820+0.507402420E-2,
+1.72869456
821call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
823+0.414376736,
+0.884575248
824statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
827call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
829+0.642189682,
+0.741084993
831nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
832sample1
= getUnifRand(
0.,
1., nsam1)
834+0.104163170,
+0.402716398E-1,
+0.769958138,
+0.656907082,
+0.405296147,
+0.925457478E-2,
+0.585104287
835sample2
= getUnifRand(
0.,
1., nsam2)
837+0.956998885,
+0.592527926,
+0.185544491,
+0.867063165,
+0.385210156,
+0.383065343,
+0.414294243
841call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
843+0.423218310,
+0.878411233
845iweight1
= getUnifRand(
1,
9, nsam1)
847+6,
+4,
+9,
+5,
+9,
+6,
+3
849statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
852call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
854+0.274354696,
+0.995960891
855statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
858call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
860+0.610514164,
+0.759908319
862iweight2
= getUnifRand(
1,
9, nsam2)
864+3,
+8,
+8,
+8,
+9,
+4,
+8
866statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
869call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
871+0.201338530E-2,
+1.85756409
872call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
874+0.319010973,
+0.956964016
875statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
878call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
880+0.635967255,
+0.744775832
882nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
883sample1
= getUnifRand(
0.,
1., nsam1)
885+0.745877802,
+0.229861557,
+0.281458616,
+0.608282864,
+0.171462417,
+0.545058608,
+0.473988771,
+0.401905775,
+0.952116907
886sample2
= getUnifRand(
0.,
1., nsam2)
888+0.763356626,
+0.928640425,
+0.496933281,
+0.160847783,
+0.144419253,
+0.689561665,
+0.962090194,
+0.854294658
892call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
894+0.394428611,
+0.898771822
896iweight1
= getUnifRand(
1,
9, nsam1)
898+9,
+4,
+2,
+6,
+7,
+6,
+8,
+8,
+6
900statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
903call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
905+0.175483108,
+1.10287869
906statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
909call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
911+0.454138935,
+0.857405424
913iweight2
= getUnifRand(
1,
9, nsam2)
915+4,
+4,
+7,
+3,
+8,
+9,
+6,
+1
917statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
920call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
922+0.184753537E-1,
+1.53043497
923call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
925+0.478657722,
+0.841282308
926statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
929call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
931+0.821491599,
+0.630513668
933nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
934sample1
= getUnifRand(
0.,
1., nsam1)
936+0.652330875,
+0.508425772,
+0.105229914,
+0.137712717,
+0.406633317,
+0.340529203,
+0.374366641
937sample2
= getUnifRand(
0.,
1., nsam2)
939+0.671436191E-1,
+0.237999558,
+0.600627184,
+0.991357327,
+0.594674885,
+0.767674208
943call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
945+0.232980132,
+1.03642654
947iweight1
= getUnifRand(
1,
9, nsam1)
949+1,
+6,
+5,
+6,
+8,
+6,
+2
951statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
954call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
956+0.167223215E-1,
+1.54663455
957statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
960call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
962+0.107979298,
+1.20805740
964iweight2
= getUnifRand(
1,
9, nsam2)
966+1,
+1,
+7,
+9,
+8,
+4
968statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
971call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
973+0.00000000,
+3.74199462
974call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
976+0.218808651E-3,
+2.13545394
977statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
980call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
982+0.132895112E-1,
+1.58333957
984nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
985sample1
= getUnifRand(
0.,
1., nsam1)
987+0.176995516,
+0.545700431,
+0.310479403
988sample2
= getUnifRand(
0.,
1., nsam2)
990+0.327256083,
+0.532386601,
+0.938651979,
+0.545160174E-1
994call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
996+0.821437240,
+0.630550563
998iweight1
= getUnifRand(
1,
9, nsam1)
1002statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1005call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
1007+0.555401444E-1,
+1.33861399
1008statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1011call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1013+0.329634130,
+0.948253751
1015iweight2
= getUnifRand(
1,
9, nsam2)
1019statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1022call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
1024+0.770688057E-4,
+2.25432634
1025call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1027+0.979957581E-1,
+1.22797811
1028statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
1031call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
1033+0.219936013,
+1.05029273
1035nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
1036sample1
= getUnifRand(
0.,
1., nsam1)
1039sample2
= getUnifRand(
0.,
1., nsam2)
1041+0.577523351,
+0.281509161,
+0.718317807,
+0.916631460,
+0.534032226,
+0.780738473,
+0.138916552,
+0.857305825
1045call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
1047+0.414325297,
+0.884611249
1049iweight1
= getUnifRand(
1,
9, nsam1)
1053statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1056call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
1058+0.174691439,
+1.10390520
1059statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1062call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1064+0.414325297,
+0.884611249
1066iweight2
= getUnifRand(
1,
9, nsam2)
1068+4,
+6,
+2,
+6,
+5,
+3,
+3,
+9
1070statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1073call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
1075+0.321106136,
+0.955230653
1076call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1078+0.648179233,
+0.737533748
1079statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
1082call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
1084+0.696993232,
+0.708529413