Displaying 1 result from an estimated 1 matches for "tukeymc".
Did you mean:
tukey
2000 Jan 04
0
multiple comparison function
...# na, nb = total no. of data per
# level of factor tested
# (or sample sizes for one-
# way)
# df = within-cells degrees of
# freedom (or overall df for one-
# way)
tukeymc <- function(x, ms, na, nb, df) {
x <- sort(x) # rank means
se <- sqrt(ms/2 * (1/na + 1/nb))
#
# Generate index vectors a and b for
# sequential subtraction of means
#
y <- rev(1:length(x))
y <- y[-length(x)]
a <- rep(y, times = y-1)
mb...