0

我试图在 R 中重现以下部分因子设计。每个变量有 3 个级别，并且有 10 个变量。在此设计中，前 3 个变量 (AC) 进行了全因子测试。

3 10 级因子部分因子设计

我尝试过：

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

4

1 回答 1

0

我想我已经管理了一个潜在的解决方案：

在哪里

第一步，我创建一个全阶乘 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)

有更多优雅的解决方案

于 2017-10-12T10:39:30.367 回答