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)
17iweight1
= getUnifRand(
1,
9, nsam1)
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.282225609E-1,
+0.248199701E-1,
+0.694696665,
+0.298018217,
+0.203974426
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.390291989,
+0.272128820,
+0.178313613,
+0.866233766
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.418912232,
+0.779013634E-1,
+0.603253841E-1
101iweight1
= getUnifRand(
1,
9, nsam1)
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.280966341,
+0.363789022,
+0.552643776
129iweight1
= getUnifRand(
1,
9, nsam1)
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.879443884E-1,
+0.506753325
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.748675108,
+0.990334749,
+0.304757237,
+0.194653034,
+0.750824153
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.183262825E-1,
+0.413981676,
+0.889078021,
+0.992694378
213iweight1
= getUnifRand(
1,
9, nsam1)
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)
241iweight1
= getUnifRand(
1,
9, nsam1)
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.927959740,
+0.681383073,
+0.525358558,
+0.667170703,
+0.186861336,
+0.738574564
269iweight1
= getUnifRand(
1,
9, nsam1)
271+6,
+5,
+6,
+6,
+6,
+5
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.611315429,
+0.744792461,
+0.186102405,
+1.00170028,
-0.274076432,
+0.500599146,
+0.192758575,
-0.782971442
302iweight1
= getUnifRand(
1,
9, nsam1)
304+3,
+7,
+4,
+1,
+8,
+3,
+2,
+2
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-0.415161341,
-1.20703471,
-1.36520422,
+0.316389531,
+3.42308068,
+0.348236829,
+2.69583273,
+1.60451663
330iweight1
= getUnifRand(
1,
9, nsam1)
332+9,
+6,
+6,
+2,
+5,
+5,
+8,
+3
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+0.486906022E-1,
-0.799988568,
+0.598107427E-1,
+0.771048188,
+0.106948316,
+0.416614451E-1,
-0.399641693,
-0.871242821,
-1.41021359,
-0.995726764
358iweight1
= getUnifRand(
1,
9, nsam1)
360+9,
+4,
+4,
+5,
+9,
+2,
+1,
+4,
+7,
+4
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+1.09967041,
+0.160599038,
-1.21712554,
+0.534283854E-1,
+0.278614074,
-0.451634914,
-1.94911301,
+0.205989987,
-0.610692620
386iweight1
= getUnifRand(
1,
9, nsam1)
388+8,
+3,
+1,
+2,
+9,
+6,
+1,
+9,
+8
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.55475152,
+0.636219755E-1,
-0.706179857,
-0.442644775,
-1.03599489,
-2.05279970,
-0.669263899,
-0.910537541
414iweight1
= getUnifRand(
1,
9, nsam1)
416+4,
+5,
+6,
+4,
+3,
+1,
+6,
+4
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.787099659,
+0.620190382,
+0.538606420E-1,
-0.121582575,
-0.447140932,
+0.987637699,
-0.425790995E-1,
-1.17335820
442iweight1
= getUnifRand(
1,
9, nsam1)
444+8,
+7,
+6,
+8,
+3,
+7,
+6,
+4
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-0.662077367,
-0.320136309,
-0.696508467,
+0.682200730,
+0.690784156
470iweight1
= getUnifRand(
1,
9, nsam1)
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+2.41306520,
-1.00105071,
+0.310237646,
+0.886133313,
-0.550336421
498iweight1
= getUnifRand(
1,
9, nsam1)
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+1.03488457,
-0.354379743,
-0.628104210,
-0.828207731,
-1.00397944,
-0.248476490,
-0.403357409E-1
526iweight1
= getUnifRand(
1,
9, nsam1)
528+1,
+6,
+6,
+5,
+5,
+2,
+3
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.339250684,
+0.587200522,
+1.44581890,
+0.214912698,
-1.00134695
554iweight1
= getUnifRand(
1,
9, nsam1)
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.934231997,
+0.425366282,
+0.383253515,
+0.640175104,
+0.813117325,
+0.267769456,
+0.791174710,
+0.764962852,
+0.200153947,
+0.757163405
580sample2
= getUnifRand(
0.,
1., nsam2)
582+0.591842711,
+0.341469049E-2,
+0.431611955,
+0.982068956,
+0.338265896,
+0.934400976
590iweight1
= getUnifRand(
1,
9, nsam1)
592+6,
+7,
+5,
+3,
+8,
+8,
+8,
+2,
+8,
+6
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)
609+7,
+7,
+7,
+1,
+8,
+7
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.144581914,
+0.723874748,
+0.982209921,
+0.335025728,
+0.267206848,
+0.562319458,
+0.721729457
631sample2
= getUnifRand(
0.,
1., nsam2)
633+0.208220780,
+0.759055138,
+0.915828645,
+0.621930838,
+0.998717427,
+0.330597699,
+0.994642019,
+0.983479083,
+0.240846574
641iweight1
= getUnifRand(
1,
9, nsam1)
643+3,
+3,
+9,
+4,
+8,
+4,
+7
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)
660+3,
+3,
+1,
+8,
+4,
+7,
+9,
+7,
+1
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.933533072,
+0.647750139,
+0.986048520,
+0.371065497,
+0.474206209E-1,
+0.239380002E-1,
+0.312445700,
+0.318487883
682sample2
= getUnifRand(
0.,
1., nsam2)
684+0.633330107,
+0.207078755,
+0.374108016,
+0.334971786,
+0.437211394,
+0.633070171,
+0.595945120E-1,
+0.292220712,
+0.257508337,
+0.921576738
692iweight1
= getUnifRand(
1,
9, nsam1)
694+7,
+4,
+2,
+1,
+8,
+3,
+5,
+5
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)
711+6,
+9,
+8,
+4,
+2,
+6,
+1,
+4,
+8,
+6
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.668588996,
+0.771800697,
+0.663918674,
+0.570029616E-1,
+0.446385980,
+0.914100409E-1,
+0.124847531,
+0.154986560
733sample2
= getUnifRand(
0.,
1., nsam2)
735+0.811889052,
+0.280447006,
+0.361677587,
+0.287623227,
+0.929665089,
+0.234318316,
+0.874609351
743iweight1
= getUnifRand(
1,
9, nsam1)
745+4,
+7,
+2,
+7,
+3,
+6,
+4,
+3
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)
762+3,
+5,
+2,
+4,
+9,
+7,
+2
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.330772519,
+0.872871459,
+0.582035840,
+0.617837131,
+0.488358974,
+0.866998076,
+0.575729847,
+0.446105480,
+0.115672529
784sample2
= getUnifRand(
0.,
1., nsam2)
794iweight1
= getUnifRand(
1,
9, nsam1)
796+7,
+4,
+3,
+6,
+7,
+4,
+9,
+6,
+3
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.649592698,
+0.793084264,
+0.902936220,
+0.335958004
835sample2
= getUnifRand(
0.,
1., nsam2)
837+0.190795302,
+0.690337300,
+0.326023102,
+0.330879152
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)
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.158787251,
+0.503024340,
+0.185945153
886sample2
= getUnifRand(
0.,
1., nsam2)
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.633870065,
+0.792132318,
+0.296993196,
+0.646297634,
+0.532432199E-1,
+0.398275852E-1,
+0.940089941,
+0.197900534E-1
937sample2
= getUnifRand(
0.,
1., nsam2)
939+0.362821817E-1,
+0.638007402,
+0.554875433,
+0.335331380,
+0.514576554,
+0.897355914
947iweight1
= getUnifRand(
1,
9, nsam1)
949+8,
+6,
+5,
+2,
+9,
+9,
+7,
+2
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)
966+7,
+6,
+2,
+3,
+2,
+6
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.536787868,
+0.958667457,
+0.605838418,
+0.713409603,
+0.439950526,
+0.195576012,
+0.658283830E-1,
+0.564221442,
+0.367155850,
+0.747975886
988sample2
= getUnifRand(
0.,
1., nsam2)
990+0.895007849E-1,
+0.294569075,
+0.515230775,
+0.666793823,
+0.367471397
998iweight1
= getUnifRand(
1,
9, nsam1)
1000+3,
+8,
+6,
+3,
+6,
+7,
+6,
+7,
+2,
+6
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.484630883,
+0.577707827,
+0.989745736,
+0.658459961,
+0.559964180
1039sample2
= getUnifRand(
0.,
1., nsam2)
1041+0.953786075,
+0.656581283,
+0.331942439,
+0.241036177,
+0.756387353,
+0.861041546
1049iweight1
= getUnifRand(
1,
9, nsam1)
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)
1068+9,
+1,
+7,
+2,
+2,
+7
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))