用R实现一个矩阵的计算

A <-c(1,2,3)
B <-c(2,3,4)
C <-c(4,5,6)
D <-list(A,B,C)
E <- matrix(rep((D), 4), nrow=3, ncol=4)
E
得到E如下:

现在要分别从V1的第1，2，3行 每一行随机取一个数，相加 得到一个和，这样执行3次(length(c(1,2,3))), 得到三个和，组成一个新的包含三个元素的Vector. 对于V2, V3, V4来说一样。 最后得到的结果应该是四个向量，每个向量包含三个数值元素。

这只是个简化的例子，实际情况是 "c(1,2,3) "不是三个元素，而是100000个。

多谢！

1. v1=c()
   
2. for(i in 1:3){
   
3. v1=c(v1,sum(unlist(lapply(1:3,function(x)sample(unlist(E[x,1]),1)))))
   
4. }
   
5. v1

复制代码

你再优化下吧

yywan0913 发表于 2014-2-1 23:25
你再优化下吧

多谢，
有个问题，i 似乎没用到。把 sample(unlist(E[x,1]),1) 里面的1  变成 i ?

sty1882 发表于 2014-2-3 22:37
多谢，
有个问题，i 似乎没用到。把 sample(unlist(E[x,1]),1) 里面的1  变成 i ?

1. v1=c();v2=c();v3=c();v4=c()
   
2. for(i in 1:3){
   
3. 
   
4. v1=c(v1,sum(unlist(lapply(1:3,function(x)sample(unlist(E[x,1]),1)))))
   
5. 
   
6. v2=c(v2,sum(unlist(lapply(1:3,function(x)sample(unlist(E[x,2]),1)))))
   
7. 
   
8. v3=c(v3,sum(unlist(lapply(1:3,function(x)sample(unlist(E[x,3]),1)))))
   
9. 
   
10. v4=c(v4,sum(unlist(lapply(1:3,function(x)sample(unlist(E[x,4]),1)))))
    
11. 
    
12. }
    
13. 
    
14. v=data.frame(v1=v1,v2=v2,v3=v3,v4=v4)

复制代码

就这样的，i用来执行三次的。等想到好的函数可以优化下