r - 矩阵乘法(逐行)

Question

我有一个尺寸为1000x100的矩阵。我想逐行制作一个内部乘积(每行都有一个elselg)，因此从理论上讲，我可以得到一个1000x1的 vector 。例如:

A<-matrix(c(1,2,3,4),nrow=2,ncol=2,byrow=2)
     [,1] [,2]
[1,]    1    2
[2,]    3    4

我想要一个看起来像这样的 vector :

              [,1]
[1,]    c(1,2) %*% t(c(1,2))
[2,]    c(3,4) %*% t(c(3,4))  

我尝试做一个循环，但发生错误:

U<-matrix(nrow=1000,ncol=1)
U

k=0
for(i in 1:nrow(U_hat)){
  for(j in 1:nrow(U_hat)){
    k=k+1
    U[k,1]=U_hat[i,]%*%t(U_hat[j,])
  }
}

其中U_hat是尺寸为1000x100的矩阵。

我将很高兴知道如何进行乘法运算。谢谢。

Answer 1

自身乘以A并得出行总和:

rowSums(A*A)
## [1]  5 25

这也将起作用:

apply(A, 1, crossprod)
## [1]  5 25

这也可以工作:

diag(tcrossprod(A))
## [1]  5 25