我正在创建一个散点图矩阵GGally::ggpairs
.我正在使用自定义函数(下面称为my_fn
)来创建左下角的非对角线子图.在调用该自定义函数的过程中,有关于计算的每个子图的信息,以及我想要存储以供日后使用的子图.
在下面的示例中,每个h@cID
都是具有100个值的int []结构.总共,它在my_fn中创建了10次(对于10个左下角的非对角线子图中的每一个,都会创建一次).我试图将所有这10个h@cID
结构存储到listCID
列表对象中.
我没有成功使用这种方法,并且我尝试了一些其他变体(例如尝试将listCID
输入参数作为输入参数my_fn
,或者最终尝试返回它).
我可以h@cID
有效地存储十个结构,my_fn
以便以后使用吗?我觉得有几个我不完全熟悉的语法问题可以解释为什么我被卡住了,同样如果我没有使用适当的术语,我也很乐意改变这个问题的标题.谢谢!
library(hexbin) library(GGally) library(ggplot2) set.seed(1) bindata <- data.frame( ID = paste0("ID", 1:100), A = rnorm(100), B = rnorm(100), C = rnorm(100), D = rnorm(100), E = rnorm(100)) bindata$ID <- as.character(bindata$ID ) maxVal <- max(abs(bindata[ ,2:6])) maxRange <- c(-1 * maxVal, maxVal) listCID <- c() my_fn <- function(data, mapping, ...){ x <- data[ ,c(as.character(mapping$x))] y <- data[ ,c(as.character(mapping$y))] h <- hexbin(x=x, y=y, xbins=5, shape=1, IDs=TRUE, xbnds=maxRange, ybnds=maxRange) hexdf <- data.frame(hcell2xy(h), hexID=h@cell, counts=h@count) listCID <- c(listCID, h@cID) print(listCID) p <- ggplot(hexdf, aes(x=x, y=y, fill=counts, hexID=hexID)) + geom_hex(stat="identity") p } p <- ggpairs(bindata[ ,2:6], lower=list(continuous=my_fn)) p
如果我正确理解你的问题,这很容易,虽然不是很优雅,但是使用<<-
操作员来实现.
有了它,您可以在函数范围内分配类似全局变量的内容.
设置listCID <- NULL
执行函数之前和listCID <<-c(listCID,h@cID)
在函数内部.
listCID = NULL my_fn <- function(data, mapping, ...){ x = data[,c(as.character(mapping$x))] y = data[,c(as.character(mapping$y))] h <- hexbin(x=x, y=y, xbins=5, shape=1, IDs=TRUE, xbnds=maxRange, ybnds=maxRange) hexdf <- data.frame (hcell2xy (h), hexID = h@cell, counts = h@count) if(exists("listCID")) listCID <<-c(listCID,h@cID) print(listCID) p <- ggplot(hexdf, aes(x=x, y=y, fill = counts, hexID=hexID)) + geom_hex(stat="identity") p }
有关范围的更多信息,请参阅Hadleys优秀的高级R:http://adv-r.had.co.nz/Environments.html