用只有一个变量的ggplot2绘制R中的小面网格

Question

我有一个名为mouse.data的数据框架，包含三列:特征值、DualEigenvalues和实验。这个问题与DualEigenvalues数据无关，因此可以忘记。

我们进行了5个实验，并利用每个实验的数据计算了14个特征值。因此，该数据帧的前14行是第一实验的14个特征值，实验条目的值为1，第二14行为第二实验的14个特征值，实验条目的值为2等。

然后，我将绘制每个成对实验的特征值，下面是代码的一个示例：

eigen.1 <- mouse.data$Eigenvalues[mouse.data$Experiment == 1]
eigen.2 <- mouse.data$Eigenvalues[mouse.data$Experiment == 2]
p.data <- data.frame(x = eigen.1, y = eigen.2)
ggplot(p.data, aes(x,y)) + geom_abline(slope = 1, colour = "red") + geom_point()

这给了我这样的图表：

​

​

这正是我想要的这张图的样子。

我想要做的，但不能解决，是绘制一个facet_grid，这样在第1行和jth列中绘制了y轴的第i次实验的特征值和x轴上的jth实验的特征值。

这是我到目前为止最接近的，我希望这能使我的意思更清楚。

​

​

Answer 1

如果没有可重复的数据示例，这是很棘手的，但听起来我们可以大致近似于数据框架的结构，如下所示：

library(ggplot2)

set.seed(1)

Eigen <- as.vector(sapply(runif(5, .5, 1.5), 
                          function(x) sort(rgamma(14, 2, 0.02*x))))

mouse.data <- data.frame(Experiment = rep(seq(5), each = 14), Eigenvalue = Eigen)

head(mouse.data)
#>   Experiment Eigenvalue
#> 1          1   39.61451
#> 2          1   44.48163
#> 3          1   54.57964
#> 4          1   75.06725
#> 5          1   75.50014
#> 6          1   94.41255

要想让这个计划奏效，关键是将你的数据重组成一个包含每一个实验组合的长格式的数据框架。一种方法是通过split通过Experiment对数据帧进行索引，然后使用结果列表的简单索引(使用rep)来获得所有唯一的数据帧对。每个唯一的对按列方向粘在一起，然后得到的25个数据帧都按行方向连接到绘图数据帧中。

experiments <- split(mouse.data, mouse.data$Experiment)

experiments <- mapply(cbind, 
                      experiments[rep(1:5, 5)], 
                      experiments[rep(1:5, each = 5)], 
                      SIMPLIFY = FALSE)

p.data <- do.call(rbind, lapply(experiments, setNames, 
                                nm = c("Experiment1", "x",
                                       "Experiment2", "y")))

一旦我们这样做了，我们就可以使用您的情节代码，并添加一个facet_grid调用：

ggplot(p.data, aes(x,y)) + 
  geom_abline(slope = 1, colour = "red") + 
  geom_point() +
  facet_grid(Experiment1~Experiment2)

​

​