现代统计图形 - 科学网—博客

104 第五章 图库
1 > panel.hist = function(x, ...) {
2 + usr = par("usr")
3 + on.exit(par(usr))
4 + par(usr = c(usr[1:2], 0, 1.5))
5 + h = hist(x, plot = FALSE)
6 + nB = length(breaks idx = as.integer(iris[["Species"]])
11 > pairs(iris[1:4], upper.panel = function(x, y, ...) points(x,
12 + y, pch = c(17, 16, 6)[idx], col = idx), pch = 20,
13 + oma = c(2, 2, 2, 2), lower.panel = panel.smooth,
14 + diag.panel = panel.hist)
2.0 3.0 4.0
0.5 1.5 2.5
Sepal.Length
4.5 5.5 6.5 7.5
2.0 2.5 3.0 3.5 4.0 0.5 1.0 1.5 2.0 2.5
Sepal.Width
Petal.Length
1 2 3 4 5 6 7
Petal.Width
图 5.23: 鸢尾花数据的散点图矩阵:上三角区域为不同样式的点,对应着不
同种类的鸢尾花,对角线的直方图展示了花瓣花萼长宽的一维分布,下三
角区域用平滑曲线显示了变量之间的关系。
4.5 5.5 6.5 7.5
1 2 3 4 5 6 7