r - 图中每个点的高度值

Question

我在标题为：s1m 的数据框中有蛋白质-蛋白质相互作用的数据。每个 DB 和 AD 对都进行交互，我也可以绘制它：

> head(s1m)
     DB_num AD_num
[1,]      2   8153
[2,]      7   3553
[3,]      8   4812
[4,]     13   7838
[5,]     24   3315
[6,]     24   6012

数据图如下所示：

然后我使用我在这个网站上找到的代码来绘制填充轮廓线：

## compute 2D kernel density, see MASS book, pp. 130-131
require(MASS)
z <- kde2d(s1m[,1], s1m[,2], n=50)
plot(s1m, xlab="X label", ylab="Y label", pch=19, cex=.4)
filled.contour(z, drawlabels=FALSE, add=TRUE)

它给了我生成的图像（减去涂鸦）：

我的问题：我需要在原始s1m数据框中的每一行数据上标注一个对应于等高线图高度的数字（因此我在上图中涂鸦）。我认为列表 z 具有我正在寻找的值，但我不确定。

最后，我希望我的数据看起来像这样，这样我就可以分组研究蛋白质相互作用：

         DB_num AD_num   height
    [1,]      2   8153        1
    [2,]      7   3553        1
    [3,]      8   4812        3
    [4,]     13   7838        6
    [5,]     24   3315        2
    [6,]     24   6012        etc.

Answer 1

如果您想要实际高度而不是每个分配给的箱，这是一种选择

## dummy data
DF <- data.frame(DB_num = rnorm(10000), AD_num = rnorm(10000))

require("MASS")

kde <- kde2d(DF[,1], DF[,2], n = 50)

请注意kde2d作为组件的返回，z它是一个具有（在这种情况下）50 行和列的矩阵，其中行对应于x数据，列对应于y数据。由于矩阵只是一个向量，并且数据按列填充，我们可以利用这一点并将x和y值n每次堆叠（n = 50此处），然后展开kde$z

dd <- dim(kde$z)
res <- data.frame(DB_num = rep(kde$x, times = dd[1]),
                  AD_num = rep(kde$y, times = dd[2]),
                  height = as.numeric(kde$z))

这产生

> head(res)
        DB_num      AD_num                                  height
1 -3.582508378 -3.79074271 0.0000000000000000000000000006907447484
2 -3.429230262 -3.63682706 0.0000000000000000000000002951259863229
3 -3.275952146 -3.48291141 0.0000000000000000000000558203373144190
4 -3.122674029 -3.32899576 0.0000000000000000000055565720524140235
5 -2.969395913 -3.17508011 0.0000000000000000014967010810961022503
6 -2.816117797 -3.02116446 0.0000000000000008159370528768207499471

要获得垃圾箱，您需要遵循所做的filled.contour事情，即通过

nlevels <- 20 ## default
brks <- pretty(range(res$height), nlevels)

> brks
 [1] 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14
[16] 0.15 0.16

然后使用cut将每个分配height到一个 bin 的基础上brks，例如

res <- transform(res, bin = as.numeric(cut(height, brks)))

这使

> head(res)
        DB_num      AD_num                                  height bin
1 -3.582508378 -3.79074271 0.0000000000000000000000000006907447484   1
2 -3.429230262 -3.63682706 0.0000000000000000000000002951259863229   1
3 -3.275952146 -3.48291141 0.0000000000000000000000558203373144190   1
4 -3.122674029 -3.32899576 0.0000000000000000000055565720524140235   1
5 -2.969395913 -3.17508011 0.0000000000000000014967010810961022503   1
6 -2.816117797 -3.02116446 0.0000000000000008159370528768207499471   1

您可能需要检查详细信息?cut以确定垃圾箱边界上的行为，但这应该让您足够接近。