tukey-test-tables.a68
1 COMMENT
2
3 @section Synopsis
4
5 Compute Studentized range q table.
6
7 This program requires a68g to be built with R mathlib.
8
9 Sample output:
10
11 Confidence = 0.050
12 k->
13 df 2 3 4 5 6 7 8 9 10
14 2 6.080 8.331 9.799 10.881 11.734 12.435 13.028 13.542 13.994
15 3 4.501 5.910 6.825 7.502 8.037 8.478 8.852 9.177 9.462
16 4 3.927 5.040 5.757 6.287 6.706 7.053 7.347 7.602 7.826
17 5 3.635 4.602 5.218 5.673 6.033 6.330 6.582 6.801 6.995
18 6 3.460 4.339 4.896 5.305 5.628 5.895 6.122 6.319 6.493
19 7 3.344 4.165 4.681 5.060 5.359 5.606 5.815 5.997 6.158
20 8 3.261 4.041 4.529 4.886 5.167 5.399 5.596 5.767 5.918
21 9 3.199 3.948 4.415 4.755 5.024 5.244 5.432 5.595 5.738
22 10 3.151 3.877 4.327 4.654 4.912 5.124 5.304 5.460 5.598
23 11 3.113 3.820 4.256 4.574 4.823 5.028 5.202 5.353 5.486
24 12 3.081 3.773 4.199 4.508 4.750 4.950 5.119 5.265 5.395
25 13 3.055 3.734 4.151 4.453 4.690 4.884 5.049 5.192 5.318
26 14 3.033 3.701 4.111 4.407 4.639 4.829 4.990 5.130 5.253
27 15 3.014 3.673 4.076 4.367 4.595 4.782 4.940 5.077 5.198
28 16 2.998 3.649 4.046 4.333 4.557 4.741 4.896 5.031 5.150
29 17 2.984 3.628 4.020 4.303 4.524 4.705 4.858 4.991 5.108
30 18 2.971 3.609 3.997 4.276 4.494 4.673 4.824 4.955 5.071
31 19 2.960 3.593 3.977 4.253 4.468 4.645 4.794 4.924 5.037
32 20 2.950 3.578 3.958 4.232 4.445 4.620 4.768 4.895 5.008
33
34 COMMENT
35
36 PR need mathlib PR
37
38 PROC table = (REAL alpha) VOID:
39 BEGIN printf (($"Confidence = ", d.dddl4x"k->"l"df"$, 1 - alpha));
40 FOR k FROM 2 TO 10
41 DO printf (($x6zd$, k))
42 OD;
43 print (new line);
44 FOR df FROM 2 TO 20
45 DO print (whole (df, -2));
46 FOR k FROM 2 TO 10
47 DO printf (($xzzd.ddd$, r qtukey (alpha, 1, k, df, TRUE, FALSE)))
48 OD;
49 print (new line)
50 OD
51 END;
52
53 table (0.9);
54 print (new line);
55 table (0.95);
56 print (new line);
57 table (0.99)
