具有相同拟合值的回归模型在 r 中使用相等运算符相互比较时输出 FALSE

Question

给定来自 faraway 包的 odor 数据集的二次模型：

> lmod <- lm(odor ~ I(temp) + I(gas) + I(pack)+I(temp^2)+I(gas^2)+I(pack^2)+I(temp*gas)+I(gas*pack)+I(pack*temp),odor)
> lmod6 <- lm(odor ~ polym(temp,gas,pack,degree = 2),odor)

两个模型具有相同的拟合值：

> lmod$fitted
        1         2         3         4         5         6         7         8         9        10        11        12 
 86.62500  45.87500  36.12500  28.37500  42.50000  15.25000  -3.25000 -24.50000  59.87500  29.37500  20.62500 -16.87500 
       13        14        15 
-30.66667 -30.66667 -30.66667 
> lmod6$fitted
        1         2         3         4         5         6         7         8         9        10        11        12 
 86.62500  45.87500  36.12500  28.37500  42.50000  15.25000  -3.25000 -24.50000  59.87500  29.37500  20.62500 -16.87500 
       13        14        15 
-30.66667 -30.66667 -30.66667

但是，当将这些拟合值相互比较时，它们并不相同，这是为什么？

> table(lmod6$fitted==lmod$fitted)

FALSE  TRUE 
   13     2

Answer 1

这是 R Inferno by Patric Burns or R FAQ

第 1 章中描述的浮点陷阱的经典示例

差异为零达到浮点精度：

options(digits=20)
lmod$fitted - lmod6$fitted
                        1                          2
 0.0000000000000000000e+00  0.0000000000000000000e+00
                         3                          4
 7.1054273576010018587e-15  1.0658141036401502788e-14
                         5                          6
-7.1054273576010018587e-15 -3.5527136788005009294e-15
                         7                          8
 4.4408920985006261617e-16 -3.5527136788005009294e-15
                         9                         10
-1.4210854715202003717e-14 -1.4210854715202003717e-14
                        11                         12
 7.1054273576010018587e-15  3.5527136788005009294e-15
                        13                         14
 7.1054273576010018587e-15  7.1054273576010018587e-15
                        15
 7.1054273576010018587e-15

all.equal() 函数用于测试两个向量是否为 "almost equal"（达到 fp 精度）：

all.equal(lmod$fitted,lmod6$fitted)
[1] TRUE

具有相同拟合值的回归模型在 r 中使用相等运算符相互比较时输出 FALSE

Regression models with same fitted values that output FALSE when compared to each other with equal operator in r

regression

r

equals