interaction2wt
From HH v3.1-19
by Richard Heiberger
Plot all main effects and twoway interactions in a multifactor design
The main diagonal
displays boxplots for the main effects of each factor. The
off-diagonals show the interaction plots for each pair of factors.
The i,j
panel shows the same factors as the j,i
but with
the trace- and x-factor roles interchanged.
Usage
interaction2wt(x, ...)
## S3 method for class 'formula':
interaction2wt(x, data = sys.parent(), responselab, ...)
## S3 method for class 'default':
interaction2wt(x,
response.var,
responselab = deparse(substitute(response.var)),
responselab.expression = responselab,
relation = list(x = "same", y = "same"),
x.relation = relation$x,
y.relation = relation$y,
digits = 3,
x.between=1,
y.between=1,
between,
cex = 0.75,
rot=c(0,0),
panel.input = panel.interaction2wt,
strip.input =
if (label.as.interaction.formula) strip.default
else strip.interaction2wt,
par.strip.text.input = trellis.par.get()$add.text,
scales.additional,
main.in =
paste(responselab,
": ", c("main", "simple")[1+simple],
" effects and 2-way interactions",
sep=""),
xlab = "",
ylab = "",
simple=FALSE,
box.ratio=if (simple) .32 else 1,
label.as.interaction.formula=TRUE,
...,
main.cex,
key.cex.title=.8,
key.cex.text=.7,
factor.expressions=names.x,
simple.pch=NULL
)
Arguments
- x
- The object on which method dispatch is carried out.
For the
"formula"
method, a formula describing the response variable and factors. The formula is generally of the formy ~ g1 + g2 + ...
. There may be one - data
- For the
formula
method, a data frame containing values for any variables in the formula. In the R version, if not found indata
, or ifdata
is unspecified, the variables are looked for in the env - responselab
- Character name of response variable, defaults to
the name of the response variable in the
formula
. - responselab.expression
plotmath
or character name of response variable, defaults toresponselab
.- ...
- additional arguments, primarily trellis arguments.
- response.var
- For the
"default"
method, the response variable. This is usually constructed byformula
method from the input data and the input formula. - simple
- logical.
TRUE
if simple effects are to be displayed. Argumentssimple.offset
,simple.scale
, andcol.by.row
may also be needed. Seepan
- box.ratio
- #ifndef S-Plus
xyplot
. #endif #ifdef S-Plusbwplot
. #endif - relation
trellis
argument.- x.relation
x
value ofrelation
argument.- y.relation
y
value ofrelation
argument.- digits
- doesn't do anything at the moment
- x.between
x
value ofbetween
argument.- y.between
y
value ofbetween
argument.- between
- trellis/lattice
between
argument. If used,between
has precedence over both thex.between
andy.between
arguments. - cex
- S-Plus: changes the size of the median dot in the boxplots. R: doesn't do anything.
- panel.input
- panel function. Default is
panel.interaction2wt
. - label.as.interaction.formula
- logical. If
TRUE
, each panel has a single strip label of the formy ~ a | b
. IfFALSE
, each panel has a pair of strip labels, one for thetrace
factor and one for thex
factor. - strip.input
- strip function. Default depends on the value of
label.as.interaction.formula
. - par.strip.text.input
par.strip.text
argument.- scales.additional
- additional arguments to
scales
argument ofinteraction.positioned
. - main.in
- Text of main title.
- xlab
- No effect.
- ylab
- No effect.
- main.cex
cex
for main title.- key.cex.title
cex
key title. Defaults tocex
forxlab
.- key.cex.text
cex
group names in key. Defaults tocex
foraxis.text
.- factor.expressions
- Expressions for titles of keys and
xlab
for each column. Defaults to the names of the factors in the input formula. - rot
- Rotation of x tick labels and y tick labels. Only 0 and 90 will look good.
- simple.pch
- Named list containing plotting characters for each
level of one or more of the factors.
simple.pch
is used only whensimple==TRUE
. If the argumentsimple.pch
is missing, then the integers for the lev
Value
"trellis"
object containing the plot.
References
Heiberger, Richard M. and Holland, Burt (anticipated 2015). Statistical Analysis and Data Display: An Intermediate Course with Examples in R, Second Edition. Springer Texts in Statistics. Springer. ISBN 978-1-4939-2121-8. Heiberger, Richard M. and Holland, Burt (2004b). Statistical Analysis and Data Display: An Intermediate Course with Examples in S-Plus, R, and SAS. Springer Texts in Statistics. Springer. ISBN 0-387-40270-5.
See Also
Examples
data(vulcan)
interaction2wt(wear ~ filler + pretreat + raw, data=vulcan,
par.strip.text=list(cex=.7))
interaction2wt(wear ~ filler + pretreat + raw, data=vulcan)
interaction2wt(wear ~ filler + raw, data=vulcan,
simple=TRUE)
interaction2wt(wear ~ filler + raw, data=vulcan,
simple=TRUE, simple.scale=c(filler=.15, raw=.2),
xlim=c(.3, 5.6))
ToothGrowth$dose <- positioned(ToothGrowth$dose) ## modify local copy
anova(aov(len ~ supp*dose, data=ToothGrowth))
interaction2wt(len ~ supp + dose, data=ToothGrowth)
esoph$rate=with(esoph, ncases/ncontrols) ## modify local copy
position(esoph$alcgp) <- 2:5
position(esoph$tobgp) <- 2:5
interaction2wt(rate ~ agegp + alcgp + tobgp, esoph, rot=c(90,0),
par.strip.text=list(cex=.8))
interaction2wt(rate ~ agegp + alcgp + tobgp, esoph, rot=c(90,0),
par.strip.text=list(cex=.8),
factor.expressions=c(
agegp=expression(Age~~(years)),
alcgp=expression(Alcohol~
bgroup("(",scriptstyle(frac(gm, day)),")")),
tobgp=expression(Tobacco~
bgroup("(",scriptstyle(frac(gm, day)),")"))),
par.settings=list(
par.xlab.text=list(cex=.8),
par.ylab.text=list(cex=.8)),
responselab.expression="Cancer
Rate",
main=list(
"Esophogeal Cancer Rate ~ Alcohol Consumption + Tobacco Consumption",
cex=1.2))
esoph.aov <- aov(rate ~ agegp + alcgp + tobgp, data=esoph)
anova(esoph.aov)
Community examples
Looks like there are no examples yet.