# NOT RUN {
## This example is based on chunks 1, 2, 4, 6 of HH::HHscriptnames(13)
## It defines the data(cc176.y.adj) that is used in
## demo("bwplot-lattice") ## 5 calls to latex() and 5 calls to msWord()
## demo("boxplot-ggplot") ## 2 calls to latex() and 3 calls to msWord()
data(cc176, package="HH")
cc176.aov <- aov(wt.d ~ rep + wt.n + n.treats*minutes*current,
data=cc176)
cc176.y.adj <- cc176$wt.d -
(cc176$wt.n - mean(cc176$wt.n))*coef(cc176.aov)["wt.n"]
tmp <-
sapply(split(cc176.y.adj, cc176$current),
function(x)
c(min=min(x),
"m-sd"=mean(x)-sd(x),
mean=mean(x),
"m+sd"=mean(x)+sd(x),
max=max(x)))
cc176fivenumsd <- t(tmp)
save(cc176.y.adj, cc176fivenumsd, file="cc176.y.adj.rda")
# }
Run the code above in your browser using DataLab