r - 两个理论上相同的代码不会产生相同的输出

Question

我正在尝试在 中运行例程电话 ，涉及复制随机实验。我遇到了输出维度的错误，这让我发现了以下特点:

length(replicate(100 - 100*8 / 10, 4))
# 20
length(replicate(100 *(1- 8/ 10), 4))
# 19

作为完整性检查，我运行了代码行以确保 replicate 的第一个参数中的表达式。产生相同的输出。

100 - 100*8 / 10
# 20
100 *(1- 8/ 10)
# 20

我想知道人们是否遇到同样的问题。我真正想知道的是，为什么会发生这种情况？

注意:我知道 rep 之间的区别和 replicate ，而我的例程需要后者，而不是前者。

Answer 1

不完全一样

(100 *(1- 8/ 10)) == 20
#[1] FALSE

(100 - 100*8 / 10) == 20
#[1] TRUE

因为

20 - (100 *(1- 8/ 10))
#[1] 3.552714e-15

和 n在 ?replicate是一个整数

n - integer: the number of replications.

将该输出转换为 integer floor秒到 19

as.integer((100 *(1- 8/ 10)))
#[1] 19

floor((100 *(1- 8/ 10)))
#[1] 19

一种选择是用 ceiling 包裹

length(replicate(ceiling(100 *(1- 8/ 10)), 4))
#[1] 20