R中的runif（或者我）有什么问题？

Question

我试图从均匀分布生成一个随机变量，也取决于随机变量。

理想情况下，我希望z2也是一个随机变量。

例如，

set.seed(123)
#generate first random variable z1
z1 <- runif(10,0,1)
#generate secondrandom variable z2
z2 <- ifelse(z1>= 0.7, runif(1,0.6,0.9), runif(1,0.2,0.5))
df <- data.frame(z1,z2)

但如果你看看z2它只包含两个级别。

          z1        z2
1  0.2875775 0.3360002
2  0.7883051 0.8870500
3  0.4089769 0.3360002
4  0.8830174 0.8870500
5  0.9404673 0.8870500
6  0.0455565 0.3360002
7  0.5281055 0.3360002
8  0.8924190 0.8870500
9  0.5514350 0.3360002
10 0.4566147 0.3360002

我们可以生成1000个数字并总结结果。

set.seed(123)
#generate first random variable z1
z1 <- runif(1000,0,1)
#generate secondrandom variable z2
z2 <- ifelse(z1>= 0.7, runif(1,0.6,0.9), runif(1,0.2,0.5))

table(df$z2)

0.352099318057299 0.655542095610872 
              707               293

Answer 1

您使用矢量化方法计算z2。因此，当您使用runif(1 ,0.6,0.9)时，实际上只拉动一个数字。但是，如果你要做runif(length(z1) ,0.6,0.9)，你将拉长度z1。

set.seed(123)
#generate first random variable z1
z1 <- runif(1000,0,1)
#generate second random variable z2
z2 <- ifelse(z1 >= 0.7, runif(length(z1) ,0.6,0.9), runif(length(z1),0.2,0.5))

df <- data.frame(z1,z2)
table(df$z2)

Answer 2

你的runif(1,0.6,0.9)和runif(1,0.2,0.5)函数只运行过一次。每次将z1与0.7进行比较时，它们都不会运行。因此，z2最终只有两个值。

您可以考虑使用replicate每次再次运行runif函数。更新：或像@DJV提议的length(z1)。