请问如何用R语言实现对n*p矩阵每一列求外积生成p*p矩阵，再把n个p*p矩阵加起来

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

除了循环之外有什么简单的方法吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：如何用 R语言

有点意思，搞了好久！

p=3
n=4
x=matrix(1:(n*p),nrow=p)
m_m <- function(x){
  y <- as.data.fr ame(as.matrix(x)%*%as.matrix(t(x)))
  }
(reult_list <- (result=apply(x,2,m_m)))
(output <-as.matrix(Reduce(`+`, reult_list)))

```{r}
p=3
n=4
x=matrix(1:(n*p),ncol=p)
m_m <- function(x){
  y <- as.data.frame(as.matrix(x)%*%as.matrix(t(x)))
  }
(result_list <- apply(x,2,m_m))
(output <-as.matrix(Reduce(`+`, result_list)))
```
## [[1]]
##   V1 V2 V3 V4
## 1  1  2  3  4
## 2  2  4  6  8
## 3  3  6  9 12
## 4  4  8 12 16
##
## [[2]]
##   V1 V2 V3 V4
## 1 25 30 35 40
## 2 30 36 42 48
## 3 35 42 49 56
## 4 40 48 56 64
##
## [[3]]
##    V1  V2  V3  V4
## 1  81  90  99 108
## 2  90 100 110 120
## 3  99 110 121 132
## 4 108 120 132 144
(output <-as.matrix(Reduce(`+`, result_list)))
##       V1  V2  V3  V4
## [1,] 107 122 137 152
## [2,] 122 140 158 176
## [3,] 137 158 179 200
## [4,] 152 176 200 224

caozhaowen 发表于 2020-3-5 00:15
```{r}
p=3
n=4

啊，疫情丧了一段时间，好久没看论坛了，才注意到您的解答，万分感谢。您提供了很棒很清晰的思路，这里我后来自己也想了一种思路，也分享一下吧

t <- apply(X, 1, function(x) x%*%t(x))
t <- apply(t, 2, sum )
t <- matrix(t, p, p)