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) :: statKS, probKS
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(
"probKS = getProbKS(statKS, nsam1) ! assuming unweighted samples.")
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(
"probKS = getProbKS(statKS, sum(iweight1))")
64 call disp%show(
"statKS = getDisKolm(sample1, rweight1, sum(rweight1))")
65 statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
68 call disp%show(
"probKS = getProbKS(statKS, sum(rweight1), sum(rweight1**2))")
69 probKS
= getProbKS(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) :: statKS, probKS
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(
"probKS = getProbKS(statKS, nsam1) ! assuming unweighted samples.")
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(
"probKS = getProbKS(statKS, sum(iweight1))")
115 probKS
= getProbKS(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(
"probKS = getProbKS(statKS, sum(rweight1), sum(rweight1**2))")
123 probKS
= getProbKS(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) :: statKS, probKS
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(
"probKS = getProbKS(statKS, nsam1, nsam2) ! assuming unweighted samples.")
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(
"probKS = getProbKS(statKS, sum(iweight1), nsam2)")
173 probKS
= getProbKS(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(
"probKS = getProbKS(statKS, sum(rweight1), nsam2, sum(rweight1**2))")
181 probKS
= getProbKS(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(
"probKS = getProbKS(statKS, sum(iweight1), sum(iweight2))")
196 probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
199 call disp%show(
"probKS = getProbKS(statKS, sum(rweight1), sum(iweight2), sum(rweight1**2))")
200 probKS
= getProbKS(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(
"probKS = getProbKS(statKS, sum(rweight1), sum(rweight2), sum(rweight1**2), sum(rweight2**2))")
208 probKS
= getProbKS(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)
9+0.760522485E-1,
+0.242145419,
+0.911414087,
+0.802966714,
+0.518629551E-1,
+0.532067001,
+0.607682347
17iweight1
= getUnifRand(
1,
9, nsam1)
19+2,
+4,
+9,
+7,
+2,
+7,
+5
21statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
27statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
30probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
34nsam1
= getUnifRand(
1,
10)
35sample1
= getUnifRand(
0.,
1., nsam1)
37+0.432508171,
+0.261734426,
+0.201474965,
+0.868292809
45iweight1
= getUnifRand(
1,
9, nsam1)
49statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
55statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
58probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
62nsam1
= getUnifRand(
1,
10)
63sample1
= getUnifRand(
0.,
1., nsam1)
65+0.759366393,
+0.925849974,
+0.627367437,
+0.982025325
73iweight1
= getUnifRand(
1,
9, nsam1)
77statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
83statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
86probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
90nsam1
= getUnifRand(
1,
10)
91sample1
= getUnifRand(
0.,
1., nsam1)
93+0.865392148,
+0.805033326,
+0.542861819,
+0.289520323,
+0.188381970,
+0.181042314
101iweight1
= getUnifRand(
1,
9, nsam1)
103+9,
+4,
+6,
+9,
+6,
+4
105statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
111statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
114probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
118nsam1
= getUnifRand(
1,
10)
119sample1
= getUnifRand(
0.,
1., nsam1)
121+0.611972570,
+0.661492288,
+0.283098519,
+0.871342480,
+0.266463399,
+0.518711269,
+0.691151857,
+0.214590907
129iweight1
= getUnifRand(
1,
9, nsam1)
131+1,
+4,
+7,
+7,
+3,
+1,
+6,
+9
133statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
139statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
142probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
146nsam1
= getUnifRand(
1,
10)
147sample1
= getUnifRand(
0.,
1., nsam1)
149+0.412885129,
+0.666892529,
+0.890827835
157iweight1
= getUnifRand(
1,
9, nsam1)
161statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
167statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
170probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
174nsam1
= getUnifRand(
1,
10)
175sample1
= getUnifRand(
0.,
1., nsam1)
177+0.748530746,
+0.847966194,
+0.182071030,
+0.921644092,
+0.550287366
185iweight1
= getUnifRand(
1,
9, nsam1)
189statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
195statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
198probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
202nsam1
= getUnifRand(
1,
10)
203sample1
= getUnifRand(
0.,
1., nsam1)
205+0.973555446,
+0.270028114,
+0.440104485,
+0.126195967,
+0.965398848,
+0.572558343,
+0.995203555,
+0.682451844
213iweight1
= getUnifRand(
1,
9, nsam1)
215+1,
+9,
+3,
+3,
+5,
+7,
+2,
+9
217statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
223statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
226probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
230nsam1
= getUnifRand(
1,
10)
231sample1
= getUnifRand(
0.,
1., nsam1)
233+0.660184443,
+0.403744042,
+0.866553605,
+0.448995292,
+0.596399307,
+0.745425522,
+0.661064982,
+0.442012906,
+0.948382437,
+0.588714361
241iweight1
= getUnifRand(
1,
9, nsam1)
243+9,
+9,
+6,
+4,
+2,
+7,
+2,
+5,
+9,
+6
245statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
251statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
254probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
258nsam1
= getUnifRand(
1,
10)
259sample1
= getUnifRand(
0.,
1., nsam1)
261+0.602278173,
+0.325556159,
+0.434645295
269iweight1
= getUnifRand(
1,
9, nsam1)
273statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
279statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
282probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
291nsam1
= getUnifRand(
5,
10)
292sample1
= getNormRand(mean
= getFilled(
0., nsam1))
294+0.680385292,
+1.06130874,
-1.05992270,
+0.206461400,
+1.09774375
302iweight1
= getUnifRand(
1,
9, nsam1)
306statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
312statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
315probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
319nsam1
= getUnifRand(
5,
10)
320sample1
= getNormRand(mean
= getFilled(
0., nsam1))
322-1.17966509,
+1.33958066,
+0.285816044,
-2.86359000,
+0.469453603,
+1.11978781,
-0.927450478
330iweight1
= getUnifRand(
1,
9, nsam1)
332+4,
+6,
+8,
+5,
+6,
+9,
+2
334statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
340statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
343probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
347nsam1
= getUnifRand(
5,
10)
348sample1
= getNormRand(mean
= getFilled(
0., nsam1))
350+1.58451307,
-0.373223066,
-0.614647388,
-0.849718004E-1,
+0.542600632
358iweight1
= getUnifRand(
1,
9, nsam1)
362statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
368statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
371probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
375nsam1
= getUnifRand(
5,
10)
376sample1
= getNormRand(mean
= getFilled(
0., nsam1))
378+0.128500499E-1,
-0.556562662,
-0.123393029,
-0.125062624E-1,
-1.57089305,
+0.873300016
386iweight1
= getUnifRand(
1,
9, nsam1)
388+8,
+5,
+3,
+1,
+1,
+7
390statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
396statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
399probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
403nsam1
= getUnifRand(
5,
10)
404sample1
= getNormRand(mean
= getFilled(
0., nsam1))
406-1.86525989,
+0.162649155,
+1.14998841,
+1.49993563,
+0.438667208,
-0.336715341,
-0.559750676,
-0.562392473,
-0.808587253,
-1.53268886
414iweight1
= getUnifRand(
1,
9, nsam1)
416+4,
+7,
+7,
+4,
+3,
+7,
+5,
+6,
+3,
+5
418statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
424statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
427probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
431nsam1
= getUnifRand(
5,
10)
432sample1
= getNormRand(mean
= getFilled(
0., nsam1))
434+0.634779692,
-0.337696642,
+0.101046368,
+0.282361150,
-0.940403759,
+0.265650421,
+0.499229645E-2
442iweight1
= getUnifRand(
1,
9, nsam1)
444+1,
+6,
+3,
+7,
+7,
+5,
+1
446statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
452statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
455probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
459nsam1
= getUnifRand(
5,
10)
460sample1
= getNormRand(mean
= getFilled(
0., nsam1))
462-1.01495504,
-0.971003354,
+0.576281488,
+0.493916184,
+2.62294531,
+0.699197650,
+0.397261202,
-1.45296216
470iweight1
= getUnifRand(
1,
9, nsam1)
472+8,
+5,
+8,
+6,
+2,
+6,
+8,
+2
474statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
480statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
483probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
487nsam1
= getUnifRand(
5,
10)
488sample1
= getNormRand(mean
= getFilled(
0., nsam1))
490-0.521996617,
-0.218877383E-1,
+0.680643320,
-0.454270273,
+1.23284388,
+1.27680838,
-2.36960268,
+1.23671365,
+1.67219639
498iweight1
= getUnifRand(
1,
9, nsam1)
500+3,
+3,
+2,
+1,
+6,
+8,
+2,
+2,
+8
502statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
508statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
511probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
515nsam1
= getUnifRand(
5,
10)
516sample1
= getNormRand(mean
= getFilled(
0., nsam1))
518+0.281620681,
+0.343992114,
-0.586772561,
-1.00303471,
+0.218603581,
-1.69188166
526iweight1
= getUnifRand(
1,
9, nsam1)
528+5,
+9,
+3,
+1,
+1,
+2
530statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
536statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
539probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
543nsam1
= getUnifRand(
5,
10)
544sample1
= getNormRand(mean
= getFilled(
0., nsam1))
546+0.420961350,
+0.185761973,
+0.450357907E-1,
+0.319852233,
-1.46804130,
-0.434015602,
+0.535362959
554iweight1
= getUnifRand(
1,
9, nsam1)
556+6,
+5,
+9,
+7,
+3,
+4,
+4
558statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
564statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
567probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight1
**2))
576nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
577sample1
= getUnifRand(
0.,
1., nsam1)
579+0.918609083,
+0.175046325,
+0.473087370,
+0.242241025E-1,
+0.404218793,
+0.510570526,
+0.750517607,
+0.333082676E-2
580sample2
= getUnifRand(
0.,
1., nsam2)
582+0.140347481E-1,
+0.928643167,
+0.249650657,
+0.414809823
590iweight1
= getUnifRand(
1,
9, nsam1)
592+6,
+2,
+6,
+1,
+1,
+7,
+1,
+5
594statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
597probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
600statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
603probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
607iweight2
= getUnifRand(
1,
9, nsam2)
611statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
614probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
617probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
620statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
623probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
627nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
628sample1
= getUnifRand(
0.,
1., nsam1)
630+0.964345694,
+0.455367208,
+0.454771280
631sample2
= getUnifRand(
0.,
1., nsam2)
641iweight1
= getUnifRand(
1,
9, nsam1)
645statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
648probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
651statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
654probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
658iweight2
= getUnifRand(
1,
9, nsam2)
662statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
665probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
668probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
671statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
674probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
678nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
679sample1
= getUnifRand(
0.,
1., nsam1)
681+0.677905142,
+0.123868823,
+0.968774855,
+0.474617839,
+0.152571082,
+0.811890304,
+0.416277945
682sample2
= getUnifRand(
0.,
1., nsam2)
684+0.440616727,
+0.150049627,
+0.728440046
692iweight1
= getUnifRand(
1,
9, nsam1)
694+2,
+3,
+1,
+8,
+4,
+5,
+4
696statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
699probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
702statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
705probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
709iweight2
= getUnifRand(
1,
9, nsam2)
713statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
716probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
719probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
722statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
725probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
729nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
730sample1
= getUnifRand(
0.,
1., nsam1)
732+0.348249674E-1,
+0.678867698,
+0.971817613,
+0.770888627,
+0.552349925,
+0.387537420,
+0.953372598,
+0.138764381E-1,
+0.356970549
733sample2
= getUnifRand(
0.,
1., nsam2)
735+0.866714239,
+0.792740583E-1,
+0.926134109,
+0.905420780E-1,
+0.383863330
743iweight1
= getUnifRand(
1,
9, nsam1)
745+4,
+6,
+1,
+8,
+9,
+9,
+3,
+8,
+7
747statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
750probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
753statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
756probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
760iweight2
= getUnifRand(
1,
9, nsam2)
764statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
767probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
770probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
773statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
776probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
780nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
781sample1
= getUnifRand(
0.,
1., nsam1)
783+0.909213722,
+0.135219395,
+0.868196726,
+0.901540220,
+0.256987274,
+0.681144714,
+0.737713516,
+0.835244656
784sample2
= getUnifRand(
0.,
1., nsam2)
786+0.889343023E-2,
+0.185074627,
+0.998635352,
+0.985233128,
+0.933409154
794iweight1
= getUnifRand(
1,
9, nsam1)
796+9,
+3,
+4,
+8,
+6,
+7,
+5,
+8
798statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
801probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
804statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
807probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
811iweight2
= getUnifRand(
1,
9, nsam2)
815statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
818probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
821probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
824statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
827probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
831nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
832sample1
= getUnifRand(
0.,
1., nsam1)
834+0.395591557,
+0.530764639,
+0.202414036
835sample2
= getUnifRand(
0.,
1., nsam2)
837+0.211610913,
+0.128943264,
+0.532164216,
+0.403368592,
+0.951641977,
+0.401474893,
+0.589747608,
+0.336771011,
+0.399106860
845iweight1
= getUnifRand(
1,
9, nsam1)
849statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
852probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
855statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
858probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
862iweight2
= getUnifRand(
1,
9, nsam2)
864+8,
+9,
+3,
+4,
+1,
+6,
+2,
+5,
+4
866statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
869probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
872probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
875statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
878probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
882nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
883sample1
= getUnifRand(
0.,
1., nsam1)
885+0.180443525,
+0.364023328,
+0.170207798
886sample2
= getUnifRand(
0.,
1., nsam2)
888+0.788386047,
+0.967619956,
+0.626059949,
+0.426645577
896iweight1
= getUnifRand(
1,
9, nsam1)
900statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
903probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
906statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
909probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
913iweight2
= getUnifRand(
1,
9, nsam2)
917statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
920probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
923probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
926statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
929probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
933nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
934sample1
= getUnifRand(
0.,
1., nsam1)
936+0.236142457,
+0.397683144,
+0.580960214,
+0.701507449,
+0.480057001,
+0.555683672,
+0.111082137,
+0.977445781,
+0.707425237
937sample2
= getUnifRand(
0.,
1., nsam2)
939+0.348372102,
+0.577573121,
+0.531984627
947iweight1
= getUnifRand(
1,
9, nsam1)
949+4,
+2,
+5,
+5,
+1,
+5,
+2,
+7,
+9
951statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
954probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
957statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
960probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
964iweight2
= getUnifRand(
1,
9, nsam2)
968statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
971probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
974probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
977statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
980probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
984nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
985sample1
= getUnifRand(
0.,
1., nsam1)
987+0.838498235,
+0.392711699,
+0.562042117
988sample2
= getUnifRand(
0.,
1., nsam2)
990+0.841958284,
+0.854712307
998iweight1
= getUnifRand(
1,
9, nsam1)
1002statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1005probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
1008statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1011probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1015iweight2
= getUnifRand(
1,
9, nsam2)
1019statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1022probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
1025probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1028statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
1031probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))
1035nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
1036sample1
= getUnifRand(
0.,
1., nsam1)
1038+0.844324172,
+0.180904210,
+0.193701684,
+0.103821874,
+0.237490118,
+0.282863975E-1,
+0.613570452,
+0.142697513,
+0.266947806
1039sample2
= getUnifRand(
0.,
1., nsam2)
1041+0.893250287,
+0.429074585,
+0.950456500,
+0.804113150
1049iweight1
= getUnifRand(
1,
9, nsam1)
1051+6,
+5,
+5,
+6,
+1,
+3,
+3,
+9,
+5
1053statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1056probKS
= getProbKS(statKS,
sum(iweight1), nsam2)
1059statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1062probKS
= getProbKS(statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1066iweight2
= getUnifRand(
1,
9, nsam2)
1070statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1073probKS
= getProbKS(statKS,
sum(iweight1),
sum(iweight2))
1076probKS
= getProbKS(statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1079statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2, rweight2,
sum(rweight2))
1082probKS
= getProbKS(statKS,
sum(rweight1),
sum(rweight2),
sum(rweight1
**2),
sum(rweight2
**2))