panel.interaction2wt: Plot all main effects and twoway interactions in a...

In HH: Statistical Analysis and Data Display: Heiberger and Holland

Plot all main effects and twoway interactions in a multifactor design

Description

This is the panel function for interaction2wt. 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

panel.interaction2wt(x, y, subscripts,
                     responselab, trace.values,
                     factor.levels, factor.position,
                     fun = mean,
                     se,
                     type="l",
                     ...,
                     box.ratio,
                     simple=FALSE,
                     simple.offset,
                     simple.scale,
                     simple.pch,
                     data.x,
                     col.by.row=TRUE,
                     col  =trellis.par.get("superpose.line")$col,
                     lty  =trellis.par.get("superpose.line")$lty,
                     lwd  =trellis.par.get("superpose.line")$lwd,
                     alpha=trellis.par.get("superpose.line")$alpha
)

strip.interaction2wt(which.given, which.panel, var.name,
                     factor.levels, shingle.intervals,
                     strip.names = c(TRUE, TRUE), style = 1, ...)

Arguments

panel.interaction2wt arguments:

x	levels of x-factor
y	Summary value of response variable at each level of x- and trace-factors.
subscripts	used to get the right set of response values for the summary statistics on the off-diagonals
responselab	Character name of response variable, defaults to the name of the response variable.
trace.values	levels of trace-factor
fun	Summary function, defaults to mean
se	standard errors to be passed to panel.intxplot. Missing, logical, or a numeric vector. If se is missing or FALSE, or if simple is FALSE, then standard errors are not plotted. If TRUE, the standard errors are calculated from the sufficient statistics for each group as the group's standard deviation divided by the square root of the group's observation count. If a numeric vector, it is evaluated in the environment of the sufficient statistics.

,

type	See panel.xyplot.

,

box.ratio

passed to panel.bwplot.intermediate.hh

,

...	extra arguments, primarily color, to be passed to panel.bwplot.intermediate.hh
factor.position	"position" attribute of factor.
simple	logical. If TRUE, then simple effects are to be displayed.
simple.offset, simple.scale	named list of offset and scale for the response and trace factors. See interaction.positioned for their use.
simple.pch	Named list containing plotting characters for each level of one or more of the factors. simple.pch is used only when simple==TRUE. If the argument simple.pch is missing, then the integers for the levels of the factors are used. The characters are used for the median of the box plots in the diagonal panels. They match the trace factor of the interaction panel in the same column of the display.
data.x	data.frame containing factors from the input data.frame
col.by.row	logical. If TRUE (the default), simple effects plots color the simple effects on the main diagonals in the same color as the trace levels in their row. If FALSE, then simple effects are colored to match the x levels in their column.
col, lty, lwd, alpha	Arguments to trellis.par.set(superpose.line=list()).

strip.interaction2wt arguments

which.given, which.panel, var.name, factor.levels, shingle.intervals

see documentation for strip.default

.

strip.names	Force strip.names=TRUE
style	force style=1

Author(s)

Richard M. Heiberger <[e-mail]>

References

Heiberger, Richard M. and Holland, Burt (2015). Statistical Analysis and Data Display: An Intermediate Course with Examples in R. Second Edition. Springer-Verlag, New York. https://link.springer.com/book/10.1007/978-1-4939-2122-5

See Also

interaction2wt, panel.bwplot.intermediate.hh

Examples

## Not run: 
tmp <- data.frame(y=rnorm(48),
                  A=factor(rep(1:2, 24)),
                  B=factor(rep(rep(1:3, each=2), 8)),
                  C=factor(rep(rep(1:4, each=6), 2)))
interaction2wt(y ~ A+B+C, data=tmp,
               key.in=list(x=-3), ## key.in is ignored by R
               xlim=c(.4, 4.5))
interaction2wt(y ~ B+C, data=tmp, key.in=list(x=-2), xlim=c(.4, 4.5))
position(tmp$B) <- c(1, 2.4, 3.8)
interaction2wt(y ~ B+C, data=tmp, key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ B+C, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ C+B, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ B+C, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               simple.pch=list(C=c(16,17,18,19)),
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ C+B, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               simple.pch=list(C=c(16,17,18,19)),
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ C+B, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               simple.pch=list(A=c(1:2), B=c(3:5), C=c(16,17,18,19)),
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ C+B, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               simple.pch=list(A=c(1:2)),
               key.in=list(x=-2), xlim=c(.4, 4.5))
interaction2wt(y ~ B+C, data=tmp, simple=TRUE,
               simple.scale=list(B=.18, C=.27),  box.ratio=.2,
               simple.pch=list(B=c(16,17,18)),
               key.in=list(x=-2), xlim=c(.4, 4.5),
               se=TRUE)

## End(Not run)

HH documentation built on May 29, 2024, 6:24 a.m.