绘制常规2D图形，但添加第三维作为热图

Question

找不到我在哪里问，也许错误的关键词。基本上，我在矩阵中有3个维度：

> head(info)
        [,1]   [,2]  [,3]
[1,] 8.59645 251944 22.89
[2,] 6.95160 141559 21.35
[3,] 7.43870 131532 22.99
[4,] 8.64467 126688 22.72
[5,] 8.77482 123120 22.17
[6,] 7.22364 122268 24.46

我正在密封信息[，3]与信息[，2]

plot(info[,3], info[,2], type="p", pch=20)

我希望使用基于info [，1]的热图来为点着色。

我可以做这样的事情：

plot(info[which(info[,1] <= 2),3], info[which(info[,1] <= 2),2], type="p", pch=20, col="black")
lines(info[which(info[,1] >= 2),3], info[which(info[,1] >= 2),2], type="p", pch=20, col="red")

但我相信热图会更好看。

有什么想法吗？谢谢， 阿德里安

解： 谢谢大家的精彩建议！这是有用的：

qplot(info[,3], info[,2], colour=info[,1]) + scale_colour_gradient(limits=c(0, 10), low="green", high="red")

Answer 1

使用ggplot2，您可以按第三个变量着色

## Some sample data
set.seed(0)
x <- rnorm(1000, rep(c(20, 60), each=500), 8)
y <- c(rexp(500, 1/5e4)*1/(abs(x[1:500]-mean(x[1:500]))+runif(1)),
       rexp(500, 1/5e3)*1/(abs(x[501:1000]-mean(x[501:1000]))+runif(1)))
z <- c(sort(runif(1000)))
info <- matrix(c(z,y,x), ncol=3)

## Using ggplot
ggplot(as.data.frame(info), aes(V3, V2, col=V1)) +
  geom_point(alpha=0.5) +
  scale_color_gradient(low="red", high="yellow")

如果您想制作热图，可以使用akima包在您的点之间进行插值，并且确实如此，

library(akima)
dens <- interp(x, y, z,
               xo=seq(min(x), max(x), length=100),
               yo=seq(min(y), max(y), length=100),
               duplicate="median")
filled.contour(dens, xlab="x", ylab="y", main="Colored by z",
               color.palette = heat.colors)