- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我不小心发现了 R base 中的 as.integer
或 det
函数的奇怪错误。有谁知道这里发生了什么以及如何预防?
我正在计算以下 3×3 矩阵的行列式:
mat <- matrix(c(15, 6, 116, 10, 13, 16, 14, 23, 56), ncol = 3)
看起来像这样:
[,1] [,2] [,3]
[1,] 15 10 14
[2,] 6 13 23
[3,] 116 16 56
有两件事很容易看出:所有条目都是整数,并且六组三条目中的每一组都至少包含在同一行或列中的两个条目至少一个偶数。因此行列式必须是偶数。
通过键入 det(mat)
向 R 询问此行列式的实际值,它返回看起来像偶数的东西:8952
。但是你瞧:在 R 的内心深处,它实际上是一个非整数或奇数,因为在键入 as.integer(det(mat))
时我们得到 8951
.
这是怎么回事? 8951显然是错误的。此外,不太明显的是,值 8952 是正确的,可以用笔和纸看出。
所以我的问题是:
这是怎么回事?
当被要求计算整数矩阵的行列式时,我如何强制 R 给我正确整数值?
最佳答案
根本原因:is.integer
截断而不是舍入和 float 学记录的中间值来解释第二个结果,结合 print< 中显示的默认数字级别
控制台 REPL 的一部分,用于解释 det(mat)
的初始结果:
print( det(mat), digits =16)
[1] 8951.999999999993
理论上的答案很可能是 8952,但 R 不是符号数学引擎。
您可以使用 Rmpfr 包(如@BenBolker 所建议的那样)来提高精度级别:
library(Rmpfr)
mat <- mpfr(mat, 64)
as.integer( det(mat) )
[1] 8952
as.integer
截断而不是舍入。参见 ?as.integer
。 R 可以在不损失精度的情况下处理整数的加法或乘法,但一旦发生除法,就可能会出现浮点错误。 (实际上,问题的出现是因为 det
的默认设置是使用 determinant
和 log=TRUE
,然后对复数结果的模取幂。)从帮助页面的值部分:
Non-integral numeric values are truncated towards zero (i.e.,
as.integer(x)
equalstrunc(x)
there), and imaginary parts of complex numbers are discarded (with a warning).
关于r - as.integer(8952) = 8951?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/65401271/
我正在尝试执行 JavaPairRDD 和 JavaPairRDD 的 leftOuterJoin> 并且函数签名返回类型是 JavaPairRDD>>> 这里可选的是 com.google.comm
我正在尝试按元素的频率对元素进行排序 import java.io.BufferedReader; import java.io.IOException; import java.io.InputSt
这个问题已经有答案了: Is List a subclass of List? Why are Java generics not implicitly polymorphic? (19 个回答) 已
编辑:问题已解决:请参阅 Karim SNOUSSI 的答案和我在下面的评论。 这是我在堆栈溢出时遇到的第一个问题,所以我可能不会一开始就把所有事情都做对。对此感到抱歉。此外,我对 Java 和一般
#include #include using namespace std; class Integer { public: int i; Integer (int ll
我不明白: ArrayList list = new ArrayList(); Collection list1 = new ArrayList(); 类 ArrayList扩展实现接口(interf
我编写了:。它成功了。我不知道为什么?
我编写了:。它成功了。我不知道为什么
我编写了:。它成功了。我不知道为什么?
Collectors.counting()返回 long此方法中每个键的值: private static Map countDuplicates(HashSet cards) { retur
我正在尝试通过搜索旧元素并将其替换为新元素来更新节点的元素。但是有一个我不明白的错误。是什么导致我的代码出现该错误,我该如何解决?错误; The method update(Integer, Inte
我有一个称为 client 的表,其中有一列称为created_time ,所以实际上我想绘制一个 map ,以便我可以知道在哪一年和哪一个月添加了多少客户?现在的要求是假设在 2018 年 11 月
这个问题已经有答案了: Is Java "pass-by-reference" or "pass-by-value"? (91 个回答) 已关闭 8 年前。 我对 ArrayList Collecti
我意识到下面的代码是正确的 Integer.MIN_VALUE == -Integer.MIN_VALUE == Math.abs(Integer.MIN_VALUE) 这是因为当我们取反-21474
我有以下类 AccountWebappGridRow,它扩展了 AccountGridRow: public class AccountWebappGridRow extends AccountGri
我正在学习 Haskell 并看到了函数组合。 尝试复合 map和 foldl mapd = (map.foldl) 比 test = (mapd (\x y -> x + y ) [1,2,3,4]
我有两个相同大小的数组和两个方法。 public class Client { private static int[] ints; private static final int
我喜欢 Java 8 中的 Streams 概念。现在我想借助 Java Streams 将 Java 中的 Map 转换为排序列表。我只想显示列表而不将其存储在任何地方。我希望在结果列表中有这个输出
我有一个数据库表,其中包含电视节目类型列表和关联的 ARGB 颜色值,用于在显示电视指南时突出显示 Android ListView 中的电视节目。流派表看起来像这样... id genre
我有一个 Integer 类,它应该模拟一个整数 mod n。因此,它具有如下构造函数: Integer::Integer(int x) : m(x), n(0) { } Integer::I
我是一名优秀的程序员,十分优秀!