Question

我正在尝试使用一个名为InfoTrad的软件包。我试图计算一个给定的示例，但它不会为我产生相同的结果。计算能力可能是这个原因吗？我得到了一堆NaN，这可能是导致错误结果的原因。代码如下：

library(InfoTrad)
Buy <- c(350,250,500,552,163,345,847,923,123,349)
Sell <- c(382,500,463,550,200,323,456,342,578,455)
data = cbind(Buy,Sell)
#par0 = (al,del,mu,eps_b,eps_s)
par0 = c(0.5, 0.5, 300, 400, 500)

EHO_out = EHO(data)
model = optim(par0, EHO_out, gr = NULL, method = c("Nelder-Mead") , hessian = FALSE)
print(model)

#output
$par
[1] -216494.389  631863.422     380.386     307.447     435.555

$value
[1] -44532.65

真实的结果（来自示例）应如下所示：

$par
[1] 0.9111102  0.0001231429     417.1497     336.075     466.2539

这是论文的链接，他们在其中解释了程序包：https://www.semanticscholar.org/paper/InfoTrad%3A-An-R-package-for-estimating-the-of-%C3%87elik-Tini%C3%A7/9efb3f9de4781bed286e16e9f4e2a87653bea693

是否可以更改极限值，以使我不产生任何NaN或至少不减少？先感谢您。如果您可以运行一次代码并告诉我，如果您的结果与我的不同，我将不胜感激。（我只是想知道，这是否是我的计算机的故障）

Answer 1

显然，问题出在包装内（我已经在具有不同计算能力的多台计算机上测试了代码）。遗憾的是，我不知道它到底在哪里，但是我发现了另一个名为pinbasic的软件包，它可以更好地工作并且具有更多的输出选项：

https://cran.rstudio.com/web/packages/pinbasic/vignettes/pinbasicVignette.html

由于NaN导致的结果错误（软件包InfoTrad）

1 个答案: