2

我需要使用 R 绘制二维密度图filled.contour。我有两个数据集并将它们绘制成两个fill.contour. 我在这里没有 10 名声望，也无法在这里发布我的数据。我在这里发布我的代码，希望它可以帮助找出问题。

library(MASS)
density <- kde2d(multi_ligand[,21], multi_ligand[,7])
filled.contour(density,
    color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')),
    xlab=colnames(single_ligand[21]),
    ylab=colnames(single_ligand[7])
)

density1 <- kde2d(single_ligand[,21], single_ligand[,7])
filled.contour(density1,
    color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')),
    xlab=colnames(single_ligand[21]),
    ylab=colnames(single_ligand[7])
)

问题是我发现在这两个图中，相同密度的颜色不一样。例如，在第一个图中，密度 0.06 是黄色，但在第二个图中，密度 0.06 是蓝色。我在这两个图中使用相同的色标。为了使这两个图具有可比性，我想在这两个图中以相同的密度使用相同的颜色。

任何人都可以告诉我，我应该如何更改我的设置以使其正确？

4

1 回答 1

4

默认情况下filled.contour，将调整颜色块以均匀覆盖每个数据集的 值范围z，或者在这种情况下density为 值。如果您希望在两个图上使用完全相同的级别，则需要自己指定它们。这里有一些代码将指定覆盖两个数据集范围的级别。

#sample data
set.seed(15)
ax<-rnorm(50) #like multi_ligand[,21]
ay<-rnorm(50) #like multi_ligand[,7]
bx<-rnorm(75,2, .5) #like single_ligand[,21]
by<-rnorm(75,2, .5) #like single_ligand[,7]

#calculate both densities
density <- kde2d(ax, ay)
density1 <- kde2d(bx, by)

#make levels that cover both ranges of z values
lvls <- pretty(range(density$z, density1$z),20)

#draw both plots using the same levels
filled.contour(density,
    color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')),
    levels=lvls
)

filled.contour(density1,
    color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')),
    levels=lvls
)

产生这两个地块

样本输出

于 2014-05-20T21:05:41.080 回答