adding in whole plot column to the split plot design

biometryhub · Aug 1, 2024 · 5f714f4 · 5f714f4
1 parent bea7f55
commit 5f714f4
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/R/design.R b/R/design.R
@@ -156,6 +156,11 @@ design <- function(type,
                                              trt2 = sub_treatments,
                                              r = reps,
                                              seed = ifelse(is.numeric(seed), seed, 0))
+       numsp <- max(as.numeric(outdesign$book$splots))
+       lenblk <- as.vector(table(outdesign$book$block)[1])
+       numwp <- lenblk/numsp
+       outdesign$book$wplots <- rep(rep(1:numwp, each = numsp), reps)
+       outdesign$book <- outdesign$book[, c(1, 3, 6, 2, 4, 5)]
     }
 
     else if(substr(tolower(type), 1, 7) == "crossed") {