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)
9+0.203030109,
+0.641373158,
+0.205537736,
+0.731497705,
+0.163228095
13call setProbKS(probKS, quanKS, statKS, nsam1)
15+0.328979611,
+0.948784888
17iweight1
= getUnifRand(
1,
9, nsam1)
21statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
24call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
26+0.107296288,
+1.20937049
27statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
30call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
32+0.846547425,
+0.613048792
34nsam1
= getUnifRand(
1,
10)
35sample1
= getUnifRand(
0.,
1., nsam1)
37+0.682456672,
+0.478744507E-1,
+0.288029611,
+0.375700414,
+0.901460648E-3
41call setProbKS(probKS, quanKS, statKS, nsam1)
43+0.248618126,
+1.02055144
45iweight1
= getUnifRand(
1,
9, nsam1)
49statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
52call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
54+0.161528587E-4,
+2.42132521
55statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
58call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
60+0.246435523,
+1.02272165
62nsam1
= getUnifRand(
1,
10)
63sample1
= getUnifRand(
0.,
1., nsam1)
65+0.832525790,
+0.135855794,
+0.946642458,
+0.782246590E-1
69call setProbKS(probKS, quanKS, statKS, nsam1)
71+0.557189465,
+0.792013645
73iweight1
= getUnifRand(
1,
9, nsam1)
77statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
80call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
82+0.147759914E-3,
+2.18094921
83statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
86call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
88+0.412093282,
+0.886179388
90nsam1
= getUnifRand(
1,
10)
91sample1
= getUnifRand(
0.,
1., nsam1)
93+0.944220364,
+0.763934255E-1
97call setProbKS(probKS, quanKS, statKS, nsam1)
99+0.684325099,
+0.716081142
101iweight1
= getUnifRand(
1,
9, nsam1)
105statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
108call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
110+0.615888834E-2,
+1.70044208
111statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
114call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
116+0.399799585,
+0.894908011
118nsam1
= getUnifRand(
1,
10)
119sample1
= getUnifRand(
0.,
1., nsam1)
121+0.405632257,
+0.624966621,
+0.539201975,
+0.344439685,
+0.542382598,
+0.203894913,
+0.449537039,
+0.972304344E-1,
+0.709795713,
+0.159487426
125call setProbKS(probKS, quanKS, statKS, nsam1)
127+0.312233090,
+0.962625802
129iweight1
= getUnifRand(
1,
9, nsam1)
131+5,
+7,
+6,
+7,
+3,
+9,
+7,
+2,
+7,
+2
133statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
136call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
138+0.134944916E-3,
+2.19134140
139statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
142call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
144+0.402615905,
+0.892894626
146nsam1
= getUnifRand(
1,
10)
147sample1
= getUnifRand(
0.,
1., nsam1)
149+0.897032499,
+0.134232879,
+0.987594426,
+0.567355573,
+0.534966528
153call setProbKS(probKS, quanKS, statKS, nsam1)
155+0.534925818,
+0.805682063
157iweight1
= getUnifRand(
1,
9, nsam1)
161statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
164call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
166+0.760138035E-2,
+1.66921711
167statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
170call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
172+0.766618252,
+0.666131735
174nsam1
= getUnifRand(
1,
10)
175sample1
= getUnifRand(
0.,
1., nsam1)
177+0.610787988,
+0.509104133E-1
181call setProbKS(probKS, quanKS, statKS, nsam1)
183+0.671117425,
+0.723930299
185iweight1
= getUnifRand(
1,
9, nsam1)
189statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
192call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
194+0.913235545E-1,
+1.24225569
195statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
198call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
200+0.820538759,
+0.631159425
202nsam1
= getUnifRand(
1,
10)
203sample1
= getUnifRand(
0.,
1., nsam1)
205+0.180020332,
+0.318874121,
+0.306072235,
+0.246218085,
+0.597390354,
+0.673498034,
+0.972510934,
+0.211614966E-1,
+0.642313540,
+0.572485328
209call setProbKS(probKS, quanKS, statKS, nsam1)
211+0.624943018,
+0.751321316
213iweight1
= getUnifRand(
1,
9, nsam1)
215+1,
+2,
+8,
+1,
+1,
+4,
+5,
+6,
+6,
+8
217statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
220call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
222+0.461236238E-1,
+1.37287378
223statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
226call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
228+0.881347835,
+0.586732328
230nsam1
= getUnifRand(
1,
10)
231sample1
= getUnifRand(
0.,
1., nsam1)
233+0.296568871E-2,
+0.927554488
237call setProbKS(probKS, quanKS, statKS, nsam1)
239+0.542163670,
+0.801217020
241iweight1
= getUnifRand(
1,
9, nsam1)
245statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
248call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
250+0.395025611E-1,
+1.40080869
251statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
254call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
256+0.547954321,
+0.797659993
258nsam1
= getUnifRand(
1,
10)
259sample1
= getUnifRand(
0.,
1., nsam1)
261+0.279667258,
+0.642632544,
+0.273894668E-1,
+0.171967864,
+0.447345734,
+0.561631978,
+0.896118164,
+0.216700613
265call setProbKS(probKS, quanKS, statKS, nsam1)
267+0.720941603,
+0.694155455
269iweight1
= getUnifRand(
1,
9, nsam1)
271+8,
+6,
+3,
+3,
+1,
+4,
+4,
+8
273statKS
= getDisKolm(sample1, iweight1,
sum(iweight1))
276call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
278+0.927388668E-3,
+1.95911467
279statKS
= getDisKolm(sample1, rweight1,
sum(rweight1))
282call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
284+0.471109271,
+0.846200764
291nsam1
= getUnifRand(
5,
10)
292sample1
= getNormRand(mean
= getFilled(
0., nsam1))
294-0.392691314,
-0.981695205E-1,
-0.730077684,
+0.165299803,
-0.646332622,
+0.637657270E-1,
+0.153356731,
-0.123063497,
+1.11554623
298call setProbKS(probKS, quanKS, statKS, nsam1)
300+0.248801291,
+1.02037001
302iweight1
= getUnifRand(
1,
9, nsam1)
304+7,
+8,
+6,
+7,
+6,
+6,
+8,
+6,
+6
306statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
309call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
311+0.184774399E-5,
+2.63476610
312statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
315call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
317+0.221964359,
+1.04809642
319nsam1
= getUnifRand(
5,
10)
320sample1
= getNormRand(mean
= getFilled(
0., nsam1))
322+0.975554168,
+1.60844648,
+1.93492508,
-2.42087150,
-0.724315464
326call setProbKS(probKS, quanKS, statKS, nsam1)
328+0.222844481,
+1.04714811
330iweight1
= getUnifRand(
1,
9, nsam1)
334statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
337call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
339+0.111460686E-4,
+2.45920086
340statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
343call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
345+0.192721426,
+1.08137858
347nsam1
= getUnifRand(
5,
10)
348sample1
= getNormRand(mean
= getFilled(
0., nsam1))
350-0.707451165,
-0.400844477E-1,
-0.820323765,
+0.856934190,
+0.168776348,
-0.649737239
354call setProbKS(probKS, quanKS, statKS, nsam1)
356+0.717443228,
+0.696265280
358iweight1
= getUnifRand(
1,
9, nsam1)
360+1,
+5,
+1,
+5,
+3,
+8
362statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
365call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
367+0.206914246,
+1.06477380
368statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
371call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
373+0.974857330,
+0.480917752
375nsam1
= getUnifRand(
5,
10)
376sample1
= getNormRand(mean
= getFilled(
0., nsam1))
378-0.614290178,
-0.239296019,
-1.74087512,
+0.488004059,
-0.525458232E-1
382call setProbKS(probKS, quanKS, statKS, nsam1)
384+0.590333045,
+0.771976113
386iweight1
= getUnifRand(
1,
9, nsam1)
390statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
393call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
395+0.564802885E-1,
+1.33547533
396statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
399call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
401+0.786186516,
+0.653747022
403nsam1
= getUnifRand(
5,
10)
404sample1
= getNormRand(mean
= getFilled(
0., nsam1))
406-0.773658514,
+0.160499141,
+0.289900899,
-0.862072289,
+1.74731541,
+0.530319870,
+0.447142661,
-0.497536182,
+0.557672560,
-2.20405555
410call setProbKS(probKS, quanKS, statKS, nsam1)
412+0.829003096,
+0.625379145
414iweight1
= getUnifRand(
1,
9, nsam1)
416+7,
+3,
+7,
+6,
+8,
+6,
+4,
+8,
+7,
+7
418statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
421call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
423+0.667278171E-1,
+1.30388713
424statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
427call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
429+0.948744595,
+0.521201491
431nsam1
= getUnifRand(
5,
10)
432sample1
= getNormRand(mean
= getFilled(
0., nsam1))
434-1.66792583,
-0.295180291,
+0.609067619,
-0.677218497,
-0.448818028,
+1.96987820,
+0.833300412,
-0.180877954,
-0.520074069
438call setProbKS(probKS, quanKS, statKS, nsam1)
440+0.622691154,
+0.752659798
442iweight1
= getUnifRand(
1,
9, nsam1)
444+1,
+7,
+5,
+2,
+8,
+1,
+6,
+8,
+9
446statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
449call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
451+0.111937523E-3,
+2.21251726
452statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
455call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
457+0.426217198,
+0.876337409
459nsam1
= getUnifRand(
5,
10)
460sample1
= getNormRand(mean
= getFilled(
0., nsam1))
462-0.438886225,
+0.607452929,
+1.35194397,
-0.286789656,
-0.132875949
466call setProbKS(probKS, quanKS, statKS, nsam1)
468+0.552902341,
+0.794630945
470iweight1
= getUnifRand(
1,
9, nsam1)
474statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
477call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
479+0.172734261E-3,
+2.16295910
480statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
483call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
485+0.440730393,
+0.866414964
487nsam1
= getUnifRand(
5,
10)
488sample1
= getNormRand(mean
= getFilled(
0., nsam1))
490+0.125861421,
+0.450504363,
+0.932299674,
+1.17400992,
+0.491232216
494call setProbKS(probKS, quanKS, statKS, nsam1)
496+0.603265166E-1,
+1.32308412
498iweight1
= getUnifRand(
1,
9, nsam1)
502statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
505call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
507+0.774860382E-6,
+2.71671367
508statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
511call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
513+0.138429165,
+1.15546787
515nsam1
= getUnifRand(
5,
10)
516sample1
= getNormRand(mean
= getFilled(
0., nsam1))
518+0.178229481,
-1.04232657,
+0.130884767,
-0.901564479,
+0.587696657E-1
522call setProbKS(probKS, quanKS, statKS, nsam1)
524+0.236767888,
+1.03250992
526iweight1
= getUnifRand(
1,
9, nsam1)
530statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
533call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
535+0.171005726E-3,
+2.16414309
536statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
539call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
541+0.294862688,
+0.977535605
543nsam1
= getUnifRand(
5,
10)
544sample1
= getNormRand(mean
= getFilled(
0., nsam1))
546-0.430932611,
+0.120028920,
-0.611642003,
+0.434246629,
-0.149966434,
-0.594306707
550call setProbKS(probKS, quanKS, statKS, nsam1)
552+0.438213170,
+0.868122756
554iweight1
= getUnifRand(
1,
9, nsam1)
556+4,
+9,
+1,
+6,
+1,
+3
558statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), getNormCDF_RKS)
561call setProbKS(probKS, quanKS, statKS,
sum(iweight1))
563+0.736057758E-2,
+1.67403150
564statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), getNormCDF_RKS)
567call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(rweight1
**2))
569+0.673999786,
+0.722218931
576nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
577sample1
= getUnifRand(
0.,
1., nsam1)
579+0.353881717,
+0.549148619,
+0.814085484,
+0.342118919
580sample2
= getUnifRand(
0.,
1., nsam2)
582+0.966167450E-1,
+0.418221593,
+0.310454369,
+0.143491149,
+0.200273991,
+0.994065523,
+0.173179626,
+0.617039740,
+0.273689747
586call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
588+0.953835845E-1,
+1.23346829
590iweight1
= getUnifRand(
1,
9, nsam1)
594statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
597call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
599+0.102254748E-1,
+1.62419486
600statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
603call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
605+0.187837362,
+1.08731210
607iweight2
= getUnifRand(
1,
9, nsam2)
609+2,
+7,
+3,
+5,
+6,
+5,
+7,
+8,
+7
611statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
614call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
616+0.922977924E-3,
+1.95972896
617call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
619+0.203965902,
+1.06814909
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.305891573,
+0.968000233
627nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
628sample1
= getUnifRand(
0.,
1., nsam1)
630+0.844657838,
+0.735962570,
+0.884872198,
+0.698993862,
+0.885897815,
+0.731908679,
+0.963922083,
+0.901706517,
+0.979204834
631sample2
= getUnifRand(
0.,
1., nsam2)
633+0.851076305,
+0.826614678,
+0.232144594,
+0.123154879,
+0.694765568,
+0.588098526,
+0.448798299,
+0.378759682,
+0.544264019
637call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
639+0.345075130E-2,
+1.78358042
641iweight1
= getUnifRand(
1,
9, nsam1)
643+3,
+7,
+4,
+5,
+4,
+1,
+9,
+2,
+9
645statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
648call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
650+0.796914101E-4,
+2.25064230
651statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
654call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
656+0.748175383E-2,
+1.67159152
658iweight2
= getUnifRand(
1,
9, nsam2)
660+1,
+3,
+5,
+8,
+7,
+2,
+1,
+9,
+8
662statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
665call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
667+0.00000000,
+4.39442539
668call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
670+0.275969505E-4,
+2.36552072
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.256407261E-2,
+1.82473505
678nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
679sample1
= getUnifRand(
0.,
1., nsam1)
681+0.379793465,
+0.812241614,
+0.256040931,
+0.144738615,
+0.328287482E-1,
+0.849261522,
+0.796898484,
+0.116947114,
+0.459864974,
+0.398342073
682sample2
= getUnifRand(
0.,
1., nsam2)
688call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
690+0.118398249,
+1.18883157
692iweight1
= getUnifRand(
1,
9, nsam1)
694+4,
+5,
+3,
+1,
+1,
+4,
+8,
+7,
+6,
+1
696statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
699call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
701+0.102347910,
+1.21909606
702statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
705call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
707+0.126165390,
+1.17538452
709iweight2
= getUnifRand(
1,
9, nsam2)
713statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
716call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
718+0.203108788E-2,
+1.85638511
719call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
721+0.837290287E-2,
+1.65467572
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.126165390,
+1.17538452
729nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
730sample1
= getUnifRand(
0.,
1., nsam1)
732+0.282666266,
+0.191975117,
+0.651987314,
+0.157700896,
+0.197813511E-1,
+0.118347228,
+0.544796288,
+0.972208917
733sample2
= getUnifRand(
0.,
1., nsam2)
735+0.970036328,
+0.266996980,
+0.781308711,
+0.437777638,
+0.752282739,
+0.591979980
739call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
741+0.253729820,
+1.01552355
743iweight1
= getUnifRand(
1,
9, nsam1)
745+4,
+7,
+5,
+9,
+2,
+2,
+4,
+3
747statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
750call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
752+0.512637496E-1,
+1.35349572
753statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
756call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
758+0.196088552,
+1.07735574
760iweight2
= getUnifRand(
1,
9, nsam2)
762+5,
+2,
+9,
+6,
+1,
+9
764statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
767call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
769+0.345706940E-5,
+2.57537627
770call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
772+0.232598186E-1,
+1.49234259
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.182865977,
+1.09347677
780nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
781sample1
= getUnifRand(
0.,
1., nsam1)
784sample2
= getUnifRand(
0.,
1., nsam2)
786+0.486192226,
+0.219495058
790call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
792+0.201313019,
+1.07121849
794iweight1
= getUnifRand(
1,
9, nsam1)
798statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
801call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
803+0.289170742E-1,
+1.45541525
804statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
807call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
809+0.201313019,
+1.07121849
811iweight2
= getUnifRand(
1,
9, nsam2)
815statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
818call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
820+0.206232071E-4,
+2.39613128
821call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
823+0.109576046,
+1.20501435
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.201651037,
+1.07082570
831nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
832sample1
= getUnifRand(
0.,
1., nsam1)
834+0.851658583E-1,
+0.520192385,
+0.764067113,
+0.431267858,
+0.630271137,
+0.745487034,
+0.463114023,
+0.148182452,
+0.201779604E-2,
+0.153291523
835sample2
= getUnifRand(
0.,
1., nsam2)
837+0.637083530,
+0.636703074,
+0.426450431,
+0.485931754
841call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
843+0.626907468,
+0.750154197
845iweight1
= getUnifRand(
1,
9, nsam1)
847+6,
+1,
+2,
+2,
+3,
+1,
+6,
+7,
+8,
+7
849statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
852call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
854+0.491483212E-1,
+1.36125767
855statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
858call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
860+0.129603982,
+1.16964662
862iweight2
= getUnifRand(
1,
9, nsam2)
866statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
869call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
871+0.226497650E-5,
+2.61731172
872call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
874+0.910007954E-2,
+1.64204466
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.179143071,
+1.09817970
882nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
883sample1
= getUnifRand(
0.,
1., nsam1)
885+0.361340523,
+0.894375920,
+0.824879825,
+0.315386474
886sample2
= getUnifRand(
0.,
1., nsam2)
888+0.905088484,
+0.450432897E-1,
+0.668618262,
+0.639757037,
+0.980680585,
+0.395344317
892call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
894+0.889526010,
+0.580066025
896iweight1
= getUnifRand(
1,
9, nsam1)
900statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
903call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
905+0.249995828,
+1.01918888
906statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
909call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
911+0.607688308,
+0.761593521
913iweight2
= getUnifRand(
1,
9, nsam2)
915+4,
+8,
+3,
+3,
+1,
+4
917statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
920call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
922+0.136151791,
+1.15905488
923call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
925+0.729758561,
+0.688819349
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.902596772,
+0.568898559
933nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
934sample1
= getUnifRand(
0.,
1., nsam1)
937sample2
= getUnifRand(
0.,
1., nsam2)
943call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
945+0.289041579,
+0.982670248
947iweight1
= getUnifRand(
1,
9, nsam1)
951statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
954call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
956+0.167768300,
+1.11304247
957statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
960call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
962+0.289041579,
+0.982670248
964iweight2
= getUnifRand(
1,
9, nsam2)
968statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
971call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
973+0.110656619E-1,
+1.61199534
974call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
976+0.132670283,
+1.16463375
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.289041579,
+0.982670248
984nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
985sample1
= getUnifRand(
0.,
1., nsam1)
987+0.955966115,
+0.847432613,
+0.900346220,
+0.280714631E-1,
+0.952424765,
+0.409191012,
+0.936661601
988sample2
= getUnifRand(
0.,
1., nsam2)
990+0.474887848,
+0.662188649,
+0.288487136,
+0.309700966,
+0.503137946,
+0.810715914,
+0.307850122,
+0.143339992,
+0.143366098
994call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
996+0.171362162E-1,
+1.54267728
998iweight1
= getUnifRand(
1,
9, nsam1)
1000+1,
+5,
+6,
+6,
+7,
+6,
+9
1002statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1005call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
1007+0.620365143E-3,
+2.00977421
1008statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1011call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1013+0.287149549E-1,
+1.45661974
1015iweight2
= getUnifRand(
1,
9, nsam2)
1017+9,
+8,
+7,
+3,
+6,
+1,
+4,
+2,
+1
1019statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1022call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
1024+0.00000000,
+3.25087214
1025call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1027+0.516927242E-2,
+1.72600257
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.482869744E-1,
+1.36450100
1035nsam1
= getUnifRand(
1,
10); nsam2
= getUnifRand(
1,
10)
1036sample1
= getUnifRand(
0.,
1., nsam1)
1038+0.110948205,
+0.812049687,
+0.451585591,
+0.751325905,
+0.986623645,
+0.809446633,
+0.257319748,
+0.804589391
1039sample2
= getUnifRand(
0.,
1., nsam2)
1041+0.927533746,
+0.983831167,
+0.336534858,
+0.550582349
1045call setProbKS(probKS, quanKS, statKS, nsam1, nsam2)
1047+0.739918470,
+0.682632804
1049iweight1
= getUnifRand(
1,
9, nsam1)
1051+6,
+6,
+1,
+6,
+3,
+8,
+4,
+1
1053statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2)
1056call setProbKS(probKS, quanKS, statKS,
sum(iweight1), nsam2)
1058+0.452201903,
+0.858698130
1059statKS
= getDisKolm(sample1, rweight1,
sum(rweight1), sample2)
1062call setProbKS(probKS, quanKS, statKS,
sum(rweight1), nsam2,
sum(rweight1
**2))
1064+0.670893669,
+0.724063098
1066iweight2
= getUnifRand(
1,
9, nsam2)
1070statKS
= getDisKolm(sample1, iweight1,
sum(iweight1), sample2, iweight2,
sum(iweight2))
1073call setProbKS(probKS, quanKS, statKS,
sum(iweight1),
sum(iweight2))
1075+0.892077684E-1,
+1.24696553
1076call setProbKS(probKS, quanKS, statKS,
sum(rweight1),
sum(iweight2),
sum(rweight1
**2))
1078+0.624187946,
+0.751770079
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.935733438,
+0.536488414