按名称从矩阵中删除列和行 [R]

Question

我有一个亲属关系矩阵，它本质上是一个向下的 ID 列表和相同的 ID，如下所示：

ID1  ID2  ID3  ID4  ID5  ID6
ID2
ID3
ID4
ID5
ID6

等等。我想从矩阵中删除某些 ID（即删除我不想要的 ID 的列和行）。我的矩阵是 129 x 129，我想要保留的 ID 列表是 60 个 ID。我尝试了以下代码（带有虚拟数据）：

kinstmp <- matrix(c(0,0,1,0,1,0,
     0,0,1,0,0,0,
     0,1,1,1,0,0,
     1,0,1,0,0,0,
     0,0,1,0,1,0,
     0,0,1,0,0,0,
     0,1,1,1,0,0), nrow=6) #example of what my matrix looks like
kinstmp2 <- kinstmp[rownames(kinstmp) %in% id, ]
kins <- kinstmp2[,colnames(kinstmp2) %in% id ]

但是，这不会导致任何行和列被删除。有没有办法让它工作？

id文件的一个例子是：

ID3
ID5
ID6

Answer 1

您需要将dimnames 添加到矩阵中，并确保将列条件放在逗号的右侧（字面意思）。

nm <- paste0("ID", 1:6)
kinstmp <- matrix(c(0,0,1,0,1,0,
                    0,0,1,0,0,0,
                    0,1,1,1,0,0,
                    1,0,1,0,0,0,
                    0,0,1,0,1,0,
                    0,0,1,0,0,0), nrow = 6, dimnames = list(nm, nm))
kinstmp

id <- c("ID3", "ID5", "ID6")

# subsetting rows
kinstmp2 <- kinstmp[!rownames(kinstmp) %in% id, ]
kinstmp2
#     ID1 ID2 ID3 ID4 ID5 ID6
# ID1   0   0   0   1   0   0
# ID2   0   0   1   0   0   0
# ID4   0   0   1   0   0   0

# subsetting columns
kins <- kinstmp2[ , !colnames(kinstmp2) %in% id]
kins
#     ID1 ID2 ID4
# ID1   0   0   1
# ID2   0   0   0
# ID4   0   0   0

# or subsetting rows and columns in one call
kinstmp[!rownames(kinstmp) %in% id, !colnames(kinstmp2) %in% id]