使用rep()函数生成向量长度小于times参数的问题求助

问题场景 #

循环中需按指定次数复制字符串，复制次数从seq(0.02,0.50,0.02)取值递增。当预期复制56次（对应x=0.14，y=400，x*y=56）时，rep()却生成了长度为55的向量，导致循环异常。但直接赋值x=0.14时，rep()能正常生成56长度的向量，两种场景下x均为numeric类型。将x转为字符再转回数值可解决该问题，以下是原因解析。

代码复现 #

# 生成56长度向量
x <- 0.14
y <- 400
test <- rep("ES", times = x*y)  
length(test) # 输出56
class(x)     # 输出"numeric"
  
# 生成55长度向量
x <- seq(0.02, 0.50, 0.02)[7]
y <- 400
test <- rep("ES", times = x*y)
length(test) # 输出55
class(x)     # 输出"numeric"

# 可行的解决方法
i <- seq(0.02, 0.50, 0.02)
x <- as.character(i[7])
x <- as.numeric(x)
y <- 400
test <- rep("ES", times = x*y)
length(test) # 输出56

原因解析 #

核心是浮点数的二进制精度误差：

• 直接赋值x=0.14时，R存储的是一个接近0.14的精确浮点数，计算x*y后刚好得到精确的56。
• 通过seq(0.02,0.50,0.02)[7]获取的x，是浮点数累加的结果。由于二进制无法精确表示部分十进制小数，这个x实际是略小于0.14的数值（用print(x, digits=20)可查看：0.13999999999999999左右）。
• 计算x*y后得到的是略小于56的数值（比如55.999999999999986），而rep()的times参数若为非整数，会自动向下取整（复制次数必须是正整数），因此最终生成了55个元素的向量。
• 将x转为字符再转回数值的操作，相当于对近似值做了四舍五入，得到了精确的0.14，因此x*y=56，rep()正常生成56长度的向量。

替代解决方法 #

除了转字符再转数值，还有更直接的方案：

• 对计算后的次数取整：rep("ES", times = round(x*y))，确保次数为精确整数。
• 生成序列时避免浮点数精度问题，比如用整数步长再转换：seq(2,50,2)/100，生成的数值更精确。
• 用print(x, digits=20)查看浮点数的真实值，提前排查精度问题。

内容的提问来源于stack exchange，提问作者xzaktly24