[R] Hartley's table

[Peter Ehlers]

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 ...]