| plot.SK {ScottKnott} | R Documentation |
Those are S3 methods for objects of class SK and SK.nest.
## S3 method for class 'SK':
plot(x, inches=TRUE, col=NULL, xlab=NULL, ylab=NULL,
title="Groups of means divided by colors", ...)
## S3 method for class 'SK.nest':
plot(x, inches=TRUE, col=NULL, xlab=NULL, ylab=NULL,
title="Groups of means divided by colors", ...)
x |
A SK or SK.nest object. |
inches |
'TRUE', 'FALSE' or a positive number.See symbols-package:graghics |
col |
A vector of colors for the thermometers representation. See symbols-package:graghics. |
xlab |
A label for the x axis. |
ylab |
A label for the y axis. |
title |
A title for the graphics. |
... |
Other graphical parameters may also be passed as arguments to this function. |
The plot.SK and plot.SK.nest functions generates a series of
thermometers where the means are marked by a line in the midle, and the
groups are diferentiated by diferent colors.
Enio Jelihovschi (eniojelihovs@gmail.com)
Jose Claudio Faria (joseclaudio.faria@gmail.com)
Sergio Oliveira (solive@uesc.br)
Murrell, P. (2005) R Graphics. Chapman & Hall/CRC Press.
##
## Example: Factorial Experiment (FE)
##
yeld <- c(53, 27, 40, 45, 43, 45, 32, 12,
45, 57, 29, 69, 42, 32, 61, 54,
95, 27, 24, 60, 60, 98, 11, 26)
labs <- factor(rep(1:2, 12))
runs <- factor(sort(rep(1:12, 2)))
pressure <- gl(2, 12)
catalyst <- factor(rep(sort(rep(1:2, 2)), 6))
dm <- data.frame(pressure, catalyst, labs) # Design matrix (a data.frame object)
dfm <- data.frame(labs, pressure, catalyst, yeld)
## PARAMETERS ARE THE DESIGN MATRIX AND THE RESPONSE VARIABLE
## MAIN FACTOR ANALYSIS
## Main factor = catalyst
sk1 <- SK(x=dm, y=yeld, model='yeld ~ labs * catalyst * pressure',
which='catalyst')
summary(sk1)
plot(sk1, title='catalyst')
## Main factor = pressure
sk2 <- SK(x=dm, y=yeld, model='yeld ~ labs * catalyst * pressure',
which='pressure')
summary(sk2)
plot(sk2, title='pressure')
## NESTED ANALYSIS
## Nested catalyst/pressure=1 -> SK.nest
nsk1 <- SK.nest(x=dm, y=yeld, model='yeld ~ labs * catalyst * pressure',
which='catalyst:pressure', fl2=1)
summary(nsk1)
plot(nsk1, title='catalyst/pressure=1')
## Nested catalyst/pressure=2 -> SK.nest
nsk2 <- SK.nest(x=dm, y=yeld, model='yeld ~ labs * catalyst * pressure',
which='catalyst:pressure', fl2=2)
summary(nsk2)
plot(nsk2, title = 'catalyst/pressure=2')
## Nested catalyst/pressure=2/labs=1 -> SK.nest
nsk3 <- SK.nest(x=dm, y=yeld, model='yeld ~ labs * catalyst * pressure',
which='labs:catalyst:pressure', fl2=2, fl3=1)
summary(nsk3)
plot(nsk3, title='catalyst/pressure=2/labs=1')