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 - 重複給定矩陣中的分式因子設計

回答

相關問題