[R] Tukey's table

[Silvano]

Hi,

I'm building Tukey's table using qtukey function.

It happens that I can't get the values of Tukey's one degree 
of freedom and also wanted to eliminate the first column.

The program is:

Trat <- c(1:30)                     # number of treatments
gl <- c(1:30, 40, 60, 120)     # degree freedom

tukval <- matrix(0, nr=length(gl), nc=length(Trat))

for(i in 1:length(gl))
  for(j in 1:length(Trat))
    tukval[i,j] <- qtukey(.95, Trat[j], gl[i])

rownames(tukval) <- gl
colnames(tukval) <- paste(Trat, "", sep="")
tukval

require(xtable)
xtable(tukval)


Some suggest?

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