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