gpt4 book ai didi

带有单元格成员数量的 R 热图

转载 作者:行者123 更新时间:2023-12-02 08:19:35 25 4
gpt4 key购买 nike

我有一个数据集,用于创建热图。问题之一是,某些单元格中的成员很少,并且可能存在无法被其他成员平均掉的异常值。

为此,我想在单元格中包含单元格中实际有多少个示例(在单元格中心绘制)。

下面是我的热图代码:

library(fields)
library(akima)

x1 <- round(runif(20) * 100,0)
y1 <- round(runif(20) * 100,0)
z1 <- round(runif(20) * 100,0)

s <- interp(x1,y1,z1,
xo = seq(0,100,20)
,yo = seq(0,100,20)
)

image.plot(s)

有什么建议吗?

最佳答案

计算完单元格的角点和中心后,您可以使用 findIntervaltable 来计算观察值。

library(fields)
library(akima)

x1 <- floor(runif(20) * 100)
y1 <- floor(runif(20) * 100)
z1 <- floor(runif(20) * 100)

# Corners of the cells, to count the observations
x0 <- seq(0,100,20)
y0 <- seq(0,100,20)

# Centers of the cells, for the interpolation
x00 <- x0[-length(x0)] + diff(x0) / 2
y00 <- y0[-length(y0)] + diff(y0) / 2

s <- interp(x1,y1,z1, xo=x00, yo=y00)
image.plot(x=x0, y=y0, z=s$z)

counts <- table(
findInterval(x1, x0),
findInterval(y1, y0)
)
# Plot the observations, to check that I have not confused rows and columns
points( x1, y1 )
# Number of observations
text(x=x00[row(counts)], y=y00[col(counts)], labels=counts)

image.plot with counts

关于带有单元格成员数量的 R 热图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/17947441/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com