[R] Hartley's table

[Silvano]

Hi,

I used the commands below to make Hartley's table,
but some values are NA.

 require(SuppDists)
 trat = seq(2, 15, 1)
 gl = seq(2, 40, 1)

har = matrix(0, nr=length(gl), nc=length(trat))

for(i in 1:length(gl))
 for(j in 1:length(trat))
 har[i,j] <- qmaxFratio(.95, df=gl[i], k=trat[j])

rownames(har) <- gl
colnames(har) <- trat

head(har)


The output (head):

          2         3          4          5         6 
7
2 39.000000 87.488567 142.502114 202.378010 266.16460 
333.18651
3 15.439182 27.758294  39.503063  50.885084  46.75297 
72.83358
4  9.604530 15.457596  20.559223  25.211423  29.54387 
33.62982
5  7.146382 10.751785  13.723953         NA  10.14968 
20.87902
6  5.819757  8.362843  10.380280  12.108103  13.64260 
15.03555
7  4.994909  6.939901   8.439993   9.697305  10.80480 
11.79569
          8          9        10        11        12 
13
2 403.07945 475.372389 549.84302 626.22781 704.41272 
784.22483
3  83.47794  93.943236 104.24551 114.40008 124.41965 
80.84930
4  37.51656  41.237201  44.81356  48.26765  51.61274 
54.86066
5        NA   2.173131  26.64526        NA        NA 
31.62688
6  16.31888  17.514001  18.63614  19.69659  20.70403 
21.66528
7  12.69999  13.535307  14.31403  15.04529  15.73600 
16.39159
         14        15
2 865.53893 948.24795
3  82.18823  83.40749
4  58.02112  61.10224
5  19.09581  34.64536
6  22.58582  23.46946
7  17.01639  17.61392

What's wrong?

Thanks.

--------------------------------------
Silvano Cesar da Costa
Departamento de Estatística
Universidade Estadual de Londrina
Fone: 3371-4346