Hi:
Try this:
Trat <- c(2:30) # number of treatments
gl <- c(2:30, 40, 60, 120)
# Write a one-line 2D function to get the Tukey distribution quantile:
f <- function(x,y) qtukey(0.95, x, y)
outer(Trat, gl, f)
It's slow (takes a few seconds) but it seems to work.
HTH,
Dennis
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 free
2 matches
Mail list logo