Rでバブルプロット

三つ変数組のデータを二次元平面上にバブルプロットする描き方をメモっておく。

n <- 50
x <- rnorm(n)
y <- rnorm(n)
z <- rnorm(n)
re.norm <- function (z) {
z <- abs(z)
z <- 10*z/max(z)
z
}
z <- re.norm(z)
op <- par(mar = c(3,2,4,2)+.1)
plot(x, y, cex = z, pch = 16, col = "blue",
xlab = "", ylab = "",
main = "")
points(x, y, cex = z)