Algol 68 Genie

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)
Deze website is gearchiveerd door de KB, nationale bibliotheek.

This website is archived by the National Library of the Netherlands.