6

我想使用函数(R 包)中的stat_binhex()统计数据。例如,我想在这个情节中使用而不是. 那可能吗?ggplot2ggpairs()GGallystat_binhex()geom_point()

在此处输入图像描述

谢谢你的帮助!

4

2 回答 2

4
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

截屏

于 2014-01-01T20:00:12.327 回答
0

不敢相信这仍然不是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)

在此处输入图像描述 在此处输入图像描述 在此处输入图像描述

于 2021-10-03T19:32:56.780 回答