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 用於設計實驗的“標準圖形”。