Algorithm 用R语言从位置信息生成距离矩阵

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

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.

 }
}

for（1:7中的i）{
（i:7中的j）{
mymatrix[i，j]您正在寻找的是
dist
函数。有关详细信息，请参阅
？dist

我不清楚你所说的7乘7矩阵是什么意思，然后是元素[1,3]
指的是它们之间的距离（注意到没有5,6,7）。我认为这意味着你希望指的是
Can.Count
。你可以通过命名行和列并引用这些名称来实现

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


注：这将使用x-y坐标计算二维距离，
c（'XCoord'，'YCoord'）
。您可以使用
c（'XCoord'，'YCoord'，'ZCoord'）


dist\u matrix您要查找的是
dist
函数。有关详细信息，请参阅
？dist

我不清楚你所说的7乘7矩阵是什么意思，然后是元素[1,3]
指的是它们之间的距离（注意到没有5,6,7）。我认为这意味着你希望指的是
Can.Count
。你可以通过命名行和列并引用这些名称来实现

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


注：这将使用x-y坐标计算二维距离，
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