为何shape=1的rweibull与rexp采样结果不同？概率密度却等价

当Weibull分布的shape参数为1、scale参数为1时，R中rweibull()的采样结果和rate=1的rexp()采样结果不一致，但两者的概率密度函数是等价的，具体验证如下：

采样结果不一致验证 #

n=2
set.seed(1)
rweibull(n,shape = 1,scale = 1)
# [1] 1.3261078 0.9885284
set.seed(1)
rexp(n,rate = 1)
# [1] 0.7551818 1.1816428

概率密度函数等价验证 #

x <- c(0, rlnorm(50))
all.equal(dweibull(x, shape = 1), dexp(x))
# 返回TRUE，说明PDF完全一致

请问这是逆变换采样的问题吗？如果是，原因是什么？

这不是逆变换采样本身的问题，核心原因是R中rweibull()和rexp()的随机数实现逻辑存在差异：

1. 理论等价性确认
   当Weibull分布的shape=1时，其分布本质就是参数为rate=1/scale的指数分布。这里scale=1，对应指数分布rate=1，所以两者的概率密度函数完全一致，验证结果符合预期。

2. 采样实现的差异

• rweibull()采用Weibull分布通用的逆变换采样公式：scale * (-log(u))^(1/shape)，当shape=1时简化为-log(u)（u是(0,1)的均匀随机数）。
• rexp()则使用了专门针对指数分布优化的随机数生成算法，并非直接复用Weibull的采样逻辑。

3. 种子相同但结果不同的原因
   即使设置了相同的随机种子，两个函数调用底层随机数生成器的方式不同：rweibull()直接从均匀分布随机数序列中取值计算，而rexp()的实现可能对随机数做了额外处理，或者调用的随机数序列偏移不同，导致最终生成的采样值不同。

手动模拟验证这一点：

set.seed(1)
u <- runif(2)
# 模拟rweibull(shape=1, scale=1)的采样逻辑
weib_sim <- (-log(u))
# 直接调用rexp
set.seed(1)
exp_sim <- rexp(2, rate=1)

weib_sim
# [1] 1.3261078 0.9885284  # 和rweibull的输出完全一致
exp_sim
# [1] 0.7551818 1.1816428  # 和rexp的输出完全一致

总结 #

两个分布在理论上完全等价，概率密度函数一致是正确的；采样结果不同仅仅是R中两个函数的随机数实现细节差异导致，和逆变换采样的正确性无关。

内容的提问来源于stack exchange，提问作者Sle R.