Tengo un gran conjunto de datos que me gustaría realizar cálculos post hoc:

dat = as.data.frame(matrix(runif(10000*300), ncol = 10000, nrow = 300))

dat$group = rep(letters[1:3], 100)

Aquí está mi código:

start <- Sys.time()

vars <- names(dat)[-ncol(dat)] 

aov.out <- lapply(vars, function(x) {
        lm(substitute(i ~ group, list(i = as.name(x))), data = dat)})

TukeyHSD.out <- lapply(aov.out, function(x) TukeyHSD(aov(x)))

Sys.time() - start

Diferencia horaria de 4.033335 minutos

Tarda unos 4 minutos, ¿hay formas más eficientes y elegantes de realizar post hoc usando R?

Muchas gracia