2015-09-15 57 views
1

我有一个数字矩阵,我想从中检索给定特定值的索引。 我正在尝试which()函数来查找矩阵中的值。如何使用R中的which()函数从数字矩阵中检索值?

问题是发现了一些值,有些没有。

我的矩阵如下:

x_lat <- as.double(seq(48.0 ,60.0, by=0.1)) 
y_long <- as.double(seq(-10.0 ,2.0, by=0.1)) 
xv <- as.double(rep(x_lat,each = 121)) 
yv <- as.double(rep(y_long, 121)) 

vMatrix <- as.matrix(cbind(xv,yv)) 

如果我想要检索的索引,其中值-2.3是TRUE函数正确返回与指数-2.3地方出现的向量。

xx<- which(vMatrix==-2.3,arr.ind=TRUE) 

> xx 
    [1] 78 199 320 441 562 683 804 925 1046 1167 1288 1409 1530 1651 1772 1893 2014 2135 2256 2377 2498 
[22] 2619 2740 2861 2982 3103 3224 3345 3466 3587 3708 3829 3950 4071 4192 4313 4434 4555 4676 4797 4918 5039 
[43] 5160 5281 5402 5523 5644 5765 5886 6007 6128 6249 6370 6491 6612 6733 6854 6975 7096 7217 7338 7459 7580 
[64] 7701 7822 7943 8064 8185 8306 8427 8548 8669 8790 8911 9032 9153 9274 9395 9516 9637 9758 9879 10000 10121 
[85] 10242 10363 10484 10605 10726 10847 10968 11089 11210 11331 11452 11573 11694 11815 11936 12057 12178 12299 12420 12541 12662 
[106] 12783 12904 13025 13146 13267 13388 13509 13630 13751 13872 13993 14114 14235 14356 14477 14598 

但对于一些数字(出现在矩阵)的功能不起作用,例如,

xx<- which(vMatrix==-2.2,arr.ind=TRUE) 

> xx 
integer(0) 

任何帮助,将不胜感激

回答

0

浮点数字可能会产生误导。即使控制台可能显示相同的输出,两个这样的数字通常也不是“相等的”。机器只有一定的准确度才能表示数字。

这里有一个简单的例子:

a <- 0.15 - 1/8 
b <- 0.025 
> a 
[1] 0.025 
> b 
[1] 0.025 

但是,如果我们用“==”比较这些数字,我们得到:

> a==b 
[1] FALSE 

那是因为有从浮点产生的差异算术超出机器的精度:

> a-b 
[1] -6.938894e-18 

可能你可以解决t他通过简单的矩阵四舍五入号码相关数位的必要量,例如像问题,在矩阵

xx<- which(round(vMatrix,3)==-2.2,arr.ind=TRUE) 
+1

值只有一个小数,所以我不认为这是有道理他们舍入到一个更高的小数位数。我现在看到它是如何有意义的。谢谢! – armg

+0

不客气,我很高兴我能提供帮助。 – RHertel

+0

使用'round'不是比较数字与浮点容差的好方法。 'all.equal'实现更好的方式。更多细节:http://stackoverflow.com/a/9508558/1968 –