R 语言中 package 'Brobdingnag' 的用法

Question

我遇到了有关使用 Brobdingnag package 的问题 - 设置后

a2 <- as.brob(0.1)^1000, 

a2 = exp(-2302.6)

a1 <- as.brob(0.1)^800, 

a1 = exp(-1842.1)

使用 sum(a1,a2) 和 sum(a2,a1) 得到不同的结果 - 每次结果等于给和 function 的第一个参数。 似乎总和没有被 Brobdingang package 覆盖，即使它应该是？ 或者也许我做错了什么？

我也问了这个问题作为对我写的另一个问题的答复，请参见此处

[编辑：包作者的回答]

嗨丹

这绝对是 package 中的错误； 谢谢你的报告，不幸的是。 纠正它需要我相当长的时间。

同时，请在下面找到用于计算两个 brobs 之和的常用 R 习语：

> a1 <- as.brob(0.1)^800
> a2 <- as.brob(0.1)^1000
> a1+a2


> a1 <- as.brob(0.1)^800
> a2 <- as.brob(0.1)^1000


> a1+a2
[1] +exp(-1842.1)
> a2+a1
[1] +exp(-1842.1)

> sum(cbrob(a1,a2))
[1] +exp(-1842.1)
> sum(cbrob(a2,a1))
[1] +exp(-1842.1)
> 

Answer 1

我可以使用以下代码重现您的问题。 一个答案可能是+而不是sum 。 补充：另一个是在求和之前将您的数据放在一个向量中

> library(Brobdingnag)
> (a1 <- as.brob(0.1)^800) 
[1] +exp(-1842.1)
> (a2 <- as.brob(0.1)^1000)
[1] +exp(-2302.6)
>
> a1 + a2
[1] +exp(-1842.1)
> a2 + a1
[1] +exp(-1842.1)
>
> sum(a1, a2)
[1] +exp(-1842.1)
> sum(a2, a1)
[1] +exp(-2302.6)
> 
> sum(as.brob(0.1)^c(1000,800))
[1] +exp(-1842.1)
> sum(as.brob(0.1)^c(800,1000))
[1] +exp(-1842.1)

似乎您不能像这样使用sum(,) 。 这里有一些类似的奇怪结果，有更实用的数字

> as.brob(0.1) + as.brob(1)           # OK, gives exp(ln(1.1))
[1] +exp(0.09531)     
> as.brob(1) + as.brob(0.1)           # OK, gives exp(ln(1.1))
[1] +exp(0.09531)
> sum(as.brob(c(0.1, 1)))             # OK, gives exp(ln(1.1))
[1] +exp(0.09531)
> sum(as.brob(c(1, 0.1)))             # OK, gives exp(ln(1.1))
[1] +exp(0.09531)
>
> sum(as.brob(0.1), as.brob(1))       # not OK, gives first term exp(ln(0.1))
[1] +exp(-2.3026)
> sum(as.brob(1), as.brob(0.1))       # not OK, gives first term exp(ln(1))
[1] +exp(0)