如何降低 R 中的机器精度 .Machine$double.eps？

Question

当使用名为 did 的 R 包时，我收到此错误：

Error in solve.default(preV) :    system is computationally singular:
reciprocal condition number = 4.09946e-19

solve 函数有一个名为tol 的参数，我可以通过它设置容差，从而避免问题。但是包没有给出这个选项来从它的顶层函数中指定这个tol，所以没有办法为solve设置tol。

我正在考虑另一种选择。由于solve() 默认使用.Machine$double.eps 作为tol，所以如果我可以减少这个默认值会很好。怎么样？

Answer 1

这是你不能减少的，因为它是一种 C 标准。但是，您可以手动将tol = 1e-20 设置为函数输入，例如，以覆盖默认容差。但是，这纯粹是一种数字解决方法。您需要考虑是否有可能在您的问题级别上解决您的排名不足问题。例如，由于数据变量之间的缩放不良可能会导致排名不足，这是一个简单的示例：Linear model singular because of large integer datetime in R?。

---

我不知道在这里做什么？正如我所说，我正在使用一个名为mp.spatt 的函数，其中使用了solve 函数。所以我真的不能（或者至少我不知道怎么做？）在solve 函数中添加tol = blabla 参数？

您可以编写该函数的补丁版本并重建软件包供您自己使用。这正是开源软件和软件包如此出色的原因。

转到https://cran.r-project.org/package=did，下载源文件.tar.gz for Linux，.tgz for Mac 或.zip for Windows。解压，打开/Rdirector/文件夹中的“did.R”文件。函数solve在这个文件中只使用了一次，所以你可以很容易地找到它：

W <- n*t(preatt)%*%solve(preV)%*%preatt

为简单起见，您可以添加tol = 0。或者，你们将solve 替换为MASS::ginv。实际上，包作者一直在脚本中的其他任何地方使用MASS::ginv。也可能是他忘记替换这个solve（我在2018/07/11查看最新版本1.1.0）。

修复后，再次生成.tar.gz、.tgz或.zip文件。打开你的 R，使用install.packages 安装这个特定的文件，然后愉快地使用它。