R矩阵排序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

例如我有个矩阵A
A=(1，3，
      2，3，
      4， 5)

现在打乱A的行的顺序
B=(2，3
      4， 5，
      1， 3)
如何找出B的矩阵的每一行在原来A矩阵的顺序，结果应该是(2,3,1)即原来的第一行在A矩阵中是第2行所以第一个元素是2，依次类推。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：r矩阵

1. apply(apply(B, 1, function(y) {apply(A, 1, function(x) identical(x, y))}), 2, which)

复制代码

cheetahfly 发表于 2017-8-31 16:23

谢谢，亲测之后，举了一个特殊的例子，觉得不行呢。比方说，
a <- matrix(c(1:16),4,4)
A <- which(upper.tri(a),arr.ind = TRUE)
b <- c(10,5,9, 13, 14,15)
G1 <- b %% ncol(a)
G2 <- floor(b / ncol(a)) + 1
B <- interactionOrder <- cbind(G1,G2)

colnames(A) <- NULL
colnames(B) <- NULL

apply(apply(B, 1, function(y) {apply(A, 1, function(x) identical(x, y))}), 2, which)
会报错。

子衿1219 发表于 2017-8-31 17:09
谢谢，亲测之后，举了一个特殊的例子，觉得不行呢。比方说，
a

问题出在identical()太严格了

1. apply(apply(B, 1, function(y) {apply(A, 1, function(x) all(dplyr::near(x, y)))}), 2, which)

复制代码

cheetahfly 发表于 2017-8-31 17:42
问题出在identical()太严格了

非常感谢。可以了~~~

cheetahfly 发表于 2017-8-31 17:42
问题出在identical()太严格了

刚又把这个程序用在大矩阵中比方说，5000*2,还是会显示跟刚才一样的报错。。

apply(B, 1, function(x){which.max(!colSums(t(A) != x))})可以解决。

我的方法中用了三重循环，不适合大矩阵，慢的要死，你的方法快多了！
如果针对这种“瘦长”型的大矩阵，我会转向用data.frame+left_join的思路：

1. library(dplyr)
   
2. tbl_df(B) %>%
   
3.   left_join(mutate(tbl_df(A), order = row_number()), by = names(.)) %>%
   
4.   pull(order)

复制代码

这样做的好处是，大矩阵中，如果有重复或出错的行，都不影响运行。过后检验结果便知。

cheetahfly 发表于 2017-8-31 20:54
我的方法中用了三重循环，不适合大矩阵，慢的要死，你的方法快多了！
如果针对这种“瘦长”型的大 ...

非常感谢~~~

#先并起来
A <- apply(A,1,paste,sep="_")
B <- apply(B,1,paste,sep="_")
#再grep
result <- sapply(B, function(x) grep(x,A))