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

【问题标题】：Fractional Factorial Design in R - replicate given matrixR中的部分因子设计 - 复制给定矩阵
【发布时间】：2017-10-12 06:57:28
【问题描述】：

我试图在 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 个变量不是全因子。

【问题讨论】：

标签： r experimental-design

【解决方案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)

还有更优雅的解决方案

【讨论】：