1

我想使用第二个矩阵中的数据对一个矩阵的数据进行子集化。一个矩阵的列被标记。例如，

area1 <- c(9836374,635440,23018,833696,936079,1472449,879042,220539,870581,217418,552303,269359,833696,936079,1472449,879042,220539,870581, 833696,936079,1472449,879042,220539,870581)
id <- c(1,2,5,30,31,34,1,2,5,1,2,5,1,2,5,30,31,34,51,52,55,81,82,85)
mat1 <- matrix(area1, ncol=3, byrow=T)
mat2 <- matrix(id, ncol=3, byrow=T)
dimnames(mat1) <-list(NULL, c("a1","a2","a3"))

mat2包含 的 id ，因此矩阵的维度是mat1相同的（即mat1[1,1]标识, 子矩阵 2 和 3 各有 1 行, 子矩阵 4 将有 4 行数据。后续行与 1,3,5 之间的行数各不相同。这有意义吗？mat2[1,1]mat1c(1, 2, 5)mat2mat1

最初，矩阵是从数据帧转换而来的，其中 id 在一列，area 在第二列。我找不到在数据帧中的 1 行之间对变量行进行子集化的方法，这就是我切换到矩阵的原因。

4

4 回答 4

2

我认为这涵盖了它并符合您的描述：

spl <- cumsum(apply(mat2,1, function(x) all(x==c(1,2,5))))
split(as.data.frame(mat1),spl)

#$`1`
#       a1     a2      a3
#1 9836374 635440   23018
#2  833696 936079 1472449
# 
#$`2`
#      a1     a2     a3
#3 879042 220539 870581
#
#$`3`
#      a1     a2     a3
#4 217418 552303 269359
#
#$`4`
#      a1     a2      a3
#5 833696 936079 1472449
#6 879042 220539  870581
#7 833696 936079 1472449
#8 879042 220539  870581

结果符合“子矩阵 1 将有 2 行数据，子矩阵 2 和 3 各有 1 行，子矩阵 4 将有来自 mat1 的 4 行数据”

于 2013-09-11T22:31:44.947 回答
1 
mat1[which(mat2[,1]==1 & mat2[,2]==2 & mat2[,3]==5),]
        [,1]   [,2]    [,3]
[1,] 9836374 635440   23018
[2,]  879042 220539  870581
[3,]  217418 552303  269359
[4,]  833696 936079 1472449
于 2013-09-11T22:18:50.877 回答
1 
split(as.data.frame(mat1), apply(mat2, 1, paste, collapse = " "))
#$`1 2 5`
#       a1     a2      a3
#1 9836374 635440   23018
#3  879042 220539  870581
#4  217418 552303  269359
#5  833696 936079 1472449
#
#$`30 31 34`
#      a1     a2      a3
#2 833696 936079 1472449
#6 879042 220539  870581
#
#$`51 52 55`
#      a1     a2      a3
#7 833696 936079 1472449
#
#$`81 82 85`
#      a1     a2     a3
#8 879042 220539 870581
于 2013-09-11T22:26:09.347 回答
0

我认为从您所说的来看，您想将其保留为数据框。您可以通过抓取具有某些列值的行来轻松制作子矩阵。

在这里，我将数据框重新组合在一起，并为 1 创建了一个子矩阵。您可以通过在多个“area1”列上使用 cbind 等操作轻松添加到它。

> area1 <- c(9836374,635440,23018,833696,936079,1472449,879042,220539,870581,217418,552303,269359,833696,936079,1472449,879042,220539,870581, 833696,936079,1472449,879042,220539,870581)
> id <- c(1,2,5,30,31,34,1,2,5,1,2,5,1,2,5,30,31,34,51,52,55,81,82,85)
> original<-as.data.frame(cbind(id,area1))
> original[original$id==1,]
   id   area1
1   1 9836374
7   1  879042
10  1  217418
13  1  833696

然后你可以像我之前说的那样做。

> col1<-original[original$id==1,"area1"]
> col2<-original[original$id==2,"area1"]
> col3<-original[original$id==5,"area1"]
> submat<-cbind(col1,col2,col3)
> colnames(submat)<-c("a1","a2","a3")
> submat
          a1     a2      a3
[1,] 9836374 635440   23018
[2,]  879042 220539  870581
[3,]  217418 552303  269359
[4,]  833696 936079 1472449
于 2013-09-11T22:29:40.483 回答