Question

我有一个包含200多个数据点的矩阵。这是我的对象x。在第二个对象（元数据）中，我有一个包含20个数据点的列（y）。我想针对对象元数据中的20个数据点（y）绘制矩阵（对象x）

plot(x, metadata$y)

不起作用，因为x和y的长度不同。是否可以绘制此图？

矩阵x：

    X1  X4  X7  X9
X4  0.7                                                            
X7  0.8 0.5                                                  
X9  0.6 0.6 0.7

元数据

X1 65.4
X4 9.7
X7 47.4
X9 14.5

元数据$ y：65.4 9.7 47.4 14.5

Answer 1

这是一个tidyverse解决方案：好的，这些数据框（帖子底部的结构）：

> df
    X1  X4  X7 X9
X4 0.7  NA  NA NA
X7 0.8 0.5  NA NA
X9 0.6 0.6 0.7 NA

> metadata
      y
X1 65.4
X4  9.7
X7 47.4
X9 14.5

为简单起见，首先从rownames()中拉出metadata：

metadata$x <- rownames(metadata)

> metadata
      y  x
X1 65.4 X1
X4  9.7 X4
X7 47.4 X7
X9 14.5 X9

现在使用gather将x矩阵转换为长格式数据帧，并以'x'作为键，就像上面的rownames步骤一样。接下来，使用left_join将元数据连接到长数据帧，并将两个数据帧中的x作为公共列。

long <- gather(df, key = "x") %>%
  left_join(metadata, by = "x")

> long
    x value    y
1  X1   0.7 65.4
2  X1   0.8 65.4
3  X1   0.6 65.4
4  X4    NA  9.7
5  X4   0.5  9.7
6  X4   0.6  9.7
7  X7    NA 47.4
8  X7    NA 47.4
9  X7   0.7 47.4
10 X9    NA 14.5
11 X9    NA 14.5
12 X9    NA 14.5

图：

plot(value ~ y, data = long, pch = 19)

数据：

df <-
  structure(
    list(
      X1 = c(0.7, 0.8, 0.6),
      X4 = c(NA, 0.5, 0.6),
      X7 = c(NA,
             NA, 0.7),
      X9 = c(NA, NA, NA)
    ),
    .Names = c("X1", "X4", "X7", "X9"),
    row.names = c("X4", "X7", "X9"),
    class = "data.frame"
  )

metadata <-
  structure(
    list(y = c(65.4, 9.7, 47.4, 14.5)),
    .Names = "y",
    row.names = c("X1",
                  "X4", "X7", "X9"),
    class = "data.frame"
  )

Answer 2

看起来x并不是一个完整的矩阵，而是对称相异矩阵的下三角部分。 class(x)说"dist"吗？如果是这样，您可以使用as.matrix(x)将您的差异更改为对称矩阵，然后将其与y（元数据）的长度匹配。

这与另一个答案中所建议的技巧相同，不同之处在于，它为您提供了一个对称矩阵，而不是在对角线上有NA。您应该知道哪个是正确的决定。

R：对照一列中的数据绘制矩阵中的数据

2 个答案:

数据：