r - 2因素直方图分析

Question

我找了很久这个问题的答案，但一直没能找到答案。

问题是:我有一个 data frame具有以下变量:流速 1 (CH_SONAR)、流速 2 (CH_SONAR_2T)、密度 (CH_DENSITY) 以及两个流速之间的百分比差异 (per_diff)。我为流量 1 创建了一个 5 级因子，为密度创建了另一个 5 级因子。

f.factor <- cut(p.pipeline$CH_SONAR_2T, 5, labels = c('Very Low','Low', 'Medium', 'High', 'Very High'))

d.factor <- cut(p.pipeline$CH_DENSITY, 5, labels = c('Water', 'Very Sparce', 'Sparce', 'Dense', 'Very Dense'))

我使用 ggplot2 绘制了这两个图，并使用每个因子作为填充变量:

qplot(per_diff, data = p.pipeline, geom = "histogram", binwidth = 1, xlim = c(-5, 15), fill = f.factor)

qplot(per_diff, data = p.pipeline, geom = "histogram", binwidth = 1, xlim = c(-5, 15), fill = d.factor)

现在我想用 ggplot 创建一个直方图，让我看到流量和密度之间的关系(水和非常低、非常稀疏和低、稀疏和低等所有 25 种可能的组合)。我尝试过创建新因子、将 d.factor 和 f.factor 绑定(bind)到数据框、将两个因子绑定(bind)在一起等等，但没有结果，你们知道如何解决这个问题吗？

我尝试过包含我生成的直方图，但我认为我没有足够的声誉来做到这一点。

感谢您的帮助!

Answer 1

您可以使用fill=interaction(f.factor, d.factor)。图例中未出现的组合(例如“Low.Very Sparse”)表示不存在属于这两个类别的观察结果。

如果您希望相邻级别的颜色更加突出，您可以做的一件事是使用彩虹生成颜色，然后将所有其他颜色与色轮上相反的颜色交换。

col <- rainbow(length(levels(interaction(f.factor, d.factor))), v=.75, s=.5)
col.index <- ifelse(seq(col) %% 2, 
                    seq(col), 
                    (seq(ceiling(length(col)/2), length.out=length(col)) %% length(col)) + 1)
mixed <- col[col.index]
qplot(per_diff, data = p.pipeline, 
      geom = "histogram", binwidth = 1, xlim = c(-5, 15), 
      fill = interaction(f.factor, d.factor)) + scale_fill_manual(values=mixed)