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。