r - 如何对矩阵进行复杂的循环操作

Question

我有一个样本矩阵

 5 4 3  
 2 6 8   
 1 9 7   

我想要像这样的输出

max(5*6,5*8,5*9,5*7)  // i!=j condition  
max(4*2,4*8,4*1,4*7)  
max(3*2,3*6,3*1,3*9)  

等等...

计算后得到的这个最大值应该是矩阵形式。我需要概括它，因此我需要一个通用代码。

Answer 1

这可以完成工作，但这是一个非常缺乏想象力的解决方案，因为它只是循环执行请求的计算的行和列，而不是做任何矢量化的事情。

sapply(1:ncol(m), function(j) sapply(1:nrow(m), function(i) max(m[i,j]*m[-i,-j])))
#      [,1] [,2] [,3]
# [1,]   45   32   27
# [2,]   18   42   72
# [3,]    8   72   42

数据：

(m <- matrix(c(5, 2, 1, 4, 6, 9, 3, 8, 7), nrow=3))
#      [,1] [,2] [,3]
# [1,]    5    4    3
# [2,]    2    6    8
# [3,]    1    9    7