# NOT RUN {
data(rimov)
groups <- factor(c(rep(1, times=12), rep(2, times=12), rep(3, times=12)))
# Plot data in groups
subs <- function(group, ...) {
cset <- rimov
cset$obs <- rimov$obs[, groups == group]
plot(cset, ...)
}
for(i in 1:3) subs(i, main=paste("group ", i, sep=""), ylab="Temperature")
# See example analysis in ?graph.fanova
# }
Run the code above in your browser using DataLab