R - replicate与えられた行列の分数因数設計

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) 

そして：

私の最初の3つの変数は階乗いっぱいではありませんしかし、27回のランで3^10のデザインを与える

library(DoE.base) 
    test.design <-oa.design(nlevels = c(3,3,3,3,3,3,3,3,3,3)) 

を。

Answer 1

私は潜在的な解決策に管理したと思う：私は固定直交する3^7マトリックスを作成し、私は完全要因3^3マトリックス

ステップ2を作成

ステップ1

を27回

ステップ3で2つをマージします。

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よりエレガントなソリューション

があります。