Runif 未生成均匀分布

Question

set.seed(1432)    
n_len <- 400000
jdc<- data.frame(rnd = numeric(n_len))
jdc$rnd <- runif(n_len,0,1)
ggplot(jdc,aes(x = rnd)) + geom_density()

您会注意到随机变量的分布向两个边界下降。

我正在尝试基于一些小于 rnd 的 i 进行采样，但 i 的范围在 0, .05 之间，因此这种分布是一个问题。

Answer 1

与 runif() 相比，这与密度估计器的关系更大。最好用直方图看数据:

ggplot(jdc, aes(x = rnd)) + geom_histogram(binwidth = 0.01, boundary = 0)

作为joran已经指出，还可以创建一个直方图来显示与密度估计器类似的偏差:

ggplot(jdc, aes(x = rnd)) + geom_histogram()

直方图的优点是很容易理解，为什么会这样。最左边和最右边的 bin 分别以 0 和 1 为中心。这意味着，例如，最左边的 bin 从 -0.005 变为 0.005。但是没有低于零的数据点，所以这个区间包含的点数只有其他区间的一半左右。