##
## Examples: Completely Randomized Design (CRD)
## More details: demo(package='TukeyC')
##
library(TukeyC)
data(CRD2)
## From: formula
tk1 <- with(CRD2,
TukeyC(y ~ x,
data=dfm,
which='x'))
old.par <- par(mar=c(6, 3, 6, 2))
plot(tk1,
id.las=2)
plot(tk1,
yl=FALSE,
disp='sd',
id.las=2)
## From: aov
av <- with(CRD2,
aov(y ~ x,
data=dfm))
summary(av)
tk2 <- TukeyC(x=av,
which='x')
plot(tk2,
disp='sd',
yl=FALSE,
id.las=2)
# From: lm
av_lm <- with(CRD2,
lm(y ~ x,
data=dfm))
tk3 <- TukeyC(x=av_lm,
which='x')
par(mfrow=c(2, 1))
plot(tk3,
disp='ci',
id.las=2,
yl=FALSE)
plot(tk3,
disp='cip',
id.las=2,
yl=FALSE)
par(mfrow=c(1, 1))
par(old.par)
Run the code above in your browser using DataLab