我正在尝试将方程式应用到两个矩阵。由于我是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
但是我不知道如何修改此代码来解决原始问题。我将不胜感激。谢谢!
正常的运营商喜欢+
,-
,*
,/
和^
做的元素方式操作。因此,简单地(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
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句