r - acos(1)返回NaN的某些值，而不是其他值

Question

我有一个经度和纬度值的列表，我试图找到它们之间的距离。使用标准的大圆方法，我需要找到:

acos(sin(lat1)*sin(lat2) + cos(lat1)*cos(lat2) * cos(long2-long1))

并将其乘以我使用的单位的地球半径。只要我们采用acos的值在[-1,1]范围内，这就是有效的。如果它们甚至略微超出此范围，即使差异是由于四舍五入，它也会返回 NaN。

我的问题是，有时候，当两个纬度/经度值相同时，这会给我一个 NaN错误。即使对于相同的一对数字，也并非总是如此，但列表中的数字总是相同。例如，我有一个人在沙漠的道路上停下来:

Time  |lat     |long
1:00PM|35.08646|-117.5023
1:01PM|35.08646|-117.5023
1:02PM|35.08646|-117.5023
1:03PM|35.08646|-117.5023
1:04PM|35.08646|-117.5023

当我计算连续点之间的距离时，例如，第三个值将始终为 NaN，即使其他值不是。 R舍入似乎是一个奇怪的错误。

Answer 1

在看不到您的数据的情况下无法准确分辨(尝试dput)，但这很可能是FAQ 7.31的结果。

(x1 <- 1)
## [1] 1
(x2 <- 1+1e-16)
## [1] 1
(x3 <- 1+1e-8)
## [1] 1
acos(x1)
## [1] 0
acos(x2)
## [1] 0
acos(x3)
## [1] NaN

也就是说，即使您的值如此相似以至于它们的打印表示形式相同，它们也可能会有所不同:有些将在 .Machine$double.eps中，而另一些将不会...

确保输入值受[-1,1]限制的一种方法是使用 pmax和 pmin: acos(pmin(pmax(x,-1.0),1.0))