R help - Feb 2011 - Hartley's table

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

On 2011-02-16 06:08, Silvano wrote:
> Hi,
>
> I used the commands below to make Hartley's table,
> but some values are NA.
You might want to contact the maintainer of SuppDists (cc'd).
Presumably, the C code uses inappropriate starting values
for those cases.

Meanwhile, here's a simpler way to generate your table:

  for(i in seq_len(ncol(har)))
    har[, i] <- qmaxFratio(.95, 2:(nrow(har) + 1), i+1)


Peter Ehlers
>
>   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):
[... snip ...]