R - 重复给定矩阵中的分式因子设计

我想在R中重现以下分式因子设计。每个变量都有3个等级，并且有10个变量。在这个设计中，前三个变量（A-C）被测试为全因子。R - 重复给定矩阵中的分式因子设计

$3 Level 10 Factor fractional factorial design$

我曾尝试：

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

来源

2017-10-12 Adam

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

凡

第一步我创建了一个全因子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)

有likley更优雅的解决方案

来源

2017-10-12 10:39:30 Adam

R - 重复给定矩阵中的分式因子设计

回答

相关问题