R coplot 調節圖

說明

該函數產生下麵參考文獻中討論的調節圖的兩種變體。

用法

coplot(formula, data, given.values, panel = points, rows, columns,
       show.given = TRUE, col = par("fg"), pch = par("pch"),
       bar.bg = c(num = gray(0.8), fac = gray(0.95)),
       xlab = c(x.name, paste("Given :", a.name)),
       ylab = c(y.name, paste("Given :", b.name)),
       subscripts = FALSE,
       axlabels = function(f) abbreviate(levels(f)),
       number = 6, overlap = 0.5, xlim, ylim, ...)
co.intervals(x, number = 6, overlap = 0.5)

參數

formula	說明條件圖形式的公式。 y ~ x | a 形式的公式表示應根據變量 a 生成 y 與 x 的圖。 y ~ x| a * b 形式的公式表示應根據兩個變量 a 和 b 生成 y 與 x 的圖。 所有三個或四個變量可以是數字或因子。當 x 或 y 是因子時，結果幾乎就像應用了 as.numeric() 一樣，而對於因子 a 或 b ，條件(及其圖形，如果 show.given 為 true)會進行調整。
data	包含公式中任何變量值的 DataFrame 。默認情況下，使用調用 coplot 的環境。
given.values	一個值或兩個值的列表，用於確定如何對 a 和 b 進行調節。 當沒有 b (即僅以 a 為條件)時，通常這是一個具有兩列的矩陣，每行給出一個間隔，以進行條件化，但也可以是單個數字向量或一組因子水平(如果所調節的變量是一個因子)。在這種情況下(沒有 b )，co.intervals 的結果可以直接用作 given.values 參數。
panel	function(x, y, col, pch, ...) 給出要在顯示屏的每個麵板中執行的操作。默認為 points 。
rows	該圖的麵板布置在 rows by columns 數組中。 rows 給出數組中的行數。
columns	麵板布局數組中的列數。
show.given	邏輯(對於 2 個條件變量，長度可能為 2)：是否應該顯示相應條件變量的條件圖(默認 TRUE )。
col	用於繪製點的顏色向量。如果太短，這些值將被回收。
pch	繪製符號或字符的向量。如果太短，這些值將被回收。
bar.bg	一個命名向量，其組件 "num" 和 "fac" 分別為數字和因子條件變量提供(木瓦)條的背景顏色。
xlab	特點;用於 x 軸和第一個條件變量的標簽。如果隻給出一個標簽，則該標簽用於 x 軸，默認標簽用於條件變量。
ylab	特點;用於 y 軸和任何第二個條件變量的標簽。
subscripts	邏輯：如果為 true，則為麵板函數提供一個附加(第三個)參數 subscripts，給出傳遞給該麵板的數據的下標。
axlabels	當 x 或 y 為因子時，用於創建軸(刻度)標簽的函數。
number	整數; a 和 b 的條件間隔數，長度可能為 2。僅當相應的條件變量不是 factor 時才使用它。
overlap	數字 < 1；條件變量的重疊部分，x 和 y 方向的長度可能為 2。當overlap<0時，數據切片之間就會有間隙。
xlim	x 軸的範圍。
ylim	y 軸的範圍。
...	麵板函數的附加參數。
x	一個數值向量。

細節

在單個條件變量 a 的情況下，當 rows 和 columns 均未指定時，將使用 columns >= rows 選擇“接近正方形”布局。

在多個 rows 的情況下，麵板圖的順序是從底部到左側(通常對應於增加 a )。

麵板函數不應嘗試開始新繪圖，而應僅在給定坐標係內進行繪圖：因此 plot 和 boxplot 不是麵板函數。

參數 xlab 和 ylab 的渲染不受 par 參數 cex.lab 和 font.lab 控製，即使它們是由 mtext 而不是 title 繪製的。

值

co.intervals(., number, .) 返回 ( number \times 2) matrix ，例如 ci ，其中 ci[k,] 是第 k 間隔的 x 值中的 range。

例子

## Tonga Trench Earthquakes
coplot(lat ~ long | depth, data = quakes)
given.depth <- co.intervals(quakes$depth, number = 4, overlap = .1)
coplot(lat ~ long | depth, data = quakes, given.values = given.depth, rows = 1)

## Conditioning on 2 variables:
ll.dm <- lat ~ long | depth * mag
coplot(ll.dm, data = quakes)
coplot(ll.dm, data = quakes, number = c(4, 7), show.given = c(TRUE, FALSE))
coplot(ll.dm, data = quakes, number = c(3, 7),
       overlap = c(-.5, .1)) # negative overlap DROPS values

## given two factors
Index <- seq_len(nrow(warpbreaks)) # to get nicer default labels
coplot(breaks ~ Index | wool * tension, data = warpbreaks,
       show.given = 0:1)
coplot(breaks ~ Index | wool * tension, data = warpbreaks,
       col = "red", bg = "pink", pch = 21,
       bar.bg = c(fac = "light blue"))

## Example with empty panels:
with(data.frame(state.x77), {
coplot(Life.Exp ~ Income | Illiteracy * state.region, number = 3,
       panel = function(x, y, ...) panel.smooth(x, y, span = .8, ...))
## y ~ factor -- not really sensible, but 'show off':
coplot(Life.Exp ~ state.region | Income * state.division,
       panel = panel.smooth)
})

參考

Chambers, J. M. (1992) Data for models. Chapter 3 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.

Cleveland, W. S. (1993) Visualizing Data. New Jersey: Summit Press.

也可以看看

pairs、panel.smooth、points。