Question

我希望获得三个级别和样本量n = 15的治疗随机化。我被困在哪里

志愿者＆lt; - 1:15 set.seed（1）;样本（志愿者，大小= 5，替换= F）

我想要三个不同的小组，每个小组五个，但我是R的新手。

这是ANOVA的数据设置，而不是给出特定数据集的特定问题。另外我不知道它对set.seed的含义是什么

Answer 1

我认为你正在寻找类似的东西：

set.seed(1337)
# replace with you real participants ids
volunteers <- 1:15
# set the number of groups
number.of.groups <- 1:3
# set group size
group.size <- 5
# generate data frame with participant > group order
df <- data.frame(group=sort(rep(number.of.groups,group.size)),
                 participant=sample(volunteers,length(volunteers)))

# show your groups
df[which(df$group==1),]
# group participant
# 1     1           9
# 2     1           8
# 3     1           1
# 4     1           6
# 5     1           5
df[which(df$group==2),]
# group participant
# 6      2           4
# 7      2          15
# 8      2           3
# 9      2           2
# 10     2          13
df[which(df$group==3),]
# group participant
# 11     3          11
# 12     3          10
# 13     3          14
# 14     3          12
# 15     3           7

如果您希望能够复制样本，则只需要使用set.seed（），因为此方法会导致您始终绘制相同的“随机”样本。因此，set.seed（）更多用于测试，而不是真正的分析代码。你设定的种子无关紧要。如果你想复制，请确保始终设置相同的种子。

Answer 2

怎么样：

install.packages("randomizr")
library(randomizr)
Z <- complete_ra(15, num_arms = 3)
table(Z)

这给出了

> table(Z)
Z
T1 T2 T3 
 5  5  5

R studio的治疗随机化

2 个答案: