我想使用函数(R 包)中的stat_binhex()
统计数据。例如,我想在这个情节中使用而不是. 那可能吗?ggplot2
ggpairs()
GGally
stat_binhex()
geom_point()
谢谢你的帮助!
set.seed(1)
library(GGally)
library(hexbin)
df <- as.data.frame(matrix(rnorm(20*3), ncol=3))
p <- ggpairs(df, lower="blank")
seq <- 1:ncol(df)
for (x in seq)
for (y in seq)
if (y>x)
p <- putPlot(p, ggplot(df, aes_string(x=names(df)[x],y=names(df)[y])) + stat_binhex(bins=4), y,x)
p
不敢相信这仍然不是GGally
' ggpair()
s 定制的一部分。
基于lukeA 的回答,让我们将其转换为函数:
设置要求和定义功能:
require(ggplot2)[enter image description here][1]
require(GGally)
require(hexbin)
ggpairs_hex <- function(df, hexbins = 10) {
# REF: https://stackoverflow.com/questions/20872133/using-stat-binhex-with-ggpairs
p <- ggpairs(df, lower="blank")
seq <- 1:ncol(df)
for (x in seq)
for (y in seq)
if (y>x)
p <- putPlot(p, ggplot(df, aes_string(x=names(df)[x],y=names(df)[y])) + stat_binhex(bins = hexbins), y,x)
return(p)
}
生成一些数据:
require(MASS)
# data generation from:
# https://predictivehacks.com/how-to-generate-correlated-data-in-r/
set.seed(312)
# create the variance covariance matrix
sigma<-rbind(c(1,-0.8,-0.7), c(-0.8,1, 0.9), c(-0.7,0.9,1))
# create the mean vector
mu<-c(10, 5, 2)
# generate the multivariate normal distribution
df<-as.data.frame(MASS::mvrnorm(n=10000, mu=mu, Sigma=sigma))
测试功能:
ggpairs_hex(df, hexbins = 5)
ggpairs_hex(df, hexbins = 10)
ggpairs_hex(df, hexbins = 20)