1

我想将大型成对矩阵的一部分子集化为较小的矩阵。例如

    TF1 TF2 TF3 TG1 TG2 TG3
TF1 0    2   1  450 460 450
TF2 2    0   1  452 462 462
TF3 1    2   0  451 461 451
TG1 450 452 450 0   2   0
TG2 460 462 462 2   0   1
TG3 450 452 451 1   2   0

如果我输入确切的列和行名称“TF1”等，我可以进行子集化，但我想要所有的 TF。

我认为它可能是一个 grep，但我不能将 grep 放在我的论坛中的任何位置：

TF <-T_PW[c("TF1","TF2","TF3"),c("TF1","TF2","TF3")]

这是预期的输出矩阵：

    TF2 TF3 TF1   
TF1 0   2   1
TF2 2   0   1
TF3 1   2   0

是不是因为它是成对的，所以列和行都有名称？

谢谢

4

1 回答 1

0

grep返回与提供的正则表达式匹配的输入索引：

> grep('^TF', rownames(T_PW))
[1] 1 2 3

> grep('^TF', colnames(T_PW))
[1] 1 2 3

因此，您可以在grep里面放置矩阵的和[]子集：rownamescolnames

TF = T_PW[grep('^TF', rownames(T_PW)), grep('^TF', colnames(T_PW))]

由于 与rownames对称colnames，您还可以执行以下操作：

TF_names = grep('^TF', colnames(T_PW))

T_PW[TF_names, TF_names]

结果：

    TF1 TF2 TF3
TF1   0   2   1
TF2   2   0   1
TF3   1   2   0

数据：

T_PW = read.table(text="    TF1 TF2 TF3 TG1 TG2 TG3
TF1 0    2   1  450 460 450
                 TF2 2    0   1  452 462 462
                 TF3 1    2   0  451 461 451
                 TG1 450 452 450 0   2   0
                 TG2 460 462 462 2   0   1
                 TG3 450 452 451 1   2   0", header = TRUE, row.names=1)

T_PW = as.matrix(T_PW)
于 2018-04-26T13:56:18.437 回答