随机向量组合

Question

如何在R中随机组合两个向量？例如,

individuals <- rep(1:24)
treatments <- c("T1", "T2", "T3") 

考虑到我希望每个治疗有3个个体，那么我如何能够随机地分配给一个个体(从个体1到24)，所以我最终得到了一个矩阵，其中包含了个体和治疗之间的组合。

Answer 1

这不是一个优雅的解决方案，但它将完成这项工作。我假设你希望每种治疗至少有3个人，因为一个完全平衡的设计会给你每种治疗8个人。

individuals <- rep(1:24)
treatments <- c("T1", "T2", "T3") 

while (length(treatments[treatments=="T1"])<=3 | length(treatments[treatments=="T2"])<=3 | length(treatments[treatments=="T3"])<=3) {
    treatments = sample(treatments, length(individuals), replace = TRUE)
    print("retry")
    flush.console()
}

df <- cbind(individuals, treatments)