我想生成一个包含三列(A,B和C)的大型数据框(100000行和3列).
该数据框满足以下两个条件:
>在每一行中,A B C = 1;
>所有A都有一个三角形分布(min = 0.2,mod = 0.3,max = 0.4),所有B都有一个三角形分布(min = 0.3,mod = 0.4,max = 0.5),所有C都有三角形分布(分= 0.1,MOD = 0.3,最大= 0.5).
我无法弄清楚如何生成这种数据集.
非常感谢你提前的建议.
Songchao
最佳答案
N <- 100000
library(triangle)
A <- rtriangle(N, 0.2, 0.4, 0.3)
B <- rtriangle(N, 0.3, 0.5, 0.4)
C <- 1 - A - B
d = data.frame(A, B, C)
summary(d)
稍后修改
nr <- 100000
u1 <- runif(nr)
u2 <- runif(nr)
u3 <- (2 - u1 - u2) / 2
U <- cbind(u1, u2, u3)
# shuffle, because I am not sure about the tails of u3
for (i in (1:nrow(U))) {
U[i, ] <- U[i, sample(1:3)]
}
t1 <- qtriangle(U[, 1], 0.2, 0.4, 0.3)
t2 <- qtriangle(U[, 2], 0.3, 0.5, 0.4)
t3 <- qtriangle(U[, 3], 0.1, 0.5, 0.3)
d <- cbind(t1, t2, t3)
summary(d)
cor(d)