1

假设我有一个如下形式的矩阵：

Residue Can.Count SideChain XCoord  YCoord ZCoord
1       MET         1         A 62.935  97.579 30.223
2       THR         2         A 63.155  95.525 27.079
3       GLU         3         A 65.289  96.895 24.308
4       TYR         4         A 64.899  96.220 20.615
8       LYS         8         A 67.593  96.715 18.023
9       LEU         9         A 65.898  97.863 14.816
10       VAL        10         A 67.664  98.557 11.533

请注意，数字 5-6-7 被跳过。我想要做的是在每个残基之间建立一个“距离矩阵”。在这种情况下，我想制作一个 7x7 矩阵，其中元素 (1,3) 是这些位置之间的距离。

现在，我意识到我不需要填写下半部分，对角线以上的所有内容就足够了。我还看到了如何使用 2 个 for 循环来做到这一点，如下所示：

 for(i in 1:7) {
  for(j in i:7){
    mymatrix[i,j] <- calcdistance(xyz1,xyz2) #I have the distance function already coded.

 }
}

我意识到它总是 O(n^2) 但我想知道我是否可以利用 R 的力量使用 apply 语句（或更聪明的东西）来制作这个矩阵？我试过这样做，但不知何故没有成功。谢谢你的帮助！

4

1 回答 1

4

您正在寻找的是dist功能。详情请参阅?dist。

我不清楚您期望 7 x 7 矩阵是什么意思，然后元素 [1,3] 指的是它们之间的距离（在注意到没有 5,6,7 之后）。我认为它的意思是您希望参考Can.Count. 您可以通过命名行和列并引用这些名称来做到这一点。

假设您的数据是一个名为的 data.frame residues，以下将起作用

  • 请注意，这使用 xy 坐标计算二维距离， c('XCoord','YCoord')。您可以使用 c('XCoord','YCoord', 'ZCoord').
dist_matrix <-  as.matrix(dist(residues[, c('XCoord','YCoord')], diag = T))
# this gives a 7 by 7 matrix
dist_matrix
##         1        2         3         4        5        6        7
## 1 0.000000 2.065748 2.4513613 2.3883419 4.737453 2.976579 4.829071
## 2 2.065748 0.000000 2.5359132 1.8773814 4.594774 3.604205 5.433609
## 3 2.451361 2.535913 0.0000000 0.7795672 2.311021 1.143637 2.898770
## 4 2.388342 1.877381 0.7795672 0.0000000 2.739099 1.922875 3.620331
## 5 4.737453 4.594774 2.3110206 2.7390986 0.000000 2.047176 1.843368
## 6 2.976579 3.604205 1.1436367 1.9228755 2.047176 0.000000 1.897470
## 7 4.829071 5.433609 2.8987703 3.6203306 1.843368 1.897470 0.000000

# set the dimension names to the Can.Count so we can refer to them
dimnames(dist_matrix) <- list(residues[['Can.Count']],residues[['Can.Count']] )

# now you can refer to the distance between Can.Count 1 and Can.Count 8
dist_matrix['1','8']

## [1] 4.737453

# note that you need to refer to the dimension names as characters, 
# as this is  7 by 7 matrix, so the following will give 
# an (obvious) error message
dist_matrix[1,8]

## Error: subscript out of bounds
于 2012-08-14T23:03:37.917 回答