R plot.design 绘制设计或模型的单变量效应

说明

绘制一个或多个 factor 的单变量效应，通常用于 aov() 分析的设计实验。

用法

plot.design(x, y = NULL, fun = mean, data = NULL, ...,
            ylim = NULL, xlab = "Factors", ylab = NULL,
            main = NULL, ask = NULL, xaxt = par("xaxt"),
            axes = TRUE, xtick = FALSE)

参数

x	包含设计因子和可选响应的 DataFrame ，或者 formula 或 terms 对象。
y	响应(如果 x 中未给出)。
fun	应用于每个子集的函数(或函数名称)。对于数字(向量)输入，它必须返回一个数字。
data	数据帧包含 x(当 formula-like 时)引用的变量。
...	graphical parameters 例如col ，参见par 。
ylim	y 值的范围，如plot.default 中。
xlab	x 轴标签，请参见 title 。
ylab	y 轴标签，默认值为‘smart’。
main	主标题，请参见title。
ask	逻辑指示在新页面开始之前是否应询问用户 - 在多个 y 的情况下。
xaxt	给出 x 轴类型的字符。
axes	逻辑指示是否应绘制轴。
xtick	逻辑指示是否应在 x 轴上绘制刻度线(每个因子一个)。

细节

所提供的函数将针对设计中每个因子的每个级别调用一次，并且绘图将显示这些汇总值。特定因子的水平沿垂直线显示，响应的 fun() 的总体值绘制为水平线。

注意

为了使其与 S 版本紧密兼容，我们付出了巨大的努力。但是， col (和 fg )规范具有不同的效果。

在 S 中，这是 design 对象的 plot 通用函数的方法。

例子

require(stats)
plot.design(warpbreaks)  # automatic for data frame with one numeric var.

Form <- breaks ~ wool + tension
summary(fm1 <- aov(Form, data = warpbreaks))
plot.design(       Form, data = warpbreaks, col = 2)  # same as above

## More than one y :
utils::str(esoph)
plot.design(esoph) ## two plots; if interactive you are "ask"ed

## or rather, compare mean and median:
op <- par(mfcol = 1:2)
plot.design(ncases/ncontrols ~ ., data = esoph, ylim = c(0, 0.8))
plot.design(ncases/ncontrols ~ ., data = esoph, ylim = c(0, 0.8),
            fun = median)
par(op)

作者

Roberto Frisullo and Martin Maechler

参考

Chambers, J. M. and Hastie, T. J. eds (1992) Statistical Models in S. Chapman & Hall, London, the white book, pp. 546-7 (and 163-4).

Freeny, A. E. and Landwehr, J. M. (1990) Displays for data from large designed experiments; Computer Science and Statistics: Proc.\ 22nd Symp\. Interface, 117-126, Springer Verlag.

也可以看看

interaction.plot 用于设计实验的“标准图形”。