Nothing
RManovaNov<-function(r,k,f,n){
tm = gl(k, 1, n*k, unique(factor(f))) # matching treatment
tm
blk = gl(n, k, k*n) # blocking factor
blk
av = aov(r ~ blk+tm)
summary(av)
qqnorm(av$residuals)
plot(av$fitted.values,av$residuals)
TukeyHSD(av)$tm
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.