R语言中怎么定义矩阵序列，并将n个矩阵相加

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如题~

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R语言

举个例子

bbslover 发表于 2016-7-27 21:08
举个例子

比如一年360天，每一天有43行数据，我们需要把每43行的数据看做一个矩阵，这样就有了360个43*43的矩阵。然后再以30天为一个月，每30个矩阵相加，就得到12个月的矩阵。。这中间怎么定义矩阵序列和n个序列求和，都遇到问题了
D=array(0,dim=c(43,43,360))
for(j in 1:365){                   #将每一天弄成一个矩阵，每一天是43行
    p[j]=43*j-42  
    q[j]=43*j
    D[j]=rbind(x[c(p[j]:q[j]),])  #第j天的矩阵

下面是求和
for(k in 1:12)    {                  #一共有12个月
   g[k]=sum(a[c(1:k)])-a[k]+1
   #计算出每个月的第一天和最后一天的时间
   d[k]=sum(a[c(1:k)])
   M[,,k]=apply(D[,,seq(g[k],d[k],by=1)],c(1,2),sum)       #将每个月中所有天数的矩阵相加


这里的D定义是不是对的？是应该写成D[,,j]还是D[,j]还是怎么办？运行的时候老出错

另外，求和的部分有问题吗？

感觉最主要的是怎么定矩阵序列

急求答案或者程序啊！谢谢各位大神了！

使用list看看

dqj1551007438 发表于 2016-7-27 21:53
急求答案或者程序啊！谢谢各位大神了！

1. L <- list()
   
2. 
   
3. for (i in 1:360) {
   
4.   L[[i]] <- matrix(0, 43, 43)
   
5. }
   
6. 
   
7. length(L)
   
8. 
   
9. monthSum <- list()
   
10. 
    
11. for (i in 1:12) {
    
12.   cur <- i * 30
    
13.   bef <- (i - 1) * 30 + 1
    
14.   if (i == 1) {
    
15.     theVec <- 1:cur
    
16.   } else {
    
17.     theVec <- bef:cur
    
18.   }
    
19.   monthSum[[i]] <- Reduce("+", L[theVec])
    
20. }
    
21. 
    
22. length(monthSum)

复制代码

看看能不能从这个简单的code中得到启发！