I have been trying to figure out how to run a simple simulation of the ANOVA
and I'm coming up just a bit short. The code I've got is:

cohen.f = .25
groups = 4
between.var = 19
within.var =  between.var / cohen.f ^ 2
n = 500
sim.means = rnorm(n = groups, mean = 0, sd = sqrt(between.var))
sim.data = lapply(sim.means, function(mean) rnorm(n = n, mean = mean, sd =
sqrt(within.var)))
sim.anova = ?

I'm pretty sure I've got data with the between group and within group
variances that I want. I don't really know how to run an ANOVA on it though.
All the examples I've been able to find (with aov and anova) require an
explicit model. The ANOVA just tests the equivalence of means for a set of
groups, right? Can I not just tell it to give me the likelihood the means
are equal for all the columns without writing an explicit model.

Any help would be appreciated,
Will

	[[alternative HTML version deleted]]

