显示双变量数据的散点图矩阵

Question

我正在尝试做应用介绍中的一个例子 R 书的多变量分析。我不知道如何得到这个答案。您可以使用获取数据 library(HSAUR2) data("pottery")

例如。 2.5 构建化学成分的散点图矩阵第 1 章 (Table 1.3) 中给出的罗马-不列颠陶器，识别每个以其窑炉编号为单位并显示估计的双变量密度每个面板。

我知道这段代码是我需要的代码示例，但我无法更改所有内容以适应陶器数据：

library("KernSmooth")
CYGOB1d <- bkde2D(CYGOB1, bandwidth = sapply(CYGOB1, dpik))
plot(CYGOB1, xlab = "log surface temperature", ylab = "log light intensity")
contour(x = CYGOB1d$x1, y = CYGOB1d$x2, z = CYGOB1d$fhat, add = TRUE)

Answer 1

您的代码没有生成散点图矩阵。它产生一个面板。这是一种使您包含的代码适应 pottery 数据的方法：

potteryd <- bkde2D(pottery[, 1:2], sapply(pottery[, 1:2], dpik))
plot(pottery[, 1:2], pch=as.numeric(pottery$kiln))
contour(x = potteryd$x1, y = potteryd$x2, z = potteryd$fhat, add = TRUE)
legend("topleft", levels(pottery$kiln), pch=as.numeric(levels(pottery$kiln)), title="Pottery Kiln")

显示双变量数据的散点图矩阵

Scatterplot Matrix showing bivariate data

r

graph

computer-science

scatter-plot