Дробно-факториальный дизайн в R - повторяет заданную матрицу

Question

Я пытаюсь воспроизвести следующий дробный факторный план в R. Каждая переменная имеет 3 уровня и 10 переменных. В этом проекте первые 3 переменные (AC) тестируются полностью факториально.

Я попытался:

library(AlgDesign)

levels.design = c(3,3,3,3,3,3,3,3,3,3)
f.design <- gen.factorial(levels.design)

fract.design <- optFederov(
  data=f.design,
  nTrials=sum(levels.design),
  approximate=FALSE)

А также:

    library(DoE.base)
     test.design <-oa.design(nlevels = c(3,3,3,3,3,3,3,3,3,3))

Это дает мне 3^10 дизайнов с 27 прогонами, однако первые 3 переменные не являются полными факториалами.

Answer 1

Я думаю, что мне удалось потенциальное решение:

куда

Шаг первый Я создаю полную факториальную матрицу 3^3

Шаг второй Я создаю ортогональную матрицу 3^7, фиксированную на 27 прогонах

Шаг третий объединить два.

Full_Factorial <- fac.design(nlevels=3, nfactors=3, randomize=TRUE)

Fractional_Factorial <- oa.design (ID=NULL, nruns=27, replications=1, repeat.only=FALSE,
  randomize = TRUE, seed = 27761, nlevels = c(3,3,3,3,3,3,3), factor.names = list(D=c(1,2,3),
   E=c(1,2,3), F=c(1,2,3), G=c(1,2,3), H=c(1,2,3), I=c(1,2,3), J=c(1,2,3)), columns ="order", 
    min.residual.df = 0)

Merged_matrix <- cbind(Full_Factorial, Fractional_Factorial)

Есть еще более элегантные решения