data <- data.frame("index" = 1:50)
names <- c()
for (i in 1:20) {
	names <- c(names, paste("someReallyLongVariableName", i, sep = ""))
	data[[names[i]]] <- rnorm(50)
}
names <- factor(names)
model <- as.formula(paste("cbind(", paste(names, collapse = ","), ") ~ 1 + Error(1)", sep = ""))
result <- aov(model, data)
