矩阵提取

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有一个matrix和两个vector，比如：
matrix
    1  2  3  4  5  
1  4 -1 -2 -2  0
2 -1  5  0 -2 -3   
3 -2  0  6  1 -3
4 -2 -2  1  6 -3
5  0 -3 -3 -3  9


vector1：1 2 3 4 5 6
vector2：3 4 1 5 2 1

需要让两个vector一一对应取出矩阵的数字，即把matrix[1,3]、[2,4]、[3,1]……这样提出来，但是又不想用for或while这种循环，有什么简便点的方法么？

还有个问题，如果矩阵的header不是数字而是字母，vector也是字母不是数字，这种矩阵要怎么提取？

谢谢~

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Vector matrix while Head ect matrix

diag(matrix[vector1,vector2])

你的vector1里面有数字已经越界了。还有你想得到的输出是一个向量么？若是，设你待提取的矩阵名为x, 两个向量为v1, v2 (必须要等长度)，以下代码给出所要结果 （楼上的代码似乎有问题):

1. ind <- cbind(v1, v2)
   
2. xind <- x[ind]

复制代码

hugebear 发表于 2015-4-11 12:24
你的vector1里面有数字已经越界了。还有你想得到的输出是一个向量么？若是，设你待提取的矩阵名为x, 两个向 ...

啊举例的时候没注意不好意思
原来这样就可以了，谢谢~~

“如果矩阵的header不是数字而是字母，vector也是字母不是数字，这种矩阵要怎么提取？”这个问题不知道你可以解答么？

syfzxy 发表于 2015-4-11 12:06
diag(matrix[vector1,vector2])

好像不行呃……thx anyway~

waxxyeah 发表于 2015-4-11 13:11
好像不行呃……thx anyway~

我试了一下可以的，能告诉我是哪儿的问题吗？一起学习下

syfzxy 发表于 2015-4-11 13:23
我试了一下可以的，能告诉我是哪儿的问题吗？一起学习下

Error in diag(matrix[vector1, vector2]) :
  (list) object cannot be coerced to type 'double'

waxxyeah 发表于 2015-4-11 14:38
Error in diag(matrix[vector1, vector2]) :
  (list) object cannot be coerced to type 'double'

好吧，你那个变量是list，谢谢了

waxxyeah 发表于 2015-4-11 13:10
啊举例的时候没注意不好意思
原来这样就可以了，谢谢~~

只要矩阵的dimnames(想必是你所说的header)和vector所含的字母一致，上述的方法依然有效，举例如下：

1. > x <- matrix(rnorm(5 * 5), 5, 5)
   
2. > dimnames(x) <- list(letters[1:5], letters[1:5])
   
3. > v1 <- sample(letters[1:5], size = 6, replace = T)
   
4. > v2 <- sample(letters[1:5], size = 6, replace = T)
   
5. xind <- cbind(v1, v2)
   
6. x[xind]

复制代码

hugebear 发表于 2015-4-11 21:46
只要矩阵的dimnames(想必是你所说的header)和vector所含的字母一致，上述的方法依然有效，举例如下：

原来matrix也可以用名字提取的，学习了！谢谢~