【问题标题】:Merge matrices from one list on matrices from another one in R将一个列表中的矩阵合并到 R 中另一个列表中的矩阵上
【发布时间】:2015-02-03 10:42:40
【问题描述】:

在看到 (Merge matrices from two lists in R) 之前,我已经问过一个类似的问题。但是,与我之前的问题相反,在这种情况下:

  1. 我在 list1 中的矩阵与 list2 中的标准矩阵具有更小、更大或相同大小的维度。
  2. list1 中的某些矩阵包含未包含在 list2 中的案例。

在下面,您将找到详细的问题:

假设我有两个包含多个矩阵的列表。第一个列表包含的矩阵的维度因矩阵而异:

创建list1的代码:

d<-c(0,1,0,1)
e<-c(1,0,0,0)
f<-c(0,0,0,0)
g<-c(1,0,0,0)
cn<-c(1,2,7,4)
p<-data.frame(d,e,f,g)
dimnames(p)<-list(cn,cn)
p<-as.matrix(p)

d<-c(0,0,0,1,1)
e<-c(0,0,1,0,0)
f<-c(0,1,0,0,0)
g<-c(1,0,0,0,0)
h<-c(1,0,0,0,0)
cn<-c(1,5,3,2,4)
q<-data.frame(d,e,f,g,h)
dimnames(q)<-list(cn,cn)
q<-as.matrix(q)

d<-c(0,1,0,1,0,0)
e<-c(1,0,0,0,0,0)
f<-c(0,0,0,0,0,1)
g<-c(1,0,0,0,1,0)
h<-c(0,0,0,1,0,0)
i<-c(0,0,1,0,0,0)
cn<-c(1,2,3,6,8,9)
r<-data.frame(d,e,f,g,h,i)
dimnames(r)<-list(cn,cn)
r<-as.matrix(r)

list1<-list(p,q,r)
names(list1)<-1990:1992

列表1:

$`1990`
  1 2 7 4
1 0 1 0 1
2 1 0 0 0
7 0 0 0 0
4 1 0 0 0

$`1991`
  1 5 3 2 4
1 0 0 0 1 1
5 0 0 1 0 0
3 0 1 0 0 0
2 1 0 0 0 0
4 1 0 0 0 0

$`1992`
  1 2 3 6 8 9
1 0 1 0 1 0 0
2 1 0 0 0 0 0
3 0 0 0 0 0 1
6 1 0 0 0 1 0
8 0 0 0 1 0 0
9 0 0 1 0 0 0

第二个列表包括始终具有相同维度的矩阵。 list1 中的矩阵可以比 list2 中的矩阵更小、更大或大小相同。此外,list2 中的标准矩阵有时并不包含 lst1 中各个矩阵的所有情况。

产生list2的代码:

o<-matrix(NA,nrow=5,ncol=5)
dimnames(o)<-list(1:5, 1:5)
list2<-list(o,o,o)
names(list2)<-1990:1992

列表2:

$`1990`
   1  2  3  4  5
1 NA NA NA NA NA
2 NA NA NA NA NA
3 NA NA NA NA NA
4 NA NA NA NA NA
5 NA NA NA NA NA

$`1991`
   1  2  3  4  5
1 NA NA NA NA NA
2 NA NA NA NA NA
3 NA NA NA NA NA
4 NA NA NA NA NA
5 NA NA NA NA NA

$`1992`
   1  2  3  4  5
1 NA NA NA NA NA
2 NA NA NA NA NA
3 NA NA NA NA NA
4 NA NA NA NA NA
5 NA NA NA NA NA

我想要做的是将list2 中的NA 替换为list1 中相应矩阵的值(如果可用)。 list2 中矩阵的维度在此过程中应保持不变。结果应如下所示:

$`1990`
   1  2  3  4  5 
1  0  1 NA  1 NA
2  1  0 NA  0 NA
3 NA NA NA NA NA
4  1  0 NA  0 NA
5 NA NA NA NA NA

$`1991`
   1  2  3  4  5
1  0  1  0  1  0
2  1  0  0  0  0
3  0  0  0  0  1
4  1  0  0  0  0
5  0  0  1  0  0

$`1992`
   1  2  3  4  5
1  0  1  0  NA NA
2  1  0  0  NA NA
3  0  0  0  NA NA
4  NA NA NA NA NA
5  NA NA NA NA NA

非常感谢任何帮助!

【问题讨论】:

    标签: r list matrix


    【解决方案1】:

    Map 仍然是要走的路,只需更新函数并按colnamesrow.names 对您的第一个列表进行排序:

    list11 = lapply(list1, function(u){
                m=data.matrix(u)
                m[order(as.integer(row.names(m))),order(as.integer(colnames(m)))]
         })
    
    f = function(A,B)
    {
        A[row.names(A) %in% row.names(B), colnames(A) %in% colnames(B)]=B[row.names(B) %in% row.names(A), colnames(B) %in% colnames(A)]
        A
    }
    
    Map(f, list2, list11)
    
    #$`1990`
    #   1  2  3  4  5
    #1  0  1 NA  1 NA
    #2  1  0 NA  0 NA
    #3 NA NA NA NA NA
    #4  1  0 NA  0 NA
    #5 NA NA NA NA NA
    
    #$`1991`
    #  1 2 3 4 5
    #1 0 1 0 1 0
    #2 1 0 0 0 0
    #3 0 0 0 0 1
    #4 1 0 0 0 0
    #5 0 0 1 0 0
    
    #$`1992`
    #   1  2  3  4  5
    #1  0  1  0 NA NA
    #2  1  0  0 NA NA
    #3  0  0  0 NA NA
    #4 NA NA NA NA NA
    #5 NA NA NA NA NA
    

    【讨论】:

    • 我的方法类似,但我在处理它时发布了它,而且您的解决方案的某些元素与预期输出不匹配
    • 嗯,有一个问题,因为第一个列表中的 colnames 和 row.names 是无序的。解决了!
    • 主要功能通常不是大问题,而是输入中的数据^ ^无论如何感谢您的鹰眼!
    【解决方案2】:

    这是与@Colonel Beauvel 的解决方案类似的方法,但使用match

      f1 <- function(x,y) {rInd1 <- match(row.names(x), row.names(y))
               rInd2 <- match(row.names(y), row.names(x))
               cInd1 <- match(colnames(x), colnames(y))
               cInd2 <- match(colnames(y), colnames(x))
               x[rInd1[!is.na(rInd1)], cInd1[!is.na(cInd1)]] <- 
                       y[rInd2[!is.na(rInd2)], cInd2[!is.na(cInd2)]]
           x}
    
      Map(f1, list2, list1)
      # $`1990`
      #   1  2  3  4  5
      #1  0  1 NA  1 NA
      #2  1  0 NA  0 NA
      #3 NA NA NA NA NA
      #4  1  0 NA  0 NA
      #5 NA NA NA NA NA
    
     #$`1991`
     #  1 2 3 4 5
     #1 0 1 0 1 0
     #2 1 0 0 0 0
     #3 0 0 0 0 1
     #4 1 0 0 0 0
     #5 0 0 1 0 0
    
     #$`1992`
     #  1  2  3  4  5
     #1  0  1  0 NA NA
     #2  1  0  0 NA NA
     #3  0  0  0 NA NA
     #4 NA NA NA NA NA
     #5 NA NA NA NA NA
    

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 1970-01-01
      • 2014-05-15
      • 2019-01-27
      • 1970-01-01
      • 2014-01-26
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多