m<-4; n<-5
p <- matrix(c(30,1,32,1,34,1,38,1),byrow=TRUE,ncol=2)
Xr <- invent.mxn(m,n,d=1,p,f2="random")
y <- c(Xr)
f <- factor(rep(LETTERS[1:m], rep(n,m)))
f.y <- data.frame(f, y)
boxplot(y~f, data=f.y,ylab="y", xlab="f")
summary(aov(y~f, data=f.y))
Run the code above in your browser using DataLab