GRASS GIS 7 Programmer's Manual
7.8.3(2020)-exported
shapiro1.c
Go to the documentation of this file.
1
#include <stdio.h>
2
#include <stdlib.h>
3
#include <math.h>
4
#include "local_proto.h"
5
6
7
double
*
Cdhc_shapiro_wilk
(
double
*
x
,
int
n)
8
{
9
static
double
y[2];
10
double
a[25], s2, *xcopy;
11
double
sumb = 0.0, sumx = 0.0, sumx2 = 0.0;
12
int
i, k;
13
14
if
((xcopy = (
double
*)malloc(n *
sizeof
(
double
))) ==
NULL
) {
15
fprintf(stderr,
"Memory error in Cdhc_shapiro_wilk\n"
);
16
exit(EXIT_FAILURE);
17
}
18
19
k = n / 2;
20
for
(i = 0; i < n; ++i) {
21
xcopy[i] =
x
[i];
22
sumx +=
x
[i];
23
sumx2 +=
x
[i] *
x
[i];
24
}
25
s2 = sumx2 - sumx * sumx / n;
26
27
qsort(xcopy, n,
sizeof
(
double
),
Cdhc_dcmp
);
28
29
if
(n == 3) {
30
a[0] = (double).7071;
31
}
32
else
if
(n == 4) {
33
a[0] = (double).6872;
34
a[1] = (double).1677;
35
}
36
else
if
(n == 5) {
37
a[0] = (double).6646;
38
a[1] = (double).2413;
39
}
40
else
if
(n == 6) {
41
a[0] = (double).6431;
42
a[1] = (double).2806;
43
a[2] = (double).0875;
44
}
45
else
if
(n == 7) {
46
a[0] = (double).6233;
47
a[1] = (double).3031;
48
a[2] = (double).1401;
49
}
50
else
if
(n == 8) {
51
a[0] = (double).6052;
52
a[1] = (double).3164;
53
a[2] = (double).1743;
54
a[3] = (double).0561;
55
}
56
else
if
(n == 9) {
57
a[0] = (double).5888;
58
a[1] = (double).3244;
59
a[2] = (double).1976;
60
a[3] = (double).0947;
61
}
62
else
if
(n == 10) {
63
a[0] = (double).5739;
64
a[1] = (double).3291;
65
a[2] = (double).2141;
66
a[3] = (double).1224;
67
a[4] = (double).0399;
68
}
69
else
if
(n == 11) {
70
a[0] = (double).5601;
71
a[1] = (double).3315;
72
a[2] = (double).226;
73
a[3] = (double).1429;
74
a[4] = (double).0695;
75
}
76
else
if
(n == 12) {
77
a[0] = (double).5475;
78
a[1] = (double).3325;
79
a[2] = (double).2347;
80
a[3] = (double).1586;
81
a[4] = (double).0922;
82
a[5] = (double).0303;
83
}
84
else
if
(n == 13) {
85
a[0] = (double).5359;
86
a[1] = (double).3325;
87
a[2] = (double).2412;
88
a[3] = (double).1707;
89
a[4] = (double).1099;
90
a[5] = (double).0539;
91
}
92
else
if
(n == 14) {
93
a[0] = (double).5251;
94
a[1] = (double).3318;
95
a[2] = (double).246;
96
a[3] = (double).1802;
97
a[4] = (double).124;
98
a[5] = (double).0727;
99
a[6] = (double).024;
100
}
101
else
if
(n == 15) {
102
a[0] = (double).515;
103
a[1] = (double).3306;
104
a[2] = (double).2495;
105
a[3] = (double).1878;
106
a[4] = (double).1353;
107
a[5] = (double).088;
108
a[6] = (double).0433;
109
}
110
else
if
(n == 16) {
111
a[0] = (double).5056;
112
a[1] = (double).329;
113
a[2] = (double).2521;
114
a[3] = (double).1939;
115
a[4] = (double).1447;
116
a[5] = (double).1005;
117
a[6] = (double).0593;
118
a[7] = (double).0196;
119
}
120
else
if
(n == 17) {
121
a[0] = (double).4968;
122
a[1] = (double).3273;
123
a[2] = (double).254;
124
a[3] = (double).1988;
125
a[4] = (double).1524;
126
a[5] = (double).1109;
127
a[6] = (double).0725;
128
a[7] = (double).0359;
129
}
130
else
if
(n == 18) {
131
a[0] = (double).4886;
132
a[1] = (double).3253;
133
a[2] = (double).2553;
134
a[3] = (double).2027;
135
a[4] = (double).1587;
136
a[5] = (double).1197;
137
a[6] = (double).0837;
138
a[7] = (double).0496;
139
a[8] = (double).0163;
140
}
141
else
if
(n == 19) {
142
a[0] = (double).4808;
143
a[1] = (double).3232;
144
a[2] = (double).2561;
145
a[3] = (double).2059;
146
a[4] = (double).1641;
147
a[5] = (double).1271;
148
a[6] = (double).0932;
149
a[7] = (double).0612;
150
a[8] = (double).0303;
151
}
152
else
if
(n == 20) {
153
a[0] = (double).4734;
154
a[1] = (double).3211;
155
a[2] = (double).2565;
156
a[3] = (double).2085;
157
a[4] = (double).1686;
158
a[5] = (double).1334;
159
a[6] = (double).1013;
160
a[7] = (double).0711;
161
a[8] = (double).0422;
162
a[9] = (double).014;
163
}
164
else
if
(n == 21) {
165
a[0] = (double).4643;
166
a[1] = (double).3185;
167
a[2] = (double).2578;
168
a[3] = (double).2119;
169
a[4] = (double).1736;
170
a[5] = (double).1399;
171
a[6] = (double).1092;
172
a[7] = (double).0804;
173
a[8] = (double).053;
174
a[9] = (double).0263;
175
}
176
else
if
(n == 22) {
177
a[0] = (double).459;
178
a[1] = (double).3156;
179
a[2] = (double).2571;
180
a[3] = (double).2131;
181
a[4] = (double).1764;
182
a[5] = (double).1443;
183
a[6] = (double).115;
184
a[7] = (double).0878;
185
a[8] = (double).0618;
186
a[9] = (double).0368;
187
a[10] = (double).0122;
188
}
189
else
if
(n == 23) {
190
a[0] = (double).4542;
191
a[1] = (double).3126;
192
a[2] = (double).2563;
193
a[3] = (double).2139;
194
a[4] = (double).1787;
195
a[5] = (double).148;
196
a[6] = (double).1201;
197
a[7] = (double).0941;
198
a[8] = (double).0696;
199
a[9] = (double).0459;
200
a[10] = (double).0228;
201
}
202
else
if
(n == 24) {
203
a[0] = (double).4493;
204
a[1] = (double).3098;
205
a[2] = (double).2554;
206
a[3] = (double).2145;
207
a[4] = (double).1807;
208
a[5] = (double).1512;
209
a[6] = (double).1245;
210
a[7] = (double).0997;
211
a[8] = (double).0764;
212
a[9] = (double).0539;
213
a[10] = (double).0321;
214
a[11] = (double).0107;
215
}
216
else
if
(n == 25) {
217
a[0] = (double).445;
218
a[1] = (double).3069;
219
a[2] = (double).2543;
220
a[3] = (double).2148;
221
a[4] = (double).1822;
222
a[5] = (double).1539;
223
a[6] = (double).1283;
224
a[7] = (double).1046;
225
a[8] = (double).0823;
226
a[9] = (double).061;
227
a[10] = (double).0403;
228
a[11] = (double).02;
229
}
230
else
if
(n == 26) {
231
a[0] = (double).4407;
232
a[1] = (double).3043;
233
a[2] = (double).2533;
234
a[3] = (double).2151;
235
a[4] = (double).1836;
236
a[5] = (double).1563;
237
a[6] = (double).1316;
238
a[7] = (double).1089;
239
a[8] = (double).0876;
240
a[9] = (double).0672;
241
a[10] = (double).0476;
242
a[11] = (double).0284;
243
a[12] = (double).0094;
244
}
245
else
if
(n == 27) {
246
a[0] = (double).4366;
247
a[1] = (double).3018;
248
a[2] = (double).2522;
249
a[3] = (double).2152;
250
a[4] = (double).1848;
251
a[5] = (double).1584;
252
a[6] = (double).1346;
253
a[7] = (double).1128;
254
a[8] = (double).0923;
255
a[9] = (double).0728;
256
a[10] = (double).054;
257
a[11] = (double).0358;
258
a[12] = (double).0178;
259
}
260
else
if
(n == 28) {
261
a[0] = (double).4328;
262
a[1] = (double).2992;
263
a[2] = (double).251;
264
a[3] = (double).2151;
265
a[4] = (double).1857;
266
a[5] = (double).1601;
267
a[6] = (double).1372;
268
a[7] = (double).1162;
269
a[8] = (double).0965;
270
a[9] = (double).0778;
271
a[10] = (double).0598;
272
a[11] = (double).0424;
273
a[12] = (double).0253;
274
a[13] = (double).0084;
275
}
276
else
if
(n == 29) {
277
a[0] = (double).4291;
278
a[1] = (double).2968;
279
a[2] = (double).2499;
280
a[3] = (double).215;
281
a[4] = (double).1864;
282
a[5] = (double).1616;
283
a[6] = (double).1395;
284
a[7] = (double).1192;
285
a[8] = (double).1002;
286
a[9] = (double).0822;
287
a[10] = (double).065;
288
a[11] = (double).0483;
289
a[12] = (double).032;
290
a[13] = (double).0159;
291
}
292
else
if
(n == 30) {
293
a[0] = (double).4254;
294
a[1] = (double).2944;
295
a[2] = (double).2487;
296
a[3] = (double).2148;
297
a[4] = (double).187;
298
a[5] = (double).163;
299
a[6] = (double).1415;
300
a[7] = (double).1219;
301
a[8] = (double).1036;
302
a[9] = (double).0862;
303
a[10] = (double).0697;
304
a[11] = (double).0537;
305
a[12] = (double).0381;
306
a[13] = (double).0227;
307
a[14] = (double).0076;
308
}
309
else
if
(n == 31) {
310
a[0] = (double).422;
311
a[1] = (double).2921;
312
a[2] = (double).2475;
313
a[3] = (double).2145;
314
a[4] = (double).1874;
315
a[5] = (double).1641;
316
a[6] = (double).1433;
317
a[7] = (double).1243;
318
a[8] = (double).1066;
319
a[9] = (double).0899;
320
a[10] = (double).0739;
321
a[11] = (double).0585;
322
a[12] = (double).0435;
323
a[13] = (double).0289;
324
a[14] = (double).0144;
325
}
326
else
if
(n == 32) {
327
a[0] = (double).4188;
328
a[1] = (double).2898;
329
a[2] = (double).2463;
330
a[3] = (double).2141;
331
a[4] = (double).1878;
332
a[5] = (double).1651;
333
a[6] = (double).1449;
334
a[7] = (double).1265;
335
a[8] = (double).1093;
336
a[9] = (double).0931;
337
a[10] = (double).0777;
338
a[11] = (double).0629;
339
a[12] = (double).0485;
340
a[13] = (double).0344;
341
a[14] = (double).0206;
342
a[15] = (double).0068;
343
}
344
else
if
(n == 33) {
345
a[0] = (double).4156;
346
a[1] = (double).2876;
347
a[2] = (double).2451;
348
a[3] = (double).2137;
349
a[4] = (double).188;
350
a[5] = (double).166;
351
a[6] = (double).1463;
352
a[7] = (double).1284;
353
a[8] = (double).1118;
354
a[9] = (double).0961;
355
a[10] = (double).0812;
356
a[11] = (double).0669;
357
a[12] = (double).053;
358
a[13] = (double).0395;
359
a[14] = (double).0262;
360
a[15] = (double).0131;
361
}
362
else
if
(n == 34) {
363
a[0] = (double).4127;
364
a[1] = (double).2854;
365
a[2] = (double).2439;
366
a[3] = (double).2132;
367
a[4] = (double).1882;
368
a[5] = (double).1667;
369
a[6] = (double).1475;
370
a[7] = (double).1301;
371
a[8] = (double).114;
372
a[9] = (double).0988;
373
a[10] = (double).0844;
374
a[11] = (double).0706;
375
a[12] = (double).0572;
376
a[13] = (double).0441;
377
a[14] = (double).0314;
378
a[15] = (double).0187;
379
a[16] = (double).0062;
380
}
381
else
if
(n == 35) {
382
a[0] = (double).4096;
383
a[1] = (double).2834;
384
a[2] = (double).2427;
385
a[3] = (double).2127;
386
a[4] = (double).1883;
387
a[5] = (double).1673;
388
a[6] = (double).1487;
389
a[7] = (double).1317;
390
a[8] = (double).116;
391
a[9] = (double).1013;
392
a[10] = (double).0873;
393
a[11] = (double).0739;
394
a[12] = (double).061;
395
a[13] = (double).0484;
396
a[14] = (double).0361;
397
a[15] = (double).0239;
398
a[16] = (double).0119;
399
}
400
else
if
(n == 36) {
401
a[0] = (double).4068;
402
a[1] = (double).2813;
403
a[2] = (double).2415;
404
a[3] = (double).2121;
405
a[4] = (double).1883;
406
a[5] = (double).1678;
407
a[6] = (double).1496;
408
a[7] = (double).1331;
409
a[8] = (double).1179;
410
a[9] = (double).1036;
411
a[10] = (double).09;
412
a[11] = (double).077;
413
a[12] = (double).0645;
414
a[13] = (double).0523;
415
a[14] = (double).0404;
416
a[15] = (double).0287;
417
a[16] = (double).0172;
418
a[17] = (double).0057;
419
}
420
else
if
(n == 37) {
421
a[0] = (double).404;
422
a[1] = (double).2794;
423
a[2] = (double).2403;
424
a[3] = (double).2116;
425
a[4] = (double).1883;
426
a[5] = (double).1683;
427
a[6] = (double).1505;
428
a[7] = (double).1344;
429
a[8] = (double).1196;
430
a[9] = (double).1056;
431
a[10] = (double).0924;
432
a[11] = (double).0798;
433
a[12] = (double).0677;
434
a[13] = (double).0559;
435
a[14] = (double).0444;
436
a[15] = (double).0331;
437
a[16] = (double).022;
438
a[17] = (double).011;
439
}
440
else
if
(n == 38) {
441
a[0] = (double).4015;
442
a[1] = (double).2774;
443
a[2] = (double).2391;
444
a[3] = (double).211;
445
a[4] = (double).1881;
446
a[5] = (double).1686;
447
a[6] = (double).1513;
448
a[7] = (double).1356;
449
a[8] = (double).1211;
450
a[9] = (double).1075;
451
a[10] = (double).0947;
452
a[11] = (double).0824;
453
a[12] = (double).0706;
454
a[13] = (double).0592;
455
a[14] = (double).0481;
456
a[15] = (double).0372;
457
a[16] = (double).0264;
458
a[17] = (double).0158;
459
a[18] = (double).0053;
460
}
461
else
if
(n == 39) {
462
a[0] = (double).3989;
463
a[1] = (double).2755;
464
a[2] = (double).238;
465
a[3] = (double).2104;
466
a[4] = (double).188;
467
a[5] = (double).1689;
468
a[6] = (double).152;
469
a[7] = (double).1366;
470
a[8] = (double).1225;
471
a[9] = (double).1092;
472
a[10] = (double).0967;
473
a[11] = (double).0848;
474
a[12] = (double).0733;
475
a[13] = (double).0622;
476
a[14] = (double).0515;
477
a[15] = (double).0409;
478
a[16] = (double).0305;
479
a[17] = (double).0203;
480
a[18] = (double).0101;
481
}
482
else
if
(n == 40) {
483
a[0] = (double).3964;
484
a[1] = (double).2737;
485
a[2] = (double).2368;
486
a[3] = (double).2098;
487
a[4] = (double).1878;
488
a[5] = (double).1691;
489
a[6] = (double).1526;
490
a[7] = (double).1376;
491
a[8] = (double).1237;
492
a[9] = (double).1108;
493
a[10] = (double).0986;
494
a[11] = (double).087;
495
a[12] = (double).0759;
496
a[13] = (double).0651;
497
a[14] = (double).0546;
498
a[15] = (double).0444;
499
a[16] = (double).0343;
500
a[17] = (double).0244;
501
a[18] = (double).0146;
502
a[19] = (double).0049;
503
}
504
else
if
(n == 41) {
505
a[0] = (double).394;
506
a[1] = (double).2719;
507
a[2] = (double).2357;
508
a[3] = (double).2091;
509
a[4] = (double).1876;
510
a[5] = (double).1693;
511
a[6] = (double).1531;
512
a[7] = (double).1384;
513
a[8] = (double).1249;
514
a[9] = (double).1123;
515
a[10] = (double).1004;
516
a[11] = (double).0891;
517
a[12] = (double).0782;
518
a[13] = (double).0677;
519
a[14] = (double).0575;
520
a[15] = (double).0476;
521
a[16] = (double).0379;
522
a[17] = (double).0283;
523
a[18] = (double).0188;
524
a[19] = (double).0094;
525
}
526
else
if
(n == 42) {
527
a[0] = (double).3917;
528
a[1] = (double).2701;
529
a[2] = (double).2345;
530
a[3] = (double).2085;
531
a[4] = (double).1874;
532
a[5] = (double).1694;
533
a[6] = (double).1535;
534
a[7] = (double).1392;
535
a[8] = (double).1259;
536
a[9] = (double).1136;
537
a[10] = (double).102;
538
a[11] = (double).0909;
539
a[12] = (double).0804;
540
a[13] = (double).0701;
541
a[14] = (double).0602;
542
a[15] = (double).0506;
543
a[16] = (double).0411;
544
a[17] = (double).0318;
545
a[18] = (double).0227;
546
a[19] = (double).0136;
547
a[20] = (double).0045;
548
}
549
else
if
(n == 43) {
550
a[0] = (double).3894;
551
a[1] = (double).2684;
552
a[2] = (double).2334;
553
a[3] = (double).2078;
554
a[4] = (double).1871;
555
a[5] = (double).1695;
556
a[6] = (double).1539;
557
a[7] = (double).1398;
558
a[8] = (double).1269;
559
a[9] = (double).1149;
560
a[10] = (double).1035;
561
a[11] = (double).0927;
562
a[12] = (double).0824;
563
a[13] = (double).0724;
564
a[14] = (double).0628;
565
a[15] = (double).0534;
566
a[16] = (double).0442;
567
a[17] = (double).0352;
568
a[18] = (double).0263;
569
a[19] = (double).0175;
570
a[20] = (double).0087;
571
}
572
else
if
(n == 44) {
573
a[0] = (double).3872;
574
a[1] = (double).2667;
575
a[2] = (double).2323;
576
a[3] = (double).2072;
577
a[4] = (double).1868;
578
a[5] = (double).1695;
579
a[6] = (double).1542;
580
a[7] = (double).1405;
581
a[8] = (double).1278;
582
a[9] = (double).116;
583
a[10] = (double).1049;
584
a[11] = (double).0943;
585
a[12] = (double).0842;
586
a[13] = (double).0745;
587
a[14] = (double).0651;
588
a[15] = (double).056;
589
a[16] = (double).0471;
590
a[17] = (double).0383;
591
a[18] = (double).0296;
592
a[19] = (double).0211;
593
a[20] = (double).0126;
594
a[21] = (double).0042;
595
}
596
else
if
(n == 45) {
597
a[0] = (double).385;
598
a[1] = (double).2651;
599
a[2] = (double).2313;
600
a[3] = (double).2065;
601
a[4] = (double).1865;
602
a[5] = (double).1695;
603
a[6] = (double).1545;
604
a[7] = (double).141;
605
a[8] = (double).1286;
606
a[9] = (double).117;
607
a[10] = (double).1062;
608
a[11] = (double).0959;
609
a[12] = (double).086;
610
a[13] = (double).0765;
611
a[14] = (double).0673;
612
a[15] = (double).0584;
613
a[16] = (double).0497;
614
a[17] = (double).0412;
615
a[18] = (double).0328;
616
a[19] = (double).0245;
617
a[20] = (double).0163;
618
a[21] = (double).0081;
619
}
620
else
if
(n == 46) {
621
a[0] = (double).383;
622
a[1] = (double).2635;
623
a[2] = (double).2302;
624
a[3] = (double).2058;
625
a[4] = (double).1862;
626
a[5] = (double).1695;
627
a[6] = (double).1548;
628
a[7] = (double).1415;
629
a[8] = (double).1293;
630
a[9] = (double).118;
631
a[10] = (double).1073;
632
a[11] = (double).0972;
633
a[12] = (double).0876;
634
a[13] = (double).0783;
635
a[14] = (double).0694;
636
a[15] = (double).0607;
637
a[16] = (double).0522;
638
a[17] = (double).0439;
639
a[18] = (double).0357;
640
a[19] = (double).0277;
641
a[20] = (double).0197;
642
a[21] = (double).0118;
643
a[22] = (double).0039;
644
}
645
else
if
(n == 47) {
646
a[0] = (double).3808;
647
a[1] = (double).262;
648
a[2] = (double).2291;
649
a[3] = (double).2052;
650
a[4] = (double).1859;
651
a[5] = (double).1695;
652
a[6] = (double).155;
653
a[7] = (double).142;
654
a[8] = (double).13;
655
a[9] = (double).1189;
656
a[10] = (double).1085;
657
a[11] = (double).0986;
658
a[12] = (double).0892;
659
a[13] = (double).0801;
660
a[14] = (double).0713;
661
a[15] = (double).0628;
662
a[16] = (double).0546;
663
a[17] = (double).0465;
664
a[18] = (double).0385;
665
a[19] = (double).0307;
666
a[20] = (double).0229;
667
a[21] = (double).0153;
668
a[22] = (double).0076;
669
}
670
else
if
(n == 48) {
671
a[0] = (double).3789;
672
a[1] = (double).2604;
673
a[2] = (double).2281;
674
a[3] = (double).2045;
675
a[4] = (double).1855;
676
a[5] = (double).1693;
677
a[6] = (double).1551;
678
a[7] = (double).1423;
679
a[8] = (double).1306;
680
a[9] = (double).1197;
681
a[10] = (double).1095;
682
a[11] = (double).0998;
683
a[12] = (double).0906;
684
a[13] = (double).0817;
685
a[14] = (double).0731;
686
a[15] = (double).0648;
687
a[16] = (double).0568;
688
a[17] = (double).0489;
689
a[18] = (double).0411;
690
a[19] = (double).0335;
691
a[20] = (double).0259;
692
a[21] = (double).0185;
693
a[22] = (double).0111;
694
a[23] = (double).0037;
695
}
696
else
if
(n == 49) {
697
a[0] = (double).377;
698
a[1] = (double).2589;
699
a[2] = (double).2271;
700
a[3] = (double).2038;
701
a[4] = (double).1851;
702
a[5] = (double).1692;
703
a[6] = (double).1553;
704
a[7] = (double).1427;
705
a[8] = (double).1312;
706
a[9] = (double).1205;
707
a[10] = (double).1105;
708
a[11] = (double).101;
709
a[12] = (double).0919;
710
a[13] = (double).0832;
711
a[14] = (double).0748;
712
a[15] = (double).0667;
713
a[16] = (double).0588;
714
a[17] = (double).0511;
715
a[18] = (double).0436;
716
a[19] = (double).0361;
717
a[20] = (double).0288;
718
a[21] = (double).0215;
719
a[22] = (double).0143;
720
a[23] = (double).0071;
721
}
722
else
if
(n == 50) {
723
a[0] = (double).3751;
724
a[1] = (double).2574;
725
a[2] = (double).226;
726
a[3] = (double).2032;
727
a[4] = (double).1847;
728
a[5] = (double).1691;
729
a[6] = (double).1554;
730
a[7] = (double).143;
731
a[8] = (double).1317;
732
a[9] = (double).1212;
733
a[10] = (double).1113;
734
a[11] = (double).102;
735
a[12] = (double).0932;
736
a[13] = (double).0846;
737
a[14] = (double).0764;
738
a[15] = (double).0685;
739
a[16] = (double).0608;
740
a[17] = (double).0532;
741
a[18] = (double).0459;
742
a[19] = (double).0386;
743
a[20] = (double).0314;
744
a[21] = (double).0244;
745
a[22] = (double).0174;
746
a[23] = (double).0104;
747
a[24] = (double).0035;
748
}
749
750
if
(n > 50 || n < 3) {
751
#ifdef NOISY
752
fprintf(stdout,
753
" THIS IS THE SHAPIRO-WILK TEST FOR SMALL SAMPLES\n"
);
754
fprintf(stdout,
755
" THE SAMPLE SIZE MUST BE LESS THAN OR EQUAL TO 50\n"
);
756
#endif
/* NOISY */
757
758
y[0] = y[1] = 0.0;
759
}
760
else
{
761
for
(i = 1; i <= k; ++i)
762
sumb += a[i - 1] * (
x
[n - i + 1] -
x
[i]);
763
764
y[0] = sumb * sumb / s2;
765
y[1] = s2;
766
767
#ifdef NOISY
768
fprintf(stdout,
" TEST13 SW(N) =%10.4f\n"
, y[0]);
769
#endif
/* NOISY */
770
}
771
772
free(xcopy);
773
774
return
y;
775
}
x
#define x
Cdhc_shapiro_wilk
double * Cdhc_shapiro_wilk(double *x, int n)
Definition:
shapiro1.c:7
Cdhc_dcmp
int Cdhc_dcmp(const void *i, const void *j)
Definition:
dcmp.c:1
NULL
#define NULL
Definition:
ccmath.h:32
cdhc
shapiro1.c
Generated on Tue May 5 2020 09:45:27 for GRASS GIS 7 Programmer's Manual by
1.8.17