Экспериментальный дизайн в R: сбалансированный неполный блочный дизайн

Я хотел бы создать матрицу для планирования эксперимента.

  • У меня 6 человек и 6 процедур, доставленных за 6 дней.
  • Я могу использовать 1 человека в неделю (всего 6 недель по 6 дней = 36 дней испытаний с 6 процедурами в день).
  • Каждый день я буду тестировать 6 процедур на человека в случайном порядке. Человек не должен начинать с того же лечения в любой день недели.

Например:

Monday: A C D B E F
Tuesday: C E A B F D
...
  • Если возможно, разные лица не должны иметь тот же порядок, что и другие.
  • Я хотел бы принять во внимание порядок обработки; то есть каждую пару обработок, следующих друг за другом, следует повторять настолько часто, насколько это возможно, чтобы избежать эффектов "переноса" (последовательности и порядка). Сейчас я работаю с R, но могу работать и с MATLAB.

Я пытался разработать дизайн с пакетом "Agricolae" и пакетом "Crossdes", но пока безуспешно. Моими идеями было использование сбалансированного незавершенного блочного дизайна. Ниже приведен пример одного из моих кодов.

> trt <- c("A", "B", "C", "D", "E", "F" )
> k <- 6
> outdesign<- design.bib (trt, k, seed=0)
Error in md[i, bi] : incorrect number of dimensions
Источник

1 ответ

Вы используете неправильную функцию. Вы хотите законченный дизайн блока (design.rcbd), а не неполный. Неполная конструкция блока создаст блоки, в которые включены только некоторые из обработок. Следовательно, для design.bib, k

РЕДАКТИРОВАТЬ: Это даст вам 6 предметов на 6 дней. Принуждение каждого дня начинать с другого лечения для каждого предмета является более трудоемким, чем у меня сейчас есть время, но, возможно, кто-то еще ответит вам, прежде чем я займусь этим.

trt <- c("A", "B", "C", "D", "E", "F" )
N <- 6
days <- 6

outdesign<-as.data.frame(matrix(nrow=N*length(trt),ncol=days+1))
names(outdesign)<-c("Subj", "Day1", "Day2", "Day3", "Day4", "Day5", "Day6")
outdesign$Subj<-rep(1:6,each=length(trt))
for (i in 1:days) {
  outdesign[,i+1]<- design.rcbd(trt, N, seed=0,first=TRUE)$book$trt
}

   Subj Day1 Day2 Day3 Day4 Day5 Day6
1     1    E    F    D    C    C    C
2     1    C    E    E    F    F    E
3     1    D    C    B    A    E    D
4     1    B    D    F    D    A    A
5     1    A    B    C    B    D    F
6     1    F    A    A    E    B    B
7     2    D    C    F    A    F    F
8     2    B    E    B    C    E    B
9     2    C    D    E    E    D    A
10    2    A    F    A    F    B    E
11    2    F    A    D    B    C    D
12    2    E    B    C    D    A    C
13    3    B    C    C    A    B    F
14    3    A    A    E    F    E    E
15    3    C    E    F    B    A    B
16    3    E    B    A    C    F    A
17    3    F    F    B    D    C    C
18    3    D    D    D    E    D    D
19    4    F    C    F    D    F    E
20    4    D    B    A    E    B    D
21    4    B    F    D    A    A    F
22    4    C    A    B    B    D    B
23    4    E    E    E    C    C    C
24    4    A    D    C    F    E    A
25    5    D    E    F    A    A    F
26    5    A    D    E    E    E    B
27    5    C    C    C    D    D    E
28    5    E    F    B    C    C    C
29    5    B    A    A    F    B    A
30    5    F    B    D    B    F    D
31    6    F    F    A    F    C    F
32    6    A    D    E    E    D    B
33    6    C    B    C    D    A    A
34    6    B    C    F    C    F    C
35    6    E    E    B    A    E    D
36    6    D    A    D    B    B    E
Источник
Другие вопросы по тегам