[R] Tukey's table

Silvano silvano at uel.br
Wed Nov 3 11:52:34 CET 2010


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="")


Some suggest?

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

More information about the R-help mailing list