R:从另一个矩阵的元素中减去矩阵的元素

时间:2013-11-26 13:45:06

标签: r matrix equation-solving

我正在尝试将方程应用于两个矩阵。由于我是初学R用户,对我来说似乎很难。如果你能给我一些建议,我会很高兴。

我有两个相似性矩阵:

> r
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    0    4    2    2    5    5
[2,]    4    0    8    3    3    2
[3,]    2    8    0    4    4    3
[4,]    2    3    4    0    0    3
[5,]    5    3    4    0    0    5
[6,]    5    2    3    3    5    0

> nr
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    0    4    7    2    4    3
[2,]    4    0    5    2    3    2
[3,]    7    5    0    3    2    2
[4,]    2    2    3    0    7    2
[5,]    4    3    2    7    0    5
[6,]    3    2    2    2    5    0

我愿意申请以下内容:

sum((r[i,j]-nr[i,j])^2)/6

我的首要问题是从元素nr中提取r的元素。 如果我用nr[i,j]替换0.4,例如s<-numeric() for (i in 1:nrow(r)) { for (j in 1:ncol(r)) { s[k]<-sum((r[i,j]-0.4)^2)/6 } } > s [1] 0.02666667 ,则以下内容完美无缺:

{{1}}

但我无法弄清楚如何修改此代码以解决原始问题。 我将不胜感激任何帮助/建议。 谢谢!

1 个答案:

答案 0 :(得分:7)

+-*/^等常规运算符执行元素操作。所以简单地(r - nr)^2/6就可以为你做到这一点。

r
##      [,1] [,2] [,3]
## [1,]    2    2    2
## [2,]    2    2    2
## [3,]    2    2    2

nr
##      [,1] [,2] [,3]
## [1,]    3    3    3
## [2,]    3    3    3
## [3,]    3    3    3


r * nr
##      [,1] [,2] [,3]
## [1,]    6    6    6
## [2,]    6    6    6
## [3,]    6    6    6


r - nr
##      [,1] [,2] [,3]
## [1,]   -1   -1   -1
## [2,]   -1   -1   -1
## [3,]   -1   -1   -1


(r - nr)^2/6
##           [,1]      [,2]      [,3]
## [1,] 0.1666667 0.1666667 0.1666667
## [2,] 0.1666667 0.1666667 0.1666667
## [3,] 0.1666667 0.1666667 0.1666667