为什么 R unique() 不能处理从 lm() 对象中提取的拟合值？

Question

Unique 查找向量的唯一值。

如果我有一个数据框：

test_data <- data.frame(x = c(rep(1.00050239485720394857,4),
                              1.00050239485720394854,rep(2.0002230948570293845,5),rep(3.0005903847502398475,5)),
                        y = c(rep(4.00423409872345,5),rep(2.034532039485722,5),rep(1.1234152304957,5)))
sapply(test_data,unique)

R returns:

            x        y
[1,] 1.000502 4.004234
[2,] 2.000223 2.034532
[3,] 3.000590 1.123415

不出所料。

但是假设我拟合了一个 lm() 或 aov() 对象，然后尝试找到唯一的拟合值()：

set.seed(123)

y = rf(100,50,3,3)
x1 <- factor(c(rep("blue",25),
               rep("green",25),
               rep("orange",25),
               rep("purple",25)))

bsFit <- aov(y ~ x1)
unique(bsFit$fitted.values)

R returns:

 [1] 2.709076 2.709076 2.709076 2.709076 2.709076 2.709076
 [7] 2.709076 4.060080 4.060080 4.060080 4.060080 3.314801
[13] 3.314801 3.314801 3.314801 1.960280 1.960280 1.960280
[19] 1.960280 1.960280

此处明显重复

Answer 1

正如其他人所说（特别是@Tim-Biegeleisen），RStudio 正在将输出格式化为特定的小数位数（记住打印到控制台的任何内容都是由 RStudio 格式化的）。因此 "duplicates"，如果格式正确以显示所有小数位，则不会重复。

我们可以用format来显示所有小数位：

format(unique(bsFit$fitted.values), digit = 22)
 [1] "2.7090760788376542" "2.7090760788376773" "2.7090760788376604" "2.7090760788376622" "2.7090760788376627"
 [6] "2.7090760788376649" "2.7090760788376640" "4.0600797479202155" "4.0600797479202164" "4.0600797479202200"
[11] "4.0600797479202146" "3.3148005388803132" "3.3148005388803128" "3.3148005388803146" "3.3148005388803137"
[16] "1.9602804435309986" "1.9602804435309984" "1.9602804435309982" "1.9602804435309988" "1.9602804435310004

我在抛出错误之前尝试了位数并得到 22。

为什么 R unique() 不能处理从 lm() 对象中提取的拟合值？

Why doesn't R unique() work with fitted values extracted from lm() objects?

r

unique

anova