在R中生成实验设计数据

Question

我有一个随机区组设计实验，在这个实验中我们想测试10个处理(8个基因型+2个对照)。

它的结构如下：

9个点，每个点内有4个区块(重复)，每个区块内有10个小区，每个小区内有144个特定处理(基因型)的个体。

我想生成一个数据集，如下所示：

注:对每个区块内的地块随机分配处理。

Site <- LETTERS[1:6]
Block <- LETTERS[1:4] ## For each Site
Plot <- paste(rep("P",10),seq(1,10,1),sep="_") ## For each Block 
Trt <- c(LETTERS[1:8],rep("control",2)) ## 10 Treatments within each block
id <- seq(1,144,1) ## for each Plot
observation <- runif(n=144,min=1, max=4) ## repeated for each treatment, each block,and each site.

data.frame(Site,Block,Plot,Trt,id,observation)

知道这段代码会丢弃一个错误，我希望为同一站点拥有所有的块级别，对于相同的块，所有的图级别及其相应的处理，对于相同的图，所有的id标签及其相应的观察结果。

非常感谢

Answer 1

我想你可能在找expand.grid：

Site <- LETTERS[1:6]
Block <- LETTERS[1:4] ## For each Site
Plot <- paste(rep("P",10),seq(1,10,1),sep="_") ## For each Block 
data <- expand.grid(Site = Site, Block = Block, Plot = Plot)
data$Trt <- sample(c(LETTERS[1:8],paste0("control",1:2)), nrow(data), replace = TRUE)
data