R中的正射规划（分数因子设计）

我试图在R中为联合分析实验创建一个析因设计（如SPSS orthoplan）

在过去的Stackoverflow问题中，我找到了以前的答案：

这确实是一个有用的答案，但只有在你有数字水平的因素的情况下

不幸的是，这不是我的情况，因为我想要使用的因子是名义变量，即它们的级别不是数字类型而是因子类型：例如，我必须处理一个指示产品颜色的因子，它可以是绿色、黄色或红色

我试图修改作为问题答案的代码 以这种方式：

f.design <- gen.factorial(levels.design,factors="all")

但为了在具有标称系数的设计中获得平衡的解决方案，我希望它至少应该是：

nTrials=prod(unique(levels.design))

---

无论如何，有一个软件包可以处理这样的问题，它是Ulrike Groemping教授的软件包FrF2，但它只处理二分法变量，我不知道如何使用它来解决我的问题。

经过一段时间的研究，我可以在这里分享我的发现：

是的，您可以在R中构建正交设计，其方式与SPSS orthoplan中的方式类似

只需将变量

nlevels

定义为包含变量级别的向量

然后你必须打电话：

库（DoE.base）

fract.design
nTrials=prod(unique(levels.design))