我有一个包含 100 个观察值的双变量数据集。我使用了六边形分箱,最终得到了 26 个六边形分箱。为了保存 26 个六边形箱中每个箱中 100 个观测值的行,我使用了base::attr
R 中的函数。在下面的代码中,这是在以下位置完成的:
attr(hexdf, "cID") <- h@cID
我正在尝试创建一个交互式 RPlotly
六边形分箱的对象,这样如果用户单击给定的六边形分箱,他们将获得分组到该箱中的 100 个观察值的行。我已经完成了这个目标的一部分。我的 MWE 如下:
library(plotly)
library(data.table)
library(GGally)
library(hexbin)
library(htmlwidgets)
set.seed(1)
bindata <- data.frame(ID = paste0("ID",1:100), A=rnorm(100), B=rnorm(100))
bindata$ID <- as.character(bindata$ID)
x = bindata[,c("A")]
y = bindata[,c("B")]
h <- hexbin(x=x, y=y, xbins=5, shape=1, IDs=TRUE)
hexdf <- data.frame (hcell2xy (h), hexID = h@cell, counts = h@count)
attr(hexdf, "cID") <- h@cID
pS <- ggplot(hexdf, aes(x=x, y=y, fill = counts, hexID=hexID)) + geom_hex(stat="identity")
ggPS <- ggplotly(pS)
myLength <- length(ggPS[["x"]][["data"]])
for (i in 1:myLength){
item =ggPS[["x"]][["data"]][[i]]$text[1]
if (!is.null(item))
if (!startsWith(item, "co")){
ggPS[["x"]][["data"]][[i]]$hoverinfo <- "none"
}
}
ggPS %>% onRender("
function(el, x, data) {
//console.log(el)
//console.log(x)
//console.log(data)
myGraph = document.getElementById(el.id);
el.on('plotly_click', function(e) {
cN = e.points[0].curveNumber
split1 = (x.data[cN].text).split(' ')
hexID = (x.data[cN].text).split(' ')[2]
counts = split1[1].split('<')[0]
console.log(hexID)
console.log(counts)
})}
", data = pS$data)
当我运行此代码并在 Web 浏览器中打开它时,我获得如下所示的交互式图(绿色框不在图中;出于解释目的而将其叠加):
如果我点击绿色框内的六边形,正确的hexID
40 和counts
3 个被打印到控制台。此时,我想获取放入该六边形箱中的 3 行原始数据帧。
我知道如何在 R 之外执行此操作onRender()
的功能htmlwidgets
包通过使用base::attr
功能。例如,我可以执行以下操作:
hexID=40
obsns <- which(attr(pS$data, "cID")==hexID)
dat <- bindata[obsns,]
并收到以下正确的 3 个数据点,这些数据点已放入我单击的垃圾箱中:
ID A B
47 ID47 0.3645820 2.087167
66 ID66 0.1887923 2.206102
71 ID71 0.4755095 2.307978
我正在处理比这个 MWE 中更大的数据集。因此,我打算使用base:attr
功能是防止更大的数据框四处浮动。但是,我不确定如何翻译该功能base::attr
函数,以便我可以访问出现在单击的六边形箱中的适当数据点行onRender()
JavaScript 代码。我确实包括了pS$data
对象进入onRender()
JavaScript 代码,但我仍然卡住了。
任何建议将不胜感激!