2015-09-15 55 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 –