R - 复制给定矩阵的分数因子设计

Question

我试图在R中重现以下小数阶乘设计。每个变量有3个级别，有10个变量。在该设计中，前3个变量（A-C）被测试为全因子。

我试过了：

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矩阵

第二步我创建一个固定在27次运行的正交3 ^ 7矩阵

第三步合并两者。

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)

有更优雅的解决方案