【发布时间】:2013-11-26 13:45:06
【问题描述】:
我正在尝试将方程应用于两个矩阵。由于我是 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
我最大的问题是从元素r 中提取nr 的元素。
如果我将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
但我不知道如何修改此代码以解决原始问题。 我将不胜感激任何帮助/建议。 谢谢!
【问题讨论】:
-
你试过简单的
sum((r-nr)^2)/6吗? -
不,我没有。我没想到会这么简单。谢谢!
标签: r matrix equation-solving