同じグラフ内に整列されたヒストグラムをプロットする

Question

整列されたヒストグラムを同じグラフのマトリックスの行からプロットする必要があります。次の例では、積み重なった5つのヒストグラムをプロットする必要があります。ただし、histコマンドは行列全体のヒストグラムをプロットします。私が使用できる回避策はありますか？

ありがとうございます。

アルトゥーロ

x <- 10; y <- 10; g <- 5 
dat <- matrix(rnorm((x + y) * g), ncol = x + y) 
hist (dat) 

Answer 1

これは私が同じチャートで2つのヒストグラムを生成するために使用していたアプローチです。

列xが連続変数であり、列yが係数変数であるデータセットを持つことを検討してください。各レベルは1つのヒストグラムを生成します。下記のggplot2を使用したコードの例をご覧ください：

OverlayedHist <- function(mData  , 
          featureVar , 
          grouper  , 
          mbinwidth , 
          mTitle  , 
          mxlab  , 
          mylab  , 
          mlegendTitle 
){ 

    # function name: OverlayedHist 
    #  purpose: To produce overlayed histograms against a grouping variable 
    #   Input: 
    #   mData: dataset object in data.table format 
    # featureVar: name of continuous variable to produce histogram 
    #  grouper: the grouping variable to produce the histogram 
    #  mbinwidth: binwidth (see ggplot2 parameters) 
    #  mTitle: Character to define title 
    #   mxlab: Character to define xlab name 
    #   mylab: Character to define ylab name 
    # mlegendTitle: Character to define legend title 

    library(data.table) 
    library(ggplot2) 
    library(plotly) 

    p <- ggplot(allDat, aes(eval(parse(text = featureVar)), fill = eval(parse(text = grouper)))) + 
    geom_histogram(alpha = 0.7, position = 'identity', binwidth = mbinwidth) + 
    scale_fill_manual(values=c("#377EB8","#E41A1C")) + 
    ggtitle(mTitle) + 
    xlab(mxlab) + ylab(mylab) + 
    guides(fill=guide_legend(title=mlegendTitle)) + theme(plot.title = element_text(size=10)) 


    return(ggplotly(p)) 

} 

私はこれが役立つことを願っています。

乾杯！ K.

Answer 2

ggplot簡単に（長い形式のデータが必要です）。

library(tidyr); library(dplyr); library(ggplot2) 

df <- dat %>% t() %>% as.data.frame() %>% gather(row) # chage data into a long format 
Breaks <- hist(dat, plot=F)$breaks      # get (or decide) breaks 

ggplot(df, aes(x = value, fill = row)) + geom_histogram(position = "stack", breaks = Breaks) 

[EDITED]
はこれが何をしたいですか？

## original 
ggplot(df, aes(x = value)) + 
    geom_histogram(breaks = Breaks) + 
    facet_wrap(~ row)   # make histogram par group. 

    ## modified 
ggplot(df, aes(x = value, fill = row)) +     # change fill colour 
    geom_histogram(breaks = Breaks) + 
    facet_wrap(~ row, ncol = 1) +       # bring graphs into line 
    # facet_wrap(~ row, ncol = 1, scales = "free_y") # if you don't want fixed scale 
    theme(strip.background = element_blank(), strip.text.x = element_blank()) # delete labels 

[EDITED2：base_plotアプローチ]
基本プロット時間を節約できます。

## example data 
x <- 1500; y <- 1500; g <- 30 
set.seed(1); dat <- matrix(rnorm((x + y) * g), ncol = x + y) 

## decide breaks 
Breaks <- seq(-4.5, 4.5, 0.5) 

## change par() to draw multiple graphs 
par.old <- par(mar = c(0.1, 4.0, 0.1, 0.5), oma = c(4, 0, 0, 0), mfrow = c(nrow(dat), 1)) 

for(i in 1:nrow(dat)) { 
    hist(dat[,i], breaks = Breaks, xaxt = "n", xlab = "", main = "") 
# grid(NULL) 
    } 
par(new=T) 
hist(dat[,nrow(dat)], breaks = Breaks, main = "", yaxt = "n", 
    xlab = "x", ylab = "", border = NA)   # add x-axis 

par(par.old)