6

関数(Rパッケージ)でstat_binhex()統計を使用したいと思います。たとえば、このプロットでは の代わりにを使用したいと思います。それは可能ですか?ggplot2ggpairs()GGallystat_binhex()geom_point()

ここに画像の説明を入力

ご協力いただきありがとうございます!

4

2 に答える 2

0

これがまだのカスタマイズGGallyの一部ではないことが信じられません。ggpair()

lukeA の answerに基づいて、それを関数に変換しましょう。

要件を設定し、関数を定義します。

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 に答える