R 简单编程

1个2乘100的矩阵
现在要求每一列的元素乘以它的转置，得到100个2乘2的矩阵
最后将得到的矩阵再求和
怎么编程？

简化的话 你看看这个如何

1. x = matrix(1:200,2,100)
   
2. dat=apply(x,2,function(i)i%*%t(i))
   
3. dat=apply(dat,1,sum)
   
4. matrix(dat,2,2)

复制代码

写好后已附件形式，我可以购买

1. x = matrix(rnorm(2*100),2,100)
   
2. tt = list()
   
3. for (i in 1:100){
   
4.   tt[[i]] = x[,i] %*% t(x[,i])
   
5. }
   
6. t=matrix(0,2,2)
   
7. for (i in 1:100) {
   
8.   t=tt[[i]]+t
   
9. }
   
10. t

复制代码

怎么用do.call命令简化呢

jmpamao 发表于 2014-7-4 11:01
简化的话 你看看这个如何

刚准备粘贴的，刷新一下发现你就回复了。速度。。。
话说命令跟你差不多

不过第二句  可化简为   dat=apply(x,2,tcrossprod)

yywan0913 发表于 2014-7-4 11:06
刚准备粘贴的，刷新一下发现你就回复了。速度。。。
话说命令跟你差不多

哈哈 学习了这句

jmpamao 发表于 2014-7-3 22:14
简化的话 你看看这个如何

十分感谢，谢谢！不知道怎么给你币？

famousid 发表于 2014-7-3 23:19

学习了，十分感谢

famousid 发表于 2014-7-3 23:19

怎么给你币呢？