在任何地方都找不到我要问的内容,也许是错误的关键词。本质上,我在一个矩阵中有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


我正在绘制info [,3] vs info [,2]

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


r - 绘制常规2D图形,但将第3维添加为热图-LMLPHP

我想用基于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")


r - 绘制常规2D图形,但将第3维添加为热图-LMLPHP

但是我认为,热图看起来会更好。

有任何想法吗?谢谢,
阿德里安

解:
谢谢大家的精彩建议!这是有效的:

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


r - 绘制常规2D图形,但将第3维添加为热图-LMLPHP

最佳答案

使用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")


r - 绘制常规2D图形,但将第3维添加为热图-LMLPHP

如果要制作热图,可以使用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)


r - 绘制常规2D图形,但将第3维添加为热图-LMLPHP

10-05 21:10
查看更多